m = 3 n = 2 h = .1 d = 1 AA = {}; FV = []; GGAM = {}; BGAM = {}; PHI = {}; hh = []; dd = []; global GG for t=1:10000 %h = rand; %d = rand; %[A,phi,fval] = gamcheck01(m,n,h,d); [A,gam,fval] = agcheck(m,n,h,d); %[gam,fval] = gamrmlp(A,d,n); [A,fval] = gamrmlp2(gam,h); Bound = d*sum(h.^(0:n)); if abs(fval) > Bound + 1e-10 BGAM{end+1} = gam; else GGAM{end+1} = gam; end save GOODBADGAM GGAM BGAM end