d1=5.6 rr1=qext(d1) print *,rr1 stop end