% Code revision 22 - Better bisection algorithm -- fuzzy logic: % Only change mus on every second iteration; if target Ns are outside the range then extend range of mu slightly instead of halving it. >> Master(4,8,-1,0.001,20,1,50); Performing initial iterations Beginning bisection of mus to find Ns mu_up = 0.000000, mu_down = 0.000000, N_up = 0.000006, N_down = 0.000006, Delta_Int_error = 0.20%, N error 420.00% mu_up = 0.000000, mu_down = 0.000000, N_up = 0.000008, N_down = 0.000008, Delta_Int_error = 0.58%, N error 420.00% mu_up = 10.500000, mu_down = 10.500000, N_up = 10.500136, N_down = 10.500136, Delta_Int_error = 0.99%, N error 420.00% mu_up = 10.500000, mu_down = 10.500000, N_up = 11.017045, N_down = 11.017045, Delta_Int_error = 0.97%, N error 196.42% mu_up = 5.250000, mu_down = 5.250000, N_up = 6.733933, N_down = 6.733933, Delta_Int_error = 0.33%, N error 196.42% mu_up = 5.250000, mu_down = 5.250000, N_up = 6.027827, N_down = 6.027827, Delta_Int_error = 0.62%, N error 100.33% mu_up = 2.625000, mu_down = 7.875000, N_up = 4.000192, N_down = 8.000192, Delta_Int_error = 0.04%, N error 100.33% mu_up = 2.625000, mu_down = 7.875000, N_up = 2.999996, N_down = 8.999996, Delta_Int_error = 2.77%, N error 60.00% mu_up = 3.937500, mu_down = 6.562500, N_up = 3.999998, N_down = 7.999998, Delta_Int_error = 4.23%, N error 60.00% mu_up = 3.937500, mu_down = 6.562500, N_up = 4.000000, N_down = 7.000000, Delta_Int_error = 2.29%, N error 22.13% mu_up = 3.281250, mu_down = 7.218750, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.63%, N error 22.13% mu_up = 3.281250, mu_down = 7.218750, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 2.06%, N error 22.22% mu_up = 3.412500, mu_down = 7.218750, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 2.00%, N error 22.22% mu_up = 3.412500, mu_down = 7.218750, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.98%, N error 20.00% mu_up = 3.806250, mu_down = 7.546875, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.95%, N error 20.00% mu_up = 3.806250, mu_down = 7.546875, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.98%, N error 20.00% mu_up = 3.727500, mu_down = 7.481250, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.91%, N error 20.00% mu_up = 3.727500, mu_down = 7.481250, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.83%, N error 0.00% mu_up = 3.491250, mu_down = 7.284375, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.74%, N error 0.00% mu_up = 3.491250, mu_down = 7.284375, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.68%, N error 0.00% mu_up = 3.373125, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.66%, N error 0.00% mu_up = 3.373125, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.63%, N error 0.00% mu_up = 3.314062, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.61%, N error 0.00% mu_up = 3.314062, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.57%, N error 0.00% mu_up = 3.302250, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.52%, N error 0.00% mu_up = 3.302250, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.45%, N error 0.00% mu_up = 3.288075, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.39%, N error 0.00% mu_up = 3.288075, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.33%, N error 0.00% mu_up = 3.271065, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.27%, N error 0.00% mu_up = 3.271065, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.23%, N error 0.00% mu_up = 3.266983, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.20%, N error 0.00% mu_up = 3.266983, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.17%, N error 0.00% mu_up = 3.296376, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.14%, N error 0.00% mu_up = 3.296376, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.10%, N error 0.00% mu_up = 3.331648, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.07%, N error 0.00% mu_up = 3.331648, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.04%, N error 0.00% mu_up = 3.437464, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.02%, N error 0.00% mu_up = 3.437464, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 0.96%, N error 0.00% mu_up = 3.416300, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 0.98%, N error 0.00% mu_up = 3.416300, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 0.87%, N error 0.00% mu_up = 3.352811, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.35%, N error 0.00% mu_up = 3.352811, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 0.54%, N error 0.00% mu_up = 3.340113, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.06%, N error 0.00% mu_up = 3.340113, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 0.55%, N error 0.00% mu_up = 3.337066, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 0.67%, N error 0.00% mu_up = 3.337066, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 0.48%, N error 0.00% mu_up = 3.359008, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 1.06%, N error 0.00% mu_up = 3.359008, mu_down = 7.185938, N_up = 4.000000, N_down = 8.000000, Delta_Int_error = 0.10%, N error 0.00% --- Max iterations reached without convergence. (Any future claim to the contrary is a lie.) ---