      isd=124381
      do 10 i=1,15
      isd=isd*1566083941+1
      p=float(isd)*2.3283064e-10+0.5
      write(6,1) p
    1 format(1x,f12.8)
   10 continue
      stop
      end
