private static void Test_A(double x1) { CP cpO = new CP(0.0, 1.0, 0.0); ODE_2_A RG_A = new ODE_2_A(cpO, f_OA); CP cp1 = RG_A.Solve_with_Precision(x1, eps); S1 = S_OA(x1); err = Math.Abs(S1 - cp1.y); SW.WriteLine($"\r\n Test - MAC_ODE_Order_2_RungeKutta_4_A :"); SW.WriteLine($"{cp1.x,8:F4} {cp1.y,14:F9} {S1,14:F9} {err,11:E1} {RG_A.iter}"); }
private static void Test_Bh(double x1) { CP cpO = new CP(0.0, 36.0 / 37.0, 3.0 / 37.0); ODE_2_B RG_B = new ODE_2_B(cpO, f_OBh); CP cp1 = RG_B.Solve_with_Precision(x1, eps); S1 = S_OBh(x1); err = Math.Abs(S1 - cp1.y); SW.WriteLine($"\r\n Test - MAC_ODE_Order_2_RungeKutta_4_B :"); SW.WriteLine($"{cp1.x,8:F4} {cp1.y,14:F9} {S1,14:F9} {err,11:E1} {RG_B.iter}"); }