g0 has to be initialized

matlab/numgrad2.m

@@ -32,6 +32,7 @@ g=zeros(n,1);
 badg=0;
 goog=1;
 scale=1;
+g0 = 0;
 for i=1:n
     if size(x,1)>size(x,2)
         tvecv=tvec(i,:);
@@ -57,4 +58,4 @@ for i=1:n
         g(i) = 0;
         badg = 1;
     end
-end
+end