求 牛顿切线法解方程组的matlab程序

来源:百度知道 编辑:UC知道 时间:2024/06/08 10:50:01
我们写论文报告,对于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;