Advertisement
Guest User

Untitled

a guest
Apr 19th, 2015
212
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C# 1.22 KB | None | 0 0
  1.         private void GetK4(double x, ref double y, double h)
  2.         {
  3.             //считам коэф. K по формулам
  4.             double k1 = h * Func(x, y);
  5.             double k2 = h * Func(x + h / 2, y + h / 2 * k1);
  6.             double k3 = h * Func(x + h / 2, y + h / 2 * k2);
  7.             double k4 = h * Func(x + h, y + h * k3);
  8.  
  9.             //получаем значение y
  10.             y = y + (k1 + 2 * k2 + 2 * k3 + k4) / 6;
  11.         }
  12.  
  13.         private void GetK5(double x, ref double y, double h)
  14.         {
  15.             //считам коэф. K по формулам
  16.             double k1 = h * Func(x, y);
  17.             double k2 = h * Func(x + h / 4, y + h / 4 * k1);
  18.             double k3 = h * Func(x + h / 4, y + h / 8 * k1 + h / 8 * k2);
  19.             double k4 = h * Func(x + h / 2, y - h / 2 * k2 + h * k3);
  20.             double k5 = h * Func(x + 3 * h / 4, y + 3 * h / 16 * k1 + 9 * h / 16 * k4);
  21.             double k6 = h * Func(x + h, y - 3 * h / 7 * k1 + 2 * h / 7 * k2
  22.                 + 12 * h / 3 * k3 - 12 * h / 7 * k4 + 8 * h / 7 * k5
  23.                 );
  24.             //получаем значение y
  25.             y = y + (7 * k1 + 32 * k3 + 12 * k4 + 32 * k5 + 7 * k6) / 90;
  26.         }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement