      do 10 i=1,10000000
      isd=random()
      p=FLOAT(isd)*4.65566128e-10
   10 continue
      stop
      end
