用牛顿法解方程组

				function X=Newton(X0,tol,nmax,varargin)%输入初值X=X0+1，误差tol，和最大迭代次数，X0为行向量格式。
				X=X0'+1;%                               可以编辑F.m和J.m文件来输入方程组和对应的雅克比矩阵。
				niter=0;
				while ((X(1)-X0(1)) >= tol || (X(2)-X0(2)) >= tol || (X(3)-X0(3)) >= tol ||....
				        (X(4)-X0(4)) >= tol) && niter 				    niter=niter+1;
				    X1=X-inv(J(X))*F(X);
				    X0=X;
				    X=X1;
				end
				niter
				end
				%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
				%Newton([0 0 0 0],1.e-3,100)
				%niter = 20
				%ans =
				%    2.9937
				%    0.0317
				%   -0.3240
				%    0.2607