m = 3 n = 2 h = .1 d = 1 AA = {}; FV = []; GAM = {}; PHI = {}; hh = []; dd = []; global GG for t=1:100000 %h = rand; %d = rand; [A,phi,fval] = gamcheck01(m,n,h,d); Bound = d*sum(h.^(0:n-1)); %if abs(fval-Bound)<1e-9 if 1 AA{end+1} = A; PHI{end+1} = phi; FV(end+1)=fval; GAM{end+1}=GG; %hh(end+1)=h; %dd(end+1)=d; else if abs(fval) > Bound + 1e-10 'baaaaaaaaaaad' keyboard end end save GAMDAT AA PHI FV GAM sofar = length(GAM) end