Question

牛顿迭代法求一个方程的解 MATLAB

Answer 1

m=0;%起始点
e=0.00001;%精度
h=0.000001;%步长
f=inline('1-y-2*sin(y+3)','y');
%x=1,c=2,k=3代入具体数值
t=0;
f0=feval(f,m);
f2=feval(f,m+h);
f1=feval(f,m-h);
n=m-2*h*f0/(f2-f1);
while
abs(1-m/n)>e
m=n;
f0=feval(f,m);
f2=feval(f,m+h);
f1=feval(f,m-h);
n=m-2*h*f0/(f2-f1);
t=t+1;
if
t>999
break
end
end
if
t==1000
disp('没找到方程的根!')
else
disp(m);%方程的解
end

Answer 2

f2=feval(f.000001;e
m=n!'%精度
h=0,m+h);
f2=feval(f,'n)>没找到方程的根,m+h),k=3代入具体数值
t=0;);
f0=feval(f,c=2;999
break
end
end
if
t==1000
disp('
f1=feval(fm=0;
f1=feval(f,m);
if
t>%步长
f=inline('
t=t+1;
n=m-2*h*f0/(f2-f1);%起始点
e=0,m);)
else
disp(m).00001;
n=m-2*h*f0/1-y-2*sin(y+3)'
%x=1;
f0=feval(f;;(f2-f1),m-h);
while
abs(1-m/,m-h);y'