alpha = 0.05;
NM = @(x,f)x-alpha*prod(([1,0]*f([x,1;0,x])).^[1,-1]);
f = @(x)32*x.^6-48*x.^4+18*x.^2-x.^0;
x = 0.5;
eps = 1;
while eps > 1e-5
    xn = NM(x,f);
    eps = abs(xn-x)
    x = xn;
end