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}");
}
