数值分析实验一：牛顿迭代法matlab实现

特普

于 2022-11-12 20:01:22 发布

阅读量788

点赞数 1

文章标签： matlab

本文链接：https://blog.csdn.net/qq_57945409/article/details/127824579

版权

1.求，精度为10-6

clc,clear
f=inline('x^2-113');
df=inline('2*x');
[y,k]=newton(f,df,10,10^(-6),500,1)
function[y,k]=newton(fname,dfname,x0,tol,N,m)
y=x0;
x0=y+2*tol;
k=0;
while abs(x0-y)>tol&k<N
    k=k+1
    x0=y
    y=x0-m*fname(x0)/dfname(x0);
    disp(y)
end
if k==N
    disp(warning)
end
end

clc,clear
f=inline('2*x^3-3*x^2+1');
df=inline('6*x^2-6*x');
[y,k]=newton(f,df,-0.4,10^(-3),500,1)
function[y,k]=newton(fname,dfname,x0,tol,N,m)
y=x0;
x0=y+2*tol;
k=0;
while abs(x0-y)>tol&k<N
    k=k+1
    x0=y
    y=x0-m*fname(x0)/dfname(x0);
    disp(y)
end
if k==N
    disp(warning)
end
end