Su-34

Надо описать алгоритм решения нелинейного уравнения, решенного методом Ньютона в Фортране:
Помогите если несложно.

program newton
eps=0.001
x=0
do
dx=f(x)/df(x)
if(abs(dx)<eps) exit
x=x-dx
print*,x
end do
print*,'x=',x,'f(x)=',f(x)
end
function f(x)
f=3*x-cos(x)-1 ----------функция исходная
end
function df(x)
df=3+sin(x)
end