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

[f(0), f(1)];

[is(f(0)), is(f(1)), is(f(2)),is(f(-1/2))];

[f(2),f(-1/2)];

[f(0),f(1),f(2),f(-1/2)];

load("mnewton")$

first(mnewton(f(x),x,2));

f(1.475686517795721);

float(%pi/2);