%用牛頓法求方程的根
%x(k+1)=x(k)-f(xk)/df(xk)
%不是特別好用，有時不收斂

x1=0;
e=0.0005;

for i=1:200
    y=f(x1);
    x2=x1-y/df(x1);
    dy=df(x2);
    if abs(dy)<e
        break;
    else
        x1=x2;
    end
end

i
x2
x2-x1
f(x2)