el de arriba es newton rhapson con do y este con do while

program newton_r
implicit none
real::xo,tol,E,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
E=1
do while(I<=max.and.E>tol)
x=xo-(f(xo)/df(xo))
E=abs(x-xo)
xo=x
I=I+1
end do
if(E<=tol)then
write(*,*)x
else
write(*,*)'no converge'
end if
end program