Advertisement
Guest User

Untitled

a guest
Dec 13th, 2017
73
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.37 KB | None | 0 0
  1. #include <iostream>
  2. #include <fstream>
  3. #include <math.h>
  4. using namespace std;
  5. double function_1(double t, double x, double y) { return (3*x+4*y) ; }
  6. double function_2(double t, double x, double y) { return (2*x+y); }
  7. double y_n = 1;
  8. double x_n = -1;
  9. double t_n = 0;
  10. double t_k = 1000;
  11. double n = 5000, h = 0.01;
  12. int main()
  13. {
  14.     string path = "E:\\university\\Math\\laba3\\laba3.txt";
  15.     ofstream file(path);
  16.     double k1 = 0, k2 = 0, k3 = 0, k4 = 0;
  17.     double l1 = 0, l2 = 0, l3 = 0, l4 = 0;
  18.     double y = y_n;
  19.     double x = x_n;
  20.     file << t_n <<'\t' << x << '\t' << y << '\n';
  21.     for (double i = t_n + h; i <= t_k; i += h)
  22.     {
  23.         k1 = function_1(i, x, y);
  24.         l1 = function_2(i, x, y);
  25.         k2 = function_1(i + (h / 2.0), x + ((h / 2.0)*k1), y + ((h / 2.0)*l1));
  26.         l2 = function_2(i + (h / 2.0), x + ((h / 2.0)*k1), y + ((h / 2.0)*l1));
  27.         k3 = function_1(i + (h / 2.0), x + ((h / 2.0)*k2), y + ((h / 2.0)*l2));
  28.         l3 = function_2(i + (h / 2.0), x + ((h / 2.0)*k2), y + ((h / 2.0)*l2));
  29.         k4 = function_1(i + h, x + h*k3, y + h*l3);
  30.         l4 = function_2(i + h, x + h*k3, y + h*l3);
  31.         x = x + (h / 6.0)*(k1 + 2.0*k2 + 2.0*k3 + k4);
  32.         y = y + (h / 6.0)*(l1 + 2.0*l2 + 2.0*l3 + l4);
  33.         file << i << '\t'<< x << '\t' << y << '\n';
  34.     }
  35.     file.close();
  36.     //system("notepad file.txt");
  37.     return 0;
  38. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement