newton := proc (f, x0, s) local x; x := x0; while abs(f(x)) > 10^(-s) do x := x-f(x)/(D(f)(x)); print(evalf(x,s)); end do; end proc;

Zio2JUkiZkc2IkkjeDBHRiVJInNHRiU2I0kieEdGJUYlRiVDJD44JDklPyhGJSIiIkYvRiUyKSIjNSwkOSYhIiItSSRhYnNHJSpwcm90ZWN0ZWRHNiMtOSQ2I0YsQyQ+RiwsJkYsRi8qJkY6Ri8tLUkiREdGJTYjRjtGPEY1RjUtSSZwcmludEdGODYjLUkmZXZhbGZHRjg2JEYsRjRGJUYlRiU=

f:=x->(x-1)*(x+3);

Zio2I0kieEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiYsJjkkIiIiISIiRixGLCwmRitGLCIiJEYsRixGJUYlRiU=

newton(f,5,50);

JCJTTExMTExMTExMTExMTExMTExMTExMTExMQiEjXA==

JCJTbm1tbW1tbW1tbW1tbW1tbW1tbW1tbW1tNyEjXA==

JCJTI1IhKTRYRidvOiNSISk0WEYnbzojUiEpNFhGJ286NSEjXA==

JCJTbWlUIz5hcipRYygpKjNZUHFvJnknZSgzTzUxKzUhI1w=

JCJTOEJmLVJaWil6ZXMsJio9QiRbZEE4JDQrKysrIiEjXA==

JCJTQW8uYkM4bykzNShcTS8lb0ArKysrKysrKysiISNc

JCJTIzNOJVx2NisrKysrKysrKysrKysrKysrKzUhI1w=

JCJTKysrKysrKysrKysrKysrKysrKysrKysrNSEjXA==

g:=x->(x-1)^2;

Zio2I0kieEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiQsJjkkIiIiISIiRiwiIiNGJUYlRiU=

newton(g,5,10);

JCIiJCIiIQ==

JCIiIyIiIQ==

JCIrKysrKzohIio=

JCIrKysrXTchIio=

JCIrKysrRDYhIio=

JCIrKytdaTUhIio=

JCIrKytESjUhIio=

JCIrK11pOjUhIio=

JCIrK0QieSsiISIq

JCIrXWkhUisiISIq

JCIrREomPisiISIq

JCIraWwoNCsiISIq

JCIrIkcpWys1ISIq

JCIrVFRDKzUhIio=

JCIrcT83KzUhIio=

JCIrTjUxKzUhIio=

JCIrPTAuKzUhIio=

JCIrZl8sKzUhIio=

JCIrSHcrKzUhIio=

save newton,"newton.mpl"; read "newton.mpl";

Zio2JUkiZkc2IkkjeDBHRiVJInNHRiU2I0kieEdGJUYlRiVDJD44JDklPyhGJSIiIkYvRiUyKSIjNSwkOSYhIiItSSRhYnNHJSpwcm90ZWN0ZWRHNiMtOSQ2I0YsQyQ+RiwsJkYsRi8qJkY6Ri8tLUkiREc2JEY4SShfc3lzbGliR0YlNiNGO0Y8RjVGNS1JJnByaW50R0Y4NiMtSSZldmFsZkdGODYkRixGNEYlRiVGJQ==