      REAL(8) :: x,t,x0=0,t0=1,te=11,h
      REAL(8) :: rk1,f
      INTEGER :: n=10000,i

      h=(te-t0)/n

      t=t0
      x=x0
      DO i=1,n
         rk1=h*f(x,t)
         x=x+rk1
         t=t+h
      END DO

      PRINT *,'t=',t,'  x=',x,'  60-x=',60-x

      END
!--------------------------------------
      FUNCTION f(x,t)

      REAL(8) :: x,t,f

      f=(x+SQRT(x*x+t*t))/t

      END



      
      
