nevot_n.m

来自「数值分析各种公式的解法」· M 代码 · 共 22 行

function x=nevot_n(f_name x0)
x=x0;
xb=x_999
n=0;
del_x=0.01
while abs(x-xb)>0.000001
    n=n+1;
    xb=x;
    if n>300 break;
    end
    y=feval(f_name,x)
    y_driv=(feval(f_name,x+del_x)-y)/del_x
    x=xb-y/y_driv;
    fprintf('n=%3.of,x=%12.5e,y=%12.5e',n,x,y)
    fprintf('yd=%12.5e\n',y_driv)
end
fprintf('\n'final answer=%12.6e\n',x)
eqn_1.m
function y=eqn-1(x)
y=(0.01*x+1).*sin(x)-(x-0.01).*(x^2+1).^(-1)-0.0096;
f_name=eqn.1