function?[xminf]?=?minGN（fx0vareps）
format?long;
if?nargin?==?3
????eps?=?1.0e-6;
end
S?=?transpose（f）*f;
k?=?length（f）;
n?=?length（x0）;
x0?=?transpose（x0）;
tol?=?1;
A?=?jacobian（fvar）;

while?tol>eps
????Fx?=?zeros（k1）;
????for?i=1:k
?????