%牛頓法 
%測試函數 fun(x)= x^5+2*x^3-x^2-3*x^2+x-1;
%定義變量
syms x;
%賦初值
a0=input('請輸入初始值 a0=');
esp=input('請輸入精度 esp=');
fun=input('請輸入函數 fun(x)=');
%轉化為多項式函數
d1=inline(diff(fun));
d2=inline(diff(fun,2));

a1=a0-d1(a0)/d2(a0);
flag=0;

while abs(a1-a0)>esp
    a0=a1;
    a1=a0-d1(a0)/d2(a0);
    flag=flag+1;
end

disp('迭代次數為:'),disp(flag)
disp('在x=:'),disp(a0)
disp('處取到最優解')
d3=inline(fun);
y=d3(a1);
disp(y);