y(x):=x^(x-1)^2=2*x+1;

load(mnewton)$

a:second(first(first(mnewton(y(x),x,2))));

b:first(y(a));

[a,b];