y(x):=x^sqrt(x)=4;

load(mnewton)$

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

b:first(y(a));

[a,b];

integrate(sec(x)*tan(x)^2,x) + C;