p=4.5 grid lam=0.05: tried 343; #(N(G)>N*)=184; #(Phi(G)<Phi*)=336; max N(G)-N*=1.539e-05 at (array([0.658, 2.818]), array([0.95, 0.05]))
p=4.5 grid lam=0.3: tried 343; #(N(G)>N*)=1; #(Phi(G)<Phi*)=342; max N(G)-N*=9.610e-16 at (array([-1., -1.]), array([0.5, 0.5]))
p=4.5 grid lam=1.0: tried 343; #(N(G)>N*)=7; #(Phi(G)<Phi*)=336; max N(G)-N*=5.201e-14 at (array([-1., -1.]), array([0.5, 0.5]))
p=4.5 grid lam=3.0: tried 343; #(N(G)>N*)=7; #(Phi(G)<Phi*)=336; max N(G)-N*=2.513e-13 at (array([-1., -1.]), array([0.05, 0.95]))
p=4.5 grid lam=10.0: tried 343; #(N(G)>N*)=7; #(Phi(G)<Phi*)=336; max N(G)-N*=1.669e-12 at (array([-1., -1.]), array([0.5, 0.5]))
