You are an expert in writing code in many different languages. Your goal is to perform code completion for the following code, keeping in mind the rest of the code and the file meta data.
    Metadata details:
    filename: Newtons.m
    path of file: ./repos/Algorithms_MathModels/《基于MATLAB的高等数学问题求解》 随书附带源程序/CH10
    
    the code of the file until where you have to start completion: function [x,fval,iter,exitflag]=Newtons(fun,x0,eps,maxiter)
%NEWTONS   牛顿法求非线性方程组的根
% X=NEWTONS(FUN,X0)  牛顿法求非线性方程组的解，初始迭代点为X0
% X=NEWTONS(FUN,X0,EPS)  牛顿法求非线性方程组的解，精度要求为EPS
% X=NEWTONS(FUN,X0,EPS,MAXITER)  牛顿法求非线性方程组的解，最大迭代次数为MAXITER
% [X,FVAL]=NEWTONS(...)  牛顿法求非线性方程组的解并返回解处的函数值
% [X,FVAL,ITER]=NEWTONS(...)  牛顿法求非线性方程组的解并返回迭代次数
% [X,FVAL,ITER,EXITFLAG]=NEWTONS(...)  牛顿法求非线性方程组的解并返回迭代成功标志
%
% 输入参数：
%     ---FUN：非线性方程组的符号表达式
%     ---X0：初始迭代点向量
%     ---EPS：精度要求，默认值为1e-6
%     ---MAXITER：最大迭代次数，默认值为1e4
% 输出参数：

    
function [x,fval,iter,exitflag]=Newtons(fun,x0,eps,maxiter)
%NEWTONS   牛顿法求非线性方程组的根
% X=NEWTONS(FUN,X0)  牛顿法求非线性方程组的解，初始迭代点为X0
% X=NEWTONS(FUN,X0,EPS)  牛顿法求非线性方程组的解，精度要求为EPS
% X=NEWTONS(FUN,X0,EPS,MAXITER)  牛顿法求非线性方程组的解，最大迭代次数为MAXITER
% [X,FVAL]=NEWTONS(...)  牛顿法求非线性方程组的解并返回解处的函数值
% [X,FVAL,ITER]=NEWTONS(...)  牛顿法求非线性方程组的解并返回迭代次数
% [X,FVAL,ITER,EXITFLAG]=NEWTONS(...)  牛顿法求非线性方程组的解并返回迭代成功标志
%
% 输入参数：
%     ---FUN：非线性方程组的符号表达式
%     ---X0：初始迭代点向量
%     ---EPS：精度要求，默认值为1e-6
%     ---MAXITER：最大迭代次数，默认值为1e4
% 输出参数：
%     ---X：非线性方程的近似解向量
%     ---FVAL：解处的函数值
%     ---ITER：迭代次数
%     ---EXITFLAG：迭代成功标志，1表示成功，0表示失败
%
% See also newton

if nargin<2
    error('输入参数至少需要2个.')
end