?? newdon.m
字號:
%牛頓法
%測試函數 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);
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -