%%%程序编写者 西北工业大学自动化学院 Email: yincwxa2013@mail.nwpu.edu.cn
%% All rights reserved
clear
clc
x1=input('输入迭代初值x1=')
x2=input('输入迭代初值x2=')
t=1;
epuc=input('输入精度要求eupc=')
if epuc>0
epuc=epuc;
else epuc=0.000001
end
while(t>epuc)
x3=(1+x2-0.1*exp(x1))/4
x4=(x1-x1^2/8)/4
t=max(abs([x3-x1,x4-x2]));
x1=x3;
x2=x4;
end
x=[x3 x4]
%%%%%
输入迭代初值x1=0.423
x1 =
0.4230
输入迭代初值x2=2.58
x2 =
2.5800
输入精度要求eupc=0.001
epuc =
1.0000e-003
x3 =
0.8568
x4 =
0.1002
x3 =
0.2161
x4 =
0.1913
x3 =
0.2668
x4 =
0.0526
x3 =
0.2305
x4 =
0.0645
x3 =
0.2346
x4 =
0.0560
x3 =
0.2324
x4 =
0.0569
x3 =
0.2327
x4 =
0.0564
x =
0.2327 0.0564
转载本文请联系原作者获取授权,同时请注明本文来自殷春武科学网博客。
收藏
分享
分享到: