programa q calcula el valor de una ecuciaon de newton raphson

program n_rdo
implicit none
real::xo,tol,x,df,f
integer::max,I
f(xo)=xo**2+10*xo-1
df(xo)=2*xo+10
write(*,*)'dame el valor de xo,tol,max'
read(*,*)xo
read(*,*)tol
read(*,*)max
I=1
do I=1,max,1
x=xo-(f(xo)/df(xo))
xo=x
if(abs(x-xo)<=tol)then
write(*,*)x
stop
else
xo=x
end if
end do
end program