DO I=1600000


If (y (2). Lt. 0) then
B=- RM (J) * Y (2) * * P/(2.0 * T (1))
F=1 - EXP (B)
WRITE (*, *) f, Y (2)
Y (2)=SQRT (2.0 * DLOG (F) * T (1)/RM (J))


endif