%%%程序編寫者 ?西北工業大學自動化學院 ? ?Email： yincwxa2013@mail.nwpu.edu.cn

%% ?All rights reserved

clear

clc

x0=input('輸入迭代初值x0=')

syms x

f=input('輸入函數f(x)=')

df=diff(f)

t=1;

epuc=input('輸入精度要求eupc=')

if epuc>0

epuc=epuc;

else epuc=0.000001

end

%牛頓迭代法

while(t>epuc)

x=x0

x1=x-eval(f)/eval(df);

t=abs(x-x1);

x0=x1;

end

x1

%%%%%%%%%

輸入迭代初值x0=1.8

x0 =

1.8000

輸入函數f(x)=x^4-x^2-3

f =

x^4-x^2-3

df =

4*x^3-2*x

輸入精度要求eupc=0.00000000001

epuc =

1.0000e-011

x =

1.8000

x =

1.5842

x =

1.5223

x =

1.5175

x =

1.5175

x =

1.5175

x1 =

1.5175

轉載本文請聯系原作者獲取授權，同時請注明本文來自殷春武科學網博客。

鏈接地址：http://blog.sciencenet.cn/blog-748007-751141.html

上一篇：數值分析--非線性方程組牛頓迭代法matlab程序

下一篇：數值分析--非線性函數拋物線迭代法matlab程序

總結

以上是生活随笔為你收集整理的数值分析牛顿法 matlab,科学网—数值分析--非线性函数牛顿迭代法matlab程序 - 殷春武的博文...的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。