Question

切线法求解 matlab实现

用切线法求下列方程的近似数值解y=x^4-3x^3+5cosx+8

Answer 1

牛顿切线法的MATLAB主程序：

function [k,xk,yk,piancha,xdpiancha]=newtonqx(x0,tol,ftol,gxmax)
x(1)=x0; 
for i=1: gxmax
   x(i+1)=x(i)-fnq(x(i))/(dfnq(x(i))+eps); piancha=abs(x(i+1)-x(i));
   xdpiancha= piancha/( abs(x(i+1))+eps); i=i+1;
xk=x(i);yk=fnq(x(i)); [(i-1) xk yk piancha xdpiancha];
if (abs(yk)<ftol)&&((piancha<tol)||(xdpiancha< tol))
     k=i-1;   xk=x(i);[(i-1) xk yk piancha xdpiancha];
     return;
end
end
 if i>gxmax
    disp('请注意：迭代次数超过给定的最大值gxmax。')
    k=i-1;   xk=x(i);[(i-1) xk yk piancha xdpiancha];
    return;
end
 [(i-1),xk,yk,piancha,xdpiancha]';

里面有可能符号是中文下的，朋友注意一下就OK了！