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