kill[all];

 /* Para empezar, eliminamos la variable secante de la lista de variables */
 kill(secante);
 /* Creamos la función de iteración (secante(arg)) */
 secante(arg):=block([a:arg[1],b:arg[2],f:arg[3],xsec], xsec:(a*ev(f,x=b)-b*ev(f,x=a))/(ev(f,x=b)-ev(f,x=a)), [b,xsec,f] ); 
/* La primera iteración en a = 1, b = 2 y f = x^2 - 2 */
ini:[1,2,x^2-2]; 
 /* Repetimos el proceso desde 1 hasta 10 iteraciones */
for i:1 thru 10 do (print(i,float(ini),float(ini:secante(ini)))); 
/* Dibujamos la función (Paso opcional)*/
plot2d([x^2-2],[x,-2,2]);