double precision u,v,x,y u=2.2 v=7.4d16 x=cot(u) y=cot(v) print *,x,y end