Advertisement
Guest User

Untitled

a guest
Dec 9th, 2018
68
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 6.47 KB | None | 0 0
  1. /home/denis/python-projects/computations/math-models-course-work/venv/bin/python /home/denis/python-projects/computations/math-models-course-work/main.py
  2. Solution:
  3. +-----+---------------+-----+---------------+----+---------------+-------------+------------+
  4. | x1 | x2 | x3 | x4 | p1 | p2 | p3 | p4 |
  5. +-----+---------------+-----+---------------+----+---------------+-------------+------------+
  6. | 0.1 | 210.656160256 | 3.9 | 10.4076526956 | 2 | 40.1017397756 | 0.000950625 | 0.00950625 |
  7. | 0.2 | 55.6621578947 | 3.8 | 5.44066204986 | 2 | 20.1974115789 | 0.00361 | 0.0361 |
  8. | 0.3 | 26.2235758258 | 3.7 | 3.80744564159 | 2 | 13.6210131644 | 0.007700625 | 0.07700625 |
  9. | 0.4 | 15.6737777778 | 3.6 | 3.00809876543 | 2 | 10.3731911111 | 0.01296 | 0.1296 |
  10. | 0.5 | 10.6794642857 | 3.5 | 2.54272959184 | 2 | 8.45457589286 | 0.019140625 | 0.19140625 |
  11. | 0.6 | 7.90733986928 | 3.4 | 2.24479469435 | 2 | 7.1991172549 | 0.02601 | 0.2601 |
  12. | 0.7 | 6.20098438466 | 3.3 | 2.04300646071 | 2 | 6.32170567641 | 0.033350625 | 0.33350625 |
  13. | 0.8 | 5.0715 | 3.2 | 1.9018125 | 2 | 5.68008 | 0.04096 | 0.4096 |
  14. | 0.9 | 4.28257755874 | 3.1 | 1.80155830154 | 2 | 5.19546577509 | 0.048650625 | 0.48650625 |
  15. | 1.0 | 3.70833333333 | 3.0 | 1.73055555556 | 2 | 4.82083333333 | 0.05625 | 0.5625 |
  16. | 1.1 | 3.27656105728 | 2.9 | 1.68148062372 | 2 | 4.52647273119 | 0.063600625 | 0.63600625 |
  17. | 1.2 | 2.94336507937 | 2.8 | 1.64957823129 | 2 | 4.2927847619 | 0.07056 | 0.7056 |
  18. | 1.3 | 2.68078309226 | 2.7 | 1.6316995621 | 2 | 4.1064033084 | 0.077000625 | 0.77000625 |
  19. | 1.4 | 2.47030926217 | 2.6 | 1.62575908707 | 2 | 3.9579843956 | 0.08281 | 0.8281 |
  20. | 1.5 | 2.29930555556 | 2.5 | 1.63041666667 | 2 | 3.84088541667 | 0.087890625 | 0.87890625 |
  21. | 1.6 | 2.15891666667 | 2.4 | 1.64488888889 | 2 | 3.75034666667 | 0.09216 | 0.9216 |
  22. | 1.7 | 2.0428094253 | 2.3 | 1.66883973646 | 2 | 3.68297036125 | 0.095550625 | 0.95550625 |
  23. | 1.8 | 1.94638496072 | 2.2 | 1.7023241506 | 2 | 3.6363840404 | 0.09801 | 0.9801 |
  24. | 1.9 | 1.86627288616 | 2.1 | 1.74577100489 | 2 | 3.60902381266 | 0.099500625 | 0.99500625 |
  25. | - | - | - | - | - | - | - | - |
  26. | 2.1 | 1.74577100489 | 1.9 | 1.86627288616 | 2 | 3.60902381266 | 0.099500625 | 0.99500625 |
  27. | 2.2 | 1.7023241506 | 1.8 | 1.94638496072 | 2 | 3.6363840404 | 0.09801 | 0.9801 |
  28. | 2.3 | 1.66883973646 | 1.7 | 2.0428094253 | 2 | 3.68297036125 | 0.095550625 | 0.95550625 |
  29. | 2.4 | 1.64488888889 | 1.6 | 2.15891666667 | 2 | 3.75034666667 | 0.09216 | 0.9216 |
  30. | 2.5 | 1.63041666667 | 1.5 | 2.29930555556 | 2 | 3.84088541667 | 0.087890625 | 0.87890625 |
  31. | 2.6 | 1.62575908707 | 1.4 | 2.47030926217 | 2 | 3.9579843956 | 0.08281 | 0.8281 |
  32. | 2.7 | 1.6316995621 | 1.3 | 2.68078309226 | 2 | 4.1064033084 | 0.077000625 | 0.77000625 |
  33. | 2.8 | 1.64957823129 | 1.2 | 2.94336507937 | 2 | 4.2927847619 | 0.07056 | 0.7056 |
  34. | 2.9 | 1.68148062372 | 1.1 | 3.27656105728 | 2 | 4.52647273119 | 0.063600625 | 0.63600625 |
  35. | 3.0 | 1.73055555556 | 1.0 | 3.70833333333 | 2 | 4.82083333333 | 0.05625 | 0.5625 |
  36. | 3.1 | 1.80155830154 | 0.9 | 4.28257755874 | 2 | 5.19546577509 | 0.048650625 | 0.48650625 |
  37. | 3.2 | 1.9018125 | 0.8 | 5.0715 | 2 | 5.68008 | 0.04096 | 0.4096 |
  38. | 3.3 | 2.04300646071 | 0.7 | 6.20098438466 | 2 | 6.32170567641 | 0.033350625 | 0.33350625 |
  39. | 3.4 | 2.24479469435 | 0.6 | 7.90733986928 | 2 | 7.1991172549 | 0.02601 | 0.2601 |
  40. | 3.5 | 2.54272959184 | 0.5 | 10.6794642857 | 2 | 8.45457589286 | 0.019140625 | 0.19140625 |
  41. | 3.6 | 3.00809876543 | 0.4 | 15.6737777778 | 2 | 10.3731911111 | 0.01296 | 0.1296 |
  42. | 3.7 | 3.80744564159 | 0.3 | 26.2235758258 | 2 | 13.6210131644 | 0.007700625 | 0.07700625 |
  43. | 3.8 | 5.44066204986 | 0.2 | 55.6621578947 | 2 | 20.1974115789 | 0.00361 | 0.0361 |
  44. | 3.9 | 10.4076526956 | 0.1 | 210.656160256 | 2 | 40.1017397756 | 0.000950625 | 0.00950625 |
  45. +-----+---------------+-----+---------------+----+---------------+-------------+------------+
  46.  
  47. Verification:
  48. +-----+--------+--------+--------+--------+----------+
  49. | x1 | dx1/dt | dx2/dt | dx3/dt | dx4/dt | Jacobian |
  50. +-----+--------+--------+--------+--------+----------+
  51. | 0.1 | 6e-16 | 0e+00 | 7e-12 | -7e-12 | -6e-16 |
  52. | 0.2 | -2e-16 | 0e+00 | 9e-14 | -8e-14 | 3e-15 |
  53. | 0.3 | -2e-17 | 4e-16 | 4e-13 | -4e-13 | -3e-15 |
  54. | 0.4 | -5e-16 | 0e+00 | 1e-13 | -1e-13 | 9e-15 |
  55. | 0.5 | -2e-16 | 0e+00 | -8e-15 | 1e-14 | 2e-15 |
  56. | 0.6 | -3e-16 | 0e+00 | -9e-15 | 4e-15 | -6e-15 |
  57. | 0.7 | 7e-17 | 2e-16 | 2e-14 | -2e-14 | 2e-15 |
  58. | 0.8 | 3e-16 | 0e+00 | 3e-14 | -3e-14 | -9e-15 |
  59. | 0.9 | -1e-16 | -2e-16 | 5e-15 | -7e-15 | -2e-14 |
  60. | 1.0 | 2e-16 | -2e-16 | 2e-15 | -2e-15 | -5e-15 |
  61. | 1.1 | 3e-16 | 0e+00 | 6e-15 | -7e-15 | -6e-15 |
  62. | 1.2 | 1e-17 | -1e-16 | -9e-16 | 1e-16 | 5e-15 |
  63. | 1.3 | 6e-16 | -2e-16 | 7e-15 | -7e-15 | -3e-15 |
  64. | 1.4 | 3e-17 | -1e-16 | 9e-16 | -3e-15 | -1e-14 |
  65. | 1.5 | -1e-15 | 1e-16 | -2e-15 | -1e-16 | 3e-15 |
  66. | 1.6 | 8e-16 | 6e-16 | -8e-16 | 2e-15 | -3e-15 |
  67. | 1.7 | 9e-17 | -2e-16 | 8e-16 | 2e-16 | -1e-14 |
  68. | 1.8 | -1e-15 | -4e-16 | -1e-15 | 4e-16 | -1e-14 |
  69. | 1.9 | 5e-16 | 3e-16 | -5e-16 | 6e-16 | 7e-15 |
  70. | 2.1 | -2e-16 | -3e-16 | 1e-15 | -6e-16 | 0e+00 |
  71. | 2.2 | 1e-15 | -2e-16 | -1e-15 | 2e-16 | 8e-15 |
  72. | 2.3 | -7e-16 | 4e-16 | 7e-16 | -1e-15 | -5e-15 |
  73. | 2.4 | 9e-16 | -8e-16 | 9e-16 | -1e-16 | -1e-14 |
  74. | 2.5 | 0e+00 | 0e+00 | -9e-16 | 0e+00 | 0e+00 |
  75. | 2.6 | -1e-15 | -6e-16 | -2e-15 | 1e-15 | 6e-15 |
  76. | 2.7 | -1e-15 | 7e-16 | -5e-16 | 2e-16 | 0e+00 |
  77. | 2.8 | -8e-16 | 2e-16 | -5e-15 | 4e-15 | 1e-14 |
  78. | 2.9 | -7e-16 | 7e-16 | -1e-15 | 1e-15 | 0e+00 |
  79. | 3.0 | 7e-16 | -4e-16 | 2e-16 | -4e-16 | -4e-15 |
  80. | 3.1 | -2e-15 | 0e+00 | -3e-15 | 3e-15 | 9e-15 |
  81. | 3.2 | 2e-15 | -2e-15 | 7e-16 | -7e-16 | 2e-15 |
  82. | 3.3 | -4e-16 | 0e+00 | -2e-15 | 1e-15 | 1e-15 |
  83. | 3.4 | -2e-15 | 2e-16 | -2e-15 | 2e-15 | 1e-15 |
  84. | 3.5 | -7e-16 | -3e-15 | 7e-16 | -9e-16 | -1e-15 |
  85. | 3.6 | -1e-15 | -2e-16 | -5e-15 | 5e-15 | 3e-15 |
  86. | 3.7 | 6e-15 | -1e-15 | -8e-15 | 8e-15 | -1e-15 |
  87. | 3.8 | 1e-15 | 2e-15 | 7e-15 | -6e-15 | 1e-15 |
  88. | 3.9 | -1e-14 | 6e-15 | 2e-14 | -2e-14 | -2e-16 |
  89. +-----+--------+--------+--------+--------+----------+
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement