      isd=124381
      do 10 i=1,10000000
      isd=isd*1566083941+1
      p=float(isd)*2.3283064e-10+0.5
   10 continue
      stop
      end
