f(x):=x^x^2-9;

 newton(t):=subst(t,x,x-f(x)/diff(f(x),x))$ newton(t);

a:newton(1.875);

for i : 1 thru 3 do (a:newton(a), print("        a:", a));