load(mnewton);

mnewton([x**2-8*x+15],[x],[1]);

mnewton([x**2-8*x+15],[x],[10]);