kill[all];

/* Definimos en la variable g la función a resolver */
g:x^2-2;
/* Definimos en la variable gd la derivada de la función */
gd:diff(g,x);
 /* Para empezar, eliminamos las asignaciones a la variable newton */
 kill(newton);
 /* Creamos la función de iteración (newton(pn)) con x = pn */
 newton(pn):=ev(x-g/gd,x=pn);
/* La primera iteración en pn = 0.2 */
pn:.2; 
 /* Repetimos el proceso desde 1 hasta 10 iteraciones */
i:0$
while i<10 do (
    i:i+1,
    print(i,float(pn)),
    pn:newton(pn)
);
/* Dibujamos la función (Paso opcional)*/
plot2d([x^2-2],[x,-2,2]);