share/numeric/newton1.mac

   1
   2 newton(exp,var,x0,eps):=
   3         block([xn,s,numer],
   4         numer:true,
   5         s:diff(exp,var),
   6         xn:x0,
   7   loop, if abs(subst(xn,var,exp))<eps then return(xn),
   8         xn:xn-subst(xn,var,exp)/subst(xn,var,s),
   9         go(loop) )$