matlab 高斯牛顿迭代求最小二乘问题

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