matlab数值分析程序--高等数学,数值代数的matlab实现-文字版, matlab电子书, 和matlab 有关的电子书:

2.7.2 割线法的MATLAB程序

2.7.2  割线法的MATLAB程序

割线法的MATLAB主程序

现提供名为gexian.m的M文件:

function [k,piancha,xdpiancha,xk,yk]=gexian (x01,x02,tol,ftol,gxmax)

x(1)=x01;x(2)=x02;

for i=2: gxmax

   u(i)= fnq(x(i))*(x(i)-x(i-1)); v(i)= fnq(x(i))-fnq(x(i-1));

   x(i+1)=x(i)- u(i)/( v(i)); 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-2) piancha xdpiancha xk yk]

 if (abs(yk)<ftol)&(( piancha <tol)|(xdpiancha< tol))

     k=i-2;  xk=x(i);yk=fnq(x(i));

 [(i-2) piancha xdpiancha xk yk];

 return;

end

end

 if i>gxmax

disp('请注意:迭代次数超过给定的最大值gxmax.')

k=i-2;  xk=x(i);yk=fnq(x(i));

    return;

end

欢迎转载,转载请注明来自一手册:http://yishouce.com/book/3/3020702.html
友情链接It题库(ittiku.com)| 版权归yishouce.com所有| 友链等可联系 admin#yishouce.com|粤ICP备16001685号-1