kill[all];
/* guardamos en la variable g la función a iterar */
g:(3*x-1)^(1/3);
/* Para empezar, eliminamos la variable iter de la lista */
kill(iter);
/* Creamos la función de iteración con x = pn */
iter(pn):=ev(g,x=pn);
/* La primera iteración en pn = -3 */
pn:-3; 
/* Repetimos el proceso desde 1 hasta 15 iteraciones */
for i:1 thru 15 do (print(i,float(pn)),pn:iter(pn));
/* Dibujamos la función a iterar incluída y=x */
plot2d([x,(3*x-1)^(1/3)],[x,-3,3]);