real(8) function f1(x, h)
! met. Newton
! Derivata calculata numeric
! "f(x)" definita in "fun.f90"

!real h/1E-3/
implicit real(8) (a-h, o-z)

f1 =(f(x+h) -f(x))/h

END
