Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- private void GetK4(double x, ref double y, double h)
- {
- //считам коэф. K по формулам
- double k1 = h * Func(x, y);
- double k2 = h * Func(x + h / 2, y + h / 2 * k1);
- double k3 = h * Func(x + h / 2, y + h / 2 * k2);
- double k4 = h * Func(x + h, y + h * k3);
- //получаем значение y
- y = y + (k1 + 2 * k2 + 2 * k3 + k4) / 6;
- }
- private void GetK5(double x, ref double y, double h)
- {
- //считам коэф. K по формулам
- double k1 = h * Func(x, y);
- double k2 = h * Func(x + h / 4, y + h / 4 * k1);
- double k3 = h * Func(x + h / 4, y + h / 8 * k1 + h / 8 * k2);
- double k4 = h * Func(x + h / 2, y - h / 2 * k2 + h * k3);
- double k5 = h * Func(x + 3 * h / 4, y + 3 * h / 16 * k1 + 9 * h / 16 * k4);
- double k6 = h * Func(x + h, y - 3 * h / 7 * k1 + 2 * h / 7 * k2
- + 12 * h / 3 * k3 - 12 * h / 7 * k4 + 8 * h / 7 * k5
- );
- //получаем значение y
- y = y + (7 * k1 + 32 * k3 + 12 * k4 + 32 * k5 + 7 * k6) / 90;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement