Слідкуюча система літакового витратоміра
Д. Каханер, К. Моулер, С. Нэш, Численные методы и программное обеспечение, М., ``Мир'', 1998. Лук’яненко С.О. Чиселові методи в інформатиці: Навч. посіб. — К.: НТУУ «КПІ», 2007. — 140с. Підпрограма розв’язування диференціальних рівнянь системи методом Рунге-Кутта}. Задаємо число рівнянь, похибку, початковий крок інтегрування, початкові умови. Дж. Форсайт, М. Малькольм, К. Моулер, Машинные методы… Читати ще >
Слідкуюча система літакового витратоміра (реферат, курсова, диплом, контрольна)
Міністерство освіти і науки України Національний університет «Львівська політехніка»
ІКТА Кафедра БІТ Курсова робота з курсу: «Комп'ютерні методи дослідження інформаційних процесів та систем»
на тему: «Слідкуюча система літакового витратоміра»
Виконав:
студент групи Перевірив:
Львів 2011
Зміст
1. Завдання
2. Теоретичні відомості
3. Зведення рівнянь
4. Лістинги програм
5. Результати
6. Графік перехідного процесу Список літератури
1. Завдання
слідкуючий витратомір літак рівняння Тема № 6
Слідкуюча система літакового витратоміра
Схема
Рівняння ланок
Варіант 11
Параметри | ||
Тм (сек) | 0,1 | |
Те (сек) | 0,04 | |
С (рад/сек) | ||
Ky | ||
S (рад/сек) | ||
i | ||
1) Звести систему алгебро-диференціальних рівнянь до системи трьох диференціальних рівнянь першого порядку, представивши її у нормальній формі та розв’язати цю систему вказаними методами. Початкові умови — иin(0)=1 радіан, решта початкових умов — нульові. Числові значення сталих параметрів, заданих в таблиці, виразити з допомогою одиниць системи СІ.
2) Побудувати графік зміни величини иout(t)
Метод Мерсона-Рунге-Кутта та Рунге-Кутта 4 порядку з автоматичним кроком
2. Теоретичні відомості
Методи з автоматичною зміною кроку
Застосовуються в тому випадку, якщо розв’язок потрібно одержати із заданою точністю. При високій точності (похибка) автоматична зміна кроку забезпечує зменшення загального числа кроків в декілька разів (особливо при розв’язках у вигляді кривих, що сильно відрізняються крутизною).
Метод Рунге-Кутта з автоматичною зміною кроку
Після обчислення з кроком всі обчислення виконуються повторно з кроком. Після цього порівнюються результати, отримані в точці хn+1 з кроком і. Якщо модуль різниці менший, то обчислення продовжуються з кроком, в іншому випадку крок зменшують. Якщо нерівність дуже сильна, то крок збільшують.
Маємо
— значення незалежної змінної в точці
— значення функції в точці
— значення функції в точці, обчислене з кроком
— значення функції в точці, обчислене з кроком
— значення функції, обчислене з кроком
1) Якщо обчислення повторюються з кроком і т.д., доки не виконається умова .
2) Якщо виконується ця умова, то можливі два варіанти, в залежності від значення K, де K — ознака поділу кроку.
Початкове значенняі залишається таким після першого поділу кроку на два. Надалі, якщо крок ділиться, то K приймає значення одиниці.
а) Якщо, то навіть коли виконалась умова, крок не змінюється, тобто лишається тим самим (обчислення далі проводяться з попереднім кроком).
б) Якщо і виконалась умова, тоді .
В обох випадках а) і б) результат виводиться на друк.
Метод Рунге-Кутта-Мерсона з автоматичною зміною кроку
Метод дозволяє оцінити похибку на кожному кроці інтегрування. При цьому не потрібно зберігати в пам’яті обчислення значень функцій на кроці і для оцінки похибки.
Алгоритм методу
1. Задаємо число рівнянь, похибку, початковий крок інтегрування, початкові умови .
2. За допомогою п’яти циклів з керуючою змінною обчислюємо коефіцієнти
3. Знаходимо значення та похибку
4. Перевіряємо виконання умов Можливі випадки:
а) Якщо перша умова не виконується, тобто, то ділимо крок на 2 та повторюємо обчислення з п. 2, встановивши початкові значення .
б) Якщо виконується перша та друга умови, значення та виводяться на друк.
Якщо друга умова не виконується, крок збільшується вдвічі і тоді обчислення знову повторюється з п. 2.
Треба відмітити, що похибка на кожному кроці методу Рунге-Кутта-Мерсона оцінюється приблизно. При розв’язуванні нелінійних ДР істинна похибка може відрізнятися в декілька разів від заданої .
де .
— крок поділити на 2 і повернутися на початок.
для всіх рівнянь: виводимо на друк, а крок збільшуємо удвічі.
3. Зведення рівнянь
З рівняння редуктора маємо:
З рівняння двигуна:
Розв’яжемо відносно :
Тобто ми отримаємо три рівняння для нашої системи.
Виконаємо заміну :
иout=y1;
y=y2;
Звідси отримаємо систему трьох диференційних рівнянь:
4. Лістинги програм
Метод Рунге-Кутта 4 порядку
program rk611 (Input, Output);
Uses crt;
const n=3;
type mas=array[1.n] of real;
mmas=array[1.4] of real;
var i, k, j:integer;
ep, TM, TE, KU, C, S, h, hh, x, x0,II, xk: real;
f, y, yy, z, y0,y1,y2,yy1:mas;
w, p: mmas;
ff:text;
label 4,5;
{Праві частини диференціальних рівнянь}
procedure drr (y:mas; var f: mas);
BEGIN
f[1]: =y[2];
f[2]:=(C*S*KU*(1-y[3])-y[1]-TM*y[2])/(TM*TE);
f[3]:=y[1]/II;
END;
{Підпрограма розв’язування диференціальних рівнянь системи методом Рунге-Кутта}
procedure dyf (yy, f: mas; w, p: mmas; var yy1: mas);
var xx: real;
BEGIN
xx:=x;
for j:=1 to n do
BEGIN
y[j]: =yy[j];
z[j]:=y[j];
END;
for i:=1 to 4 do
BEGIN
x:=xx+w[i]*h;
drr (yy, f);
for j:=1 to n do
BEGIN
yy[j]: =yy[j]+p[i]*h*f[j];
y[j]:=z[j]+w[i+1]*h*f[j];
END;
END;
for i:=1 to n do yy1[i]: =yy[i];
END;
procedure wk (var w, p: mmas);
BEGIN
for i:=1 to 4 do
BEGIN
w[i]: =trunc (i/2)/2;
p[i]:=(1+2*w[i]-2*trunc (i/4))/6;
END;
END;
{ Підпрограма автоматичного вибору кроку }
procedure run_kut (yy, f: mas; w, p: mmas; var yy1: mas);
var kk: integer;
label 1,2,6;
BEGIN
for i:=1 to n do y0[i]: =yy[i];
kk:=0; x0:=x; hh:=h;
dyf (yy, f, w, p, yy);
for i:=1 to n do
BEGIN
y1[i]: =yy[i]; yy[i]: =y0[i];
END;
6:x:=x0; h:=hh/2;
dyf (yy, f, w, p, yy);
for i:=1 to n do y2[i]: =yy[i];
dyf (yy, f, w, p, yy);
for i:=1 to n do if abs (yy[i]-y1[i])>ep then goto 1;
for i:=1 to n do yy1[i]: =yy[i];
goto 2;
1:hh:=h;
kk:=1; {Поділ кроку навпіл}
for i:=1 to n do
BEGIN
yy[i]: =y0[i];
y1[i]:=y2[i];
END;
goto 6;
2:h:=hh;
if kk=0 then h:=2*hh; {Збільшення кроку вдвічі}
END;
BEGIN {Основна програма}
assign (ff,'D:RK611.txt');
rewrite (ff);
writeln (ff,' ');
writeln (ff,' Runge_Kutta4 Automation');
{ writeln ('Порядок рівняння N=3'); {read (n);}
x:=0; xk:=9;
{write (' Введіть крок H='); read (h);}h:=1;
for i:=1 to n do
BEGIN
{ write ('Введіть початкове значення Y (', i,')='); read (yy[i])};
yy[1]: =0;yy[2]:=0;yy[3]:=0;
END;
ep:=1.0e-4;
TM:=0.1;TE:=0.04;C:=2;KU:=1;S:=100;II:=25;
k:=0;
wk (w, p);
5:run_kut (yy, f, w, p, yy);
k:=k+1;
if k<>100 then goto 4;
writeln (ff,' T=', x,' Yout=', yy[3]);
k:=0;
4:if x<=xk then goto 5;
repeat until keypressed;
close (ff);
END.
Метод Рунге-Кутта-Мерсона
program rkm611;
const n=3;a=0;b=9;
type vector=array[1.n] of real;
var i, kk, dd, j: integer;
ks:longint;
x, x0, St, Ct, dS, dC, h, hh, eps: real;
Tm, Te, C, S, Ky, II: real;
f, y, y0,y1,y2,R:vector;
ff:text;
procedure dy;
begin
f[1]: =y[2];
f[2]:=(C*S*Ky*(1-y[3]-y[1]-Tm*y[2]))/(Tm*Te);
f[3]:=y[1]/II;
end;
procedure rkm;
var
k0,k1,k2,k3,k4,z:vector;
begin
dy;
for i:=1 to n do
begin
k0[i]: =h*f[i];
z[i]:=y[i];
y[i]:=z[i]+k0[i]/3;
end;
x:=x+h/3;
for i:=1 to n do
begin
k1[i]: =h*f[i];
y[i]:=z[i]+k0[i]/6+k1[i]/6;
end;
dy;
for i:=1 to n do
begin
k2[i]: =h*f[i];
y[i]:=z[i]+k0[i]/8+3*k2[i]/8;
end;
x:=x+h/3;
dy;
for i:=1 to n do
begin
k3[i]: =h*f[i];
y[i]:=z[i]+k0[i]/2−3*k2[i]/2+2*k3[i];
end;
x:=x+0.5*h;
for i:=1 to n do
begin
k4[i]: =h*f[i];
y[i]:=z[i]+(k0[i]/2−3*k2[i]/2+k4[i])/6;
R[i]:=(-2*k0[i]+9*k2[i]+8*k3[i]+k4[i])/30;
end;
end;
procedure aut_step (var y, y0, y1,y2,f:vector);
label 1;
begin
for i:=1 to n do
y0[i]: =y[i];x0:=x;
1: x:=x0;
rkm;
for i:=1 to n do
if abs (R[i])>eps then
begin
h:=h/2;
for i:=1 to n do
y[i]: =y0[i];
goto 1;
end;
dd:=0;
for i:=1 to n do
if abs (R[i])<(eps/30) then dd:=dd+1;
if dd=n then
begin
h:=2*h;
dd:=0;
end;
end;
begin
assign (ff,'D:rkm611.txt');
rewrite (ff);
writeln (ff, ' ');
writeln (ff,' Runge-Kutta-Mersona');
y[1]: =0;y[2]:=0;y[3]:=0;x:=a;h:=1;eps:=1e-2;
Tm:=0.1;Te:=0.04;C:=2;Ky:=1;S:=100;II:=25;
ks:=0;
repeat
ks:=ks+1;
aut_step (y, y0, y1,y2,f);
if ks=100 then
begin
writeln (ff, ' x=', x,' Uout=', y[3],' h=', h);
ks:=0;
end;
until x>b;
close (ff);
end.
5. Результати
Метод Рунге-Кутта 4 порядку
T= 3.51 757 8125E-03 Yout= 1.405 491 8036E-05
T= 6.103 515 6250E-03 Yout= 1.98 822 1132E-04
T= 9.155 273 4375E-03 Yout= 3.619 936 9461E-04
T= 1.220 703 1250E-02 Yout= 8.374 488 0165E-04
T= 1.525 878 9063E-02 Yout= 1.596 460 7579E-03
T= 1.831 054 6875E-02 Yout= 2.692 948 9732E-03
T= 2.136 230 4688E-02 Yout= 4.175 042 5822E-03
T= 2.441 406 2500E-02 Yout= 6.85 539 6149E-03
T= 2.746 582 0313E-02 Yout= 8.462 325 8697E-03
T= 3.51 757 8125E-02 Yout= 1.133 875 7190E-02
T= 3.356 933 5938E-02 Yout= 1.474 400 8511E-02
T= 3.662 109 3750E-02 Yout= 1.870 339 2738E-02
T= 3.967 285 1563E-02 Yout= 2.323 865 2229E-02
T= 4.272 460 9375E-02 Yout= 2.836 822 5410E-02
T= 4.827 880 8594E-02 Yout= 3.927 680 4451E-02
T= 5.438 232 4219E-02 Yout= 5.367 281 4390E-02
T= 6.48 583 9844E-02 Yout= 7.64 282 1180E-02
T= 6.658 935 5469E-02 Yout= 9.20 681 3944E-02
T= 7.269 287 1094E-02 Yout= 1.123 550 3126E-01
T= 7.879 638 6719E-02 Yout= 1.370 518 6637E-01
T= 8.489 990 2344E-02 Yout= 1.642 392 1295E-01
T= 9.100 341 7969E-02 Yout= 1.938 393 5109E-01
T= 9.710 693 3594E-02 Yout= 2.257 574 5474E-01
T= 1.32 104 4922E-01 Yout= 2.598 837 8816E-01
T= 1.153 564 4531E-01 Yout= 3.338 657 6254E-01
T= 1.275 634 7656E-01 Yout= 4.154 297 8490E-01
T= 1.397 705 0781E-01 Yout= 5.29 840 0884E-01
T= 1.588 134 7656E-01 Yout= 6.483 469 9905E-01
T= 1.832 275 3906E-01 Yout= 8.427 552 3013E-01
T= 1.978 759 7656E-01 Yout= 9.594 168 0716E-01
T= 2.100 830 0781E-01 Yout= 1.54 514 1358E+00
T= 2.222 900 3906E-01 Yout= 1.146 285 6313E+00
T= 2.344 970 7031E-01 Yout= 1.233 503 3241E+00
T= 2.467 041 0156E-01 Yout= 1.315 035 3052E+00
T= 2.589 111 3281E-01 Yout= 1.389 860 1611E+00
T= 2.711 181 6406E-01 Yout= 1.457 079 2229E+00
T= 2.833 251 9531E-01 Yout= 1.515 926 5784E+00
T= 2.955 322 2656E-01 Yout= 1.565 776 8078E+00
T= 3.77 392 5781E-01 Yout= 1.606 150 4207E+00
T= 3.199 462 8906E-01 Yout= 1.636 717 0017E+00
T= 3.321 533 2031E-01 Yout= 1.657 296 0914E+00
T= 3.443 603 5156E-01 Yout= 1.667 855 8615E+00
T= 3.565 673 8281E-01 Yout= 1.668 509 6685E+00
T= 3.687 744 1406E-01 Yout= 1.659 510 5951E+00
T= 3.809 814 4531E-01 Yout= 1.641 244 1167E+00
T= 3.931 884 7656E-01 Yout= 1.614 219 0483E+00
T= 4.53 955 0781E-01 Yout= 1.579 056 9504E+00
T= 4.176 025 3906E-01 Yout= 1.536 480 1889E+00
T= 4.298 095 7031E-01 Yout= 1.487 298 8563E+00
T= 4.427 490 2344E-01 Yout= 1.428 941 4557E+00
T= 4.671 630 8594E-01 Yout= 1.305 342 9340E+00
T= 4.915 771 4844E-01 Yout= 1.170 893 2533E+00
T= 5.157 470 7031E-01 Yout= 1.35 138 2699E+00
T= 5.279 541 0156E-01 Yout= 9.680 442 3305E-01
T= 5.401 611 3281E-01 Yout= 9.31 880 7154E-01
T= 5.523 681 6406E-01 Yout= 8.414 452 8809E-01
T= 5.645 751 9531E-01 Yout= 7.836 233 1994E-01
T= 5.767 822 2656E-01 Yout= 7.304 517 0419E-01
T= 5.889 892 5781E-01 Yout= 6.825 736 2647E-01
T= 6.11 962 8906E-01 Yout= 6.405 389 3802E-01
T= 6.134 033 2031E-01 Yout= 6.47 987 0297E-01
T= 6.256 103 5156E-01 Yout= 5.757 013 1393E-01
T= 6.378 173 8281E-01 Yout= 5.534 901 9027E-01
T= 6.500 244 1406E-01 Yout= 5.383 030 5041E-01
T= 6.622 314 4531E-01 Yout= 5.301 727 2678E-01
T= 6.744 384 7656E-01 Yout= 5.290 294 7102E-01
T= 6.866 455 0781E-01 Yout= 5.347 046 7662E-01
T= 6.988 525 3906E-01 Yout= 5.469 359 2770E-01
T= 7.110 595 7031E-01 Yout= 5.653 732 6638E-01
T= 7.232 666 0156E-01 Yout= 5.895 865 5555E-01
T= 7.354 736 3281E-01 Yout= 6.190 738 0255E-01
T= 7.501 220 7031E-01 Yout= 6.606 200 7028E-01
T= 7.745 361 3281E-01 Yout= 7.419 687 7648E-01
T= 8.118 896 4844E-01 Yout= 8.845 535 6240E-01
T= 8.363 037 1094E-01 Yout= 9.810 461 3978E-01
T= 8.571 777 3438E-01 Yout= 1.60 601 9432E+00
T= 8.693 847 6563E-01 Yout= 1.104 490 4102E+00
T= 8.815 917 9688E-01 Yout= 1.145 710 2525E+00
T= 8.937 988 2813E-01 Yout= 1.183 738 0524E+00
T= 9.60 058 5938E-01 Yout= 1.218 110 2646E+00
T= 9.182 128 9063E-01 Yout= 1.248 428 3119E+00
T= 9.304 199 2188E-01 Yout= 1.274 362 6081E+00
T= 9.426 269 5313E-01 Yout= 1.295 655 4816E+00
T= 9.548 339 8438E-01 Yout= 1.312 122 9839E+00
T= 9.670 410 1563E-01 Yout= 1.323 655 5884E+00
T= 9.792 480 4688E-01 Yout= 1.330 217 7969E+00
T= 9.914 550 7813E-01 Yout= 1.331 846 6910E+00
T= 1.3 662 1094E+00 Yout= 1.328 649 4745E+00
T= 1.15 869 1406E+00 Yout= 1.320 800 0710E+00
T= 1.28 076 1719E+00 Yout= 1.308 534 8499E+00
T= 1.40 283 2031E+00 Yout= 1.292 147 5649E+00
T= 1.63 598 6328E+00 Yout= 1.250 684 6286E+00
T= 1.88 012 6953E+00 Yout= 1.195 628 9038E+00
T= 1.128 051 7578E+00 Yout= 1.89 410 3221E+00
T= 1.155 151 3672E+00 Yout= 1.13 975 2222E+00
T= 1.179 565 4297E+00 Yout= 9.487 757 2167E-01
T= 1.203 979 4922E+00 Yout= 8.898 858 0764E-01
T= 1.218 505 8594E+00 Yout= 8.591 037 8587E-01
T= 1.230 712 8906E+00 Yout= 8.361 835 9866E-01
T= 1.242 919 9219E+00 Yout= 8.162 316 7070E-01
T= 1.255 126 9531E+00 Yout= 7.994 569 4453E-01
T= 1.267 333 9844E+00 Yout= 7.860 170 3435E-01
T= 1.279 541 0156E+00 Yout= 7.760 173 1398E-01
T= 1.291 748 0469E+00 Yout= 7.695 107 8186E-01
T= 1.303 955 0781E+00 Yout= 7.664 986 8468E-01
T= 1.316 162 1094E+00 Yout= 7.669 318 7018E-01
T= 1.336 547 8516E+00 Yout= 7.750 466 0483E-01
T= 1.360 961 9141E+00 Yout= 7.959 787 8690E-01
T= 1.385 375 9766E+00 Yout= 8.273 821 0937E-01
T= 1.411 254 8828E+00 Yout= 8.695 038 7284E-01
T= 1.460 083 0078E+00 Yout= 9.624 986 3688E-01
T= 1.484 985 3516E+00 Yout= 1.10 790 0446E+00
T= 1.509 399 4141E+00 Yout= 1.55 037 0478E+00
T= 1.533 813 4766E+00 Yout= 1.93 852 8997E+00
T= 1.558 227 5391E+00 Yout= 1.125 321 1051E+00
T= 1.582 641 6016E+00 Yout= 1.148 041 7763E+00
T= 1.607 055 6641E+00 Yout= 1.161 183 9965E+00
T= 1.631 469 7266E+00 Yout= 1.164 503 1962E+00
T= 1.655 883 7891E+00 Yout= 1.158 324 3851E+00
T= 1.680 297 8516E+00 Yout= 1.143 493 9749E+00
T= 1.704 711 9141E+00 Yout= 1.121 304 5638E+00
T= 1.741 333 0078E+00 Yout= 1.77 885 0115E+00
T= 1.790 161 1328E+00 Yout= 1.11 223 9278E+00
T= 1.815 795 8984E+00 Yout= 9.770 348 1210E-01
T= 1.840 209 9609E+00 Yout= 9.475 196 4549E-01
T= 1.864 624 0234E+00 Yout= 9.225 211 1814E-01
T= 1.889 038 0859E+00 Yout= 9.32 202 8518E-01
T= 1.913 452 1484E+00 Yout= 8.904 135 3831E-01
T= 1.937 866 2109E+00 Yout= 8.844 875 5829E-01
T= 1.962 280 2734E+00 Yout= 8.854 187 9983E-01
T= 1.986 694 3359E+00 Yout= 8.927 962 5633E-01
T= 2.11 108 3984E+00 Yout= 9.58 649 9816E-01
T= 2.50 415 0391E+00 Yout= 9.361 533 7266E-01
T= 2.112 426 7578E+00 Yout= 9.951 766 2467E-01
T= 2.146 118 1641E+00 Yout= 1.26 094 9065E+00
T= 2.170 532 2266E+00 Yout= 1.45 539 5650E+00
T= 2.194 946 2891E+00 Yout= 1.61 387 0948E+00
T= 2.219 360 3516E+00 Yout= 1.72 927 1161E+00
T= 2.243 774 4141E+00 Yout= 1.79 730 5241E+00
T= 2.268 188 4766E+00 Yout= 1.81 659 1179E+00
T= 2.292 602 5391E+00 Yout= 1.78 858 2379E+00
T= 2.317 504 8828E+00 Yout= 1.71 551 2194E+00
T= 2.366 333 0078E+00 Yout= 1.46 914 4745E+00
T= 2.433 715 8203E+00 Yout= 1.1 894 6780E+00
T= 2.482 543 9453E+00 Yout= 9.716 397 2891E-01
T= 2.509 399 4141E+00 Yout= 9.586 389 5199E-01
T= 2.533 813 4766E+00 Yout= 9.499 058 8179E-01
T= 2.558 227 5391E+00 Yout= 9.444 685 3601E-01
T= 2.582 641 6016E+00 Yout= 9.424 629 5409E-01
T= 2.620 239 2578E+00 Yout= 9.458 736 2610E-01
T= 2.669 067 3828E+00 Yout= 9.602 693 8957E-01
T= 2.744 750 9766E+00 Yout= 9.948 223 4606E-01
T= 2.793 579 1016E+00 Yout= 1.16 827 9985E+00
T= 2.842 407 2266E+00 Yout= 1.32 935 5234E+00
T= 2.891 235 3516E+00 Yout= 1.40 238 8433E+00
T= 2.940 063 4766E+00 Yout= 1.38 043 1648E+00
T= 2.988 891 6016E+00 Yout= 1.27 791 4971E+00
T= 3.75 805 6641E+00 Yout= 9.996 151 1478E-01
T= 3.124 633 7891E+00 Yout= 9.848 610 1583E-01
T= 3.173 461 9141E+00 Yout= 9.749 006 5645E-01
T= 3.222 290 0391E+00 Yout= 9.714 145 6129E-01
T= 3.271 118 1641E+00 Yout= 9.745 035 4723E-01
T= 3.352 172 8516E+00 Yout= 9.900 723 7706E-01
T= 3.431 274 4141E+00 Yout= 1.8 279 5460E+00
T= 3.480 102 5391E+00 Yout= 1.16 315 6318E+00
T= 3.528 930 6641E+00 Yout= 1.19 984 5181E+00
T= 3.577 758 7891E+00 Yout= 1.18 931 5131E+00
T= 3.674 438 4766E+00 Yout= 1.6 461 8904E+00
T= 3.771 606 4453E+00 Yout= 9.913 604 9308E-01
T= 3.820 434 5703E+00 Yout= 9.869 540 9281E-01
T= 3.870 727 5391E+00 Yout= 9.858 525 5325E-01
T= 3.976 196 2891E+00 Yout= 9.934 373 1612E-01
T= 4.92 407 2266E+00 Yout= 1.6 243 0137E+00
T= 4.190 063 4766E+00 Yout= 1.9 979 9414E+00
T= 4.306 274 4141E+00 Yout= 1.3 733 5254E+00
T= 4.423 461 9141E+00 Yout= 9.949 261 5879E-01
T= 4.521 118 1641E+00 Yout= 9.930 874 8894E-01
T= 4.679 321 2891E+00 Yout= 1.546 5007E+00
T= 4.779 907 2266E+00 Yout= 1.4 634 2970E+00
T= 4.888 305 6641E+00 Yout= 1.3 894 4016E+00
T= 5.72 875 9766E+00 Yout= 9.971 931 3315E-01
T= 5.198 852 5391E+00 Yout= 9.970 802 6368E-01
T= 5.421 508 7891E+00 Yout= 1.2 348 6395E+00
T= 5.655 883 7891E+00 Yout= 9.994 278 1170E-01
T= 5.855 102 5391E+00 Yout= 9.987 298 9263E-01
T= 6.110 961 9141E+00 Yout= 1.1 235 1811E+00
T= 6.400 024 4141E+00 Yout= 9.991 192 5001E-01
T= 6.800 415 0391E+00 Yout= 1.498 2377E+00
T= 7.288 696 2891E+00 Yout= 1.216 5761E+00
T= 7.874 633 7891E+00 Yout= 1.29 5823E+00
T= 8.812 133 7891E+00 Yout= 1.1 6801E+00
Метод Рунге-Кутта-Мерсона
x= 1.335 144 0430E-03 Yout= 8.577 261 1741E-08
x= 2.670 288 0858E-03 Yout= 6.802 679 1949E-07
x= 4.5 432 1282E-03 Yout= 2.274 929 4689E-06
x= 5.340 576 1726E-03 Yout= 5.342 591 8950E-06
x= 6.675 720 2170E-03 Yout= 1.33 807 5534E-05
x= 8.10 864 2607E-03 Yout= 1.769 876 2403E-05
x= 9.346 008 3009E-03 Yout= 2.784 515 4122E-05
x= 1.68 115 2341E-02 Yout= 4.118 141 2361E-05
x= 1.201 629 6381E-02 Yout= 5.809 588 2336E-05
x= 1.335 144 0421E-02 Yout= 7.896 159 9903E-05
x= 1.468 658 4461E-02 Yout= 1.41 367 8276E-04
x= 1.602 172 8504E-02 Yout= 1.339 653 0622E-04
x= 1.735 687 2553E-02 Yout= 1.687 771 6414E-04
x= 1.869 201 6601E-02 Yout= 2.88 889 1532E-04
x= 2.2 716 0650E-02 Yout= 2.546 041 1594E-04
x= 2.136 230 4699E-02 Yout= 3.62 137 3845E-04
x= 2.269 744 8747E-02 Yout= 3.639 965 7730E-04
x= 2.403 259 2796E-02 Yout= 4.282 196 4198E-04
x= 2.536 773 6845E-02 Yout= 4.991 385 3767E-04
x= 2.670 288 0893E-02 Yout= 5.769 978 3405E-04
x= 2.803 802 4942E-02 Yout= 6.620 314 2242E-04
x= 2.937 316 8991E-02 Yout= 7.544 628 6176E-04
x= 3.70 831 3039E-02 Yout= 8.545 057 1381E-04
x= 3.204 345 7078E-02 Yout= 9.623 638 6764E-04
x= 3.337 860 1109E-02 Yout= 1.78 231 8542E-03
x= 3.471 374 5141E-02 Yout= 1.202 295 1508E-03
x= 3.604 888 9173E-02 Yout= 1.334 730 4761E-03
x= 3.738 403 3204E-02 Yout= 1.475 706 0758E-03
x= 3.871 917 7236E-02 Yout= 1.625 381 9998E-03
x= 4.5 432 1267E-02 Yout= 1.783 910 3697E-03
x= 4.138 946 5299E-02 Yout= 1.951 435 6389E-03
x= 4.272 460 9331E-02 Yout= 2.128 094 8440E-03
x= 4.405 975 3362E-02 Yout= 2.314 017 8484E-03
x= 4.539 489 7394E-02 Yout= 2.509 327 5783E-03
x= 4.673 004 1425E-02 Yout= 2.714 140 2512E-03
x= 4.806 518 5457E-02 Yout= 2.928 565 5974E-03
x= 4.940 032 9489E-02 Yout= 3.152 707 0744E-03
x= 5.73 547 3520E-02 Yout= 3.386 662 0746E-03
x= 5.207 061 7552E-02 Yout= 3.630 522 1262E-03
x= 5.340 576 1583E-02 Yout= 3.884 373 0885E-03
x= 5.474 090 5615E-02 Yout= 4.148 295 3403E-03
x= 5.607 604 9647E-02 Yout= 4.422 363 9628E-03
x= 5.741 119 3678E-02 Yout= 4.706 648 9167E-03
x= 5.874 633 7710E-02 Yout= 5.1 215 2134E-03
x= 6.8 148 1741E-02 Yout= 5.306 123 0817E-03
x= 6.141 662 5773E-02 Yout= 5.621 428 1277E-03
x= 6.275 176 9811E-02 Yout= 5.947 181 4914E-03
x= 6.408 691 3877E-02 Yout= 6.283 429 9976E-03
x= 6.542 205 7942E-02 Yout= 6.630 216 3012E-03
x= 6.675 720 2008E-02 Yout= 6.987 579 0298E-03
x= 6.809 234 6074E-02 Yout= 7.355 552 9201E-03
x= 6.942 749 0139E-02 Yout= 7.734 168 9511E-03
x= 7.76 263 4205E-02 Yout= 8.123 454 4722E-03
x= 7.276 980 0750E-02 Yout= 8.728 805 3140E-03
x= 7.544 008 8813E-02 Yout= 9.571 721 7104E-03
x= 7.811 037 6876E-02 Yout= 1.45 763 8844E-02
x= 8.78 066 4939E-02 Yout= 1.138 662 1413E-02
x= 8.345 095 3002E-02 Yout= 1.235 868 9523E-02
x= 8.612 124 1066E-02 Yout= 1.337 382 1774E-02
x= 8.879 152 9129E-02 Yout= 1.443 195 8157E-02
x= 9.146 181 7192E-02 Yout= 1.553 300 2763E-02
x= 9.413 210 5255E-02 Yout= 1.667 682 6334E-02
x= 9.680 239 3318E-02 Yout= 1.786 326 8646E-02
x= 9.947 268 1382E-02 Yout= 1.909 214 0748E-02
x= 1.21 429 6944E-01 Yout= 2.36 322 7061E-02
x= 1.48 132 5751E-01 Yout= 2.167 628 7348E-02
x= 1.74 835 4557E-01 Yout= 2.303 105 8553E-02
x= 1.101 538 3363E-01 Yout= 2.442 725 6535E-02
x= 1.128 241 2170E-01 Yout= 2.586 457 7687E-02
x= 1.154 944 0976E-01 Yout= 2.734 270 0458E-02
x= 1.181 646 9782E-01 Yout= 2.886 128 6775E-02
x= 1.208 349 8589E-01 Yout= 3.41 998 3379E-02
x= 1.235 052 7395E-01 Yout= 3.201 842 3074E-02
x= 1.261 755 6204E-01 Yout= 3.365 622 5904E-02
x= 1.288 458 5018E-01 Yout= 3.533 300 0247E-02
x= 1.315 161 3831E-01 Yout= 3.704 834 3843E-02
x= 1.341 864 2644E-01 Yout= 3.880 184 4766E-02
x= 1.368 567 1457E-01 Yout= 4.59 308 2316E-02
x= 1.395 270 0270E-01 Yout= 4.242 162 7874E-02
x= 1.421 972 9083E-01 Yout= 4.428 704 5687E-02
x= 1.448 675 7896E-01 Yout= 4.618 889 3614E-02
x= 1.475 378 6709E-01 Yout= 4.812 672 3819E-02
x= 1.502 081 5523E-01 Yout= 5.10 008 3418E-02
x= 1.528 784 4336E-01 Yout= 5.210 851 5093E-02
x= 1.555 487 3149E-01 Yout= 5.415 155 7653E-02
x= 1.582 190 1962E-01 Yout= 5.622 874 6574E-02
x= 1.608 893 0775E-01 Yout= 5.833 961 4491E-02
x= 1.635 595 9588E-01 Yout= 6.48 369 1667E-02
x= 1.662 298 8401E-01 Yout= 6.266 050 6425E-02
x= 1.689 001 7215E-01 Yout= 6.486 958 5560E-02
x= 1.715 704 6028E-01 Yout= 6.711 045 4712E-02
x= 1.742 407 4841E-01 Yout= 6.938 263 8727E-02
x= 1.769 110 3654E-01 Yout= 7.168 566 1982E-02
x= 1.795 813 2467E-01 Yout= 7.401 904 8704E-02
x= 1.822 516 1280E-01 Yout= 7.638 232 3250E-02
x= 1.849 219 0093E-01 Yout= 7.877 501 0380E-02
x= 1.875 921 8907E-01 Yout= 8.119 663 5507E-02
x= 1.902 624 7720E-01 Yout= 8.364 672 4930E-02
x= 1.929 327 6533E-01 Yout= 8.612 480 6058E-02
x= 1.956 030 5346E-01 Yout= 8.863 040 7603E-02
x= 1.982 733 4159E-01 Yout= 9.116 305 9780E-02
x= 2.9 436 2972E-01 Yout= 9.372 229 4471E-02
x= 2.36 139 1785E-01 Yout= 9.630 764 5398E-02
x= 2.62 842 0599E-01 Yout= 9.891 864 8268E-02
x= 2.89 544 9412E-01 Yout= 1.15 548 4092E-01
x= 2.116 247 8225E-01 Yout= 1.42 157 6344E-01
x= 2.142 950 7038E-01 Yout= 1.69 009 5831E-01
x= 2.169 653 5851E-01 Yout= 1.96 099 7049E-01
x= 2.196 356 4664E-01 Yout= 1.123 423 4753E-01
x= 2.223 059 3477E-01 Yout= 1.150 976 3969E-01
x= 2.249 762 2291E-01 Yout= 1.178 753 9998E-01
x= 2.276 465 1104E-01 Yout= 1.206 751 8430E-01
x= 2.303 167 9917E-01 Yout= 1.234 965 5144E-01
x= 2.329 870 8730E-01 Yout= 1.263 390 6322E-01
x= 2.356 573 7543E-01 Yout= 1.292 022 8452E-01
x= 2.383 276 6356E-01 Yout= 1.320 857 8333E-01
x= 2.409 979 5169E-01 Yout= 1.349 891 3081E-01
x= 2.436 682 3983E-01 Yout= 1.379 119 0136E-01
x= 2.463 385 2796E-01 Yout= 1.408 536 7260E-01
x= 2.490 088 1609E-01 Yout= 1.438 140 2550E-01
x= 2.516 791 0413E-01 Yout= 1.467 925 4432E-01
x= 2.543 493 9213E-01 Yout= 1.497 888 1670E-01
x= 2.570 196 8012E-01 Yout= 1.528 024 3369E-01
x= 2.596 899 6812E-01 Yout= 1.558 329 8974E-01
x= 2.623 602 5611E-01 Yout= 1.588 800 8272E-01
x= 2.650 305 4411E-01 Yout= 1.619 433 1400E-01
x= 2.677 008 3210E-01 Yout= 1.650 222 8838E-01
x= 2.703 711 2010E-01 Yout= 1.681 166 1417E-01
x= 2.730 414 0809E-01 Yout= 1.712 259 0317E-01
x= 2.757 116 9609E-01 Yout= 1.743 497 7068E-01
x= 2.783 819 8408E-01 Yout= 1.774 878 3552E-01
x= 2.810 522 7208E-01 Yout= 1.806 397 2002E-01
x= 2.837 225 6007E-01 Yout= 1.838 050 5004E-01
x= 2.863 928 4807E-01 Yout= 1.869 834 5494E-01
x= 2.890 631 3606E-01 Yout= 1.901 745 6762E-01
x= 2.917 334 2406E-01 Yout= 1.933 780 2449E-01
x= 2.944 037 1205E-01 Yout= 1.965 934 6548E-01
x= 2.970 740 0005E-01 Yout= 1.998 205 3402E-01
x= 2.997 442 8804E-01 Yout= 2.30 588 7706E-01
x= 3.24 145 7604E-01 Yout= 2.63 081 4506E-01
x= 3.50 848 6403E-01 Yout= 2.95 679 9193E-01
x= 3.77 551 5203E-01 Yout= 2.128 380 7511E-01
x= 3.104 254 4002E-01 Yout= 2.161 180 5548E-01
x= 3.130 957 2802E-01 Yout= 2.194 075 9740E-01
x= 3.157 660 1601E-01 Yout= 2.227 063 6868E-01
x= 3.184 363 0401E-01 Yout= 2.260 140 4056E-01
x= 3.211 065 9200E-01 Yout= 2.293 302 8772E-01
x= 3.237 768 8000E-01 Yout= 2.326 547 8825E-01
x= 3.264 471 6799E-01 Yout= 2.359 872 2363E-01
x= 3.291 174 5599E-01 Yout= 2.393 272 7874E-01
x= 3.317 877 4398E-01 Yout= 2.426 746 4182E-01
x= 3.344 580 3198E-01 Yout= 2.460 290 0446E-01
x= 3.371 283 1997E-01 Yout= 2.493 900 6162E-01
x= 3.397 986 0797E-01 Yout= 2.527 575 1154E-01
x= 3.424 688 9596E-01 Yout= 2.561 310 5579E-01
x= 3.451 391 8396E-01 Yout= 2.595 103 9924E-01
x= 3.478 094 7195E-01 Yout= 2.628 952 5000E-01
x= 3.504 797 5995E-01 Yout= 2.662 853 1946E-01
x= 3.531 500 4794E-01 Yout= 2.696 803 2224E-01
x= 3.558 203 3594E-01 Yout= 2.730 799 7618E-01
x= 3.584 906 2393E-01 Yout= 2.764 840 0231E-01
x= 3.611 609 1193E-01 Yout= 2.798 921 2485E-01
x= 3.638 311 9992E-01 Yout= 2.833 040 7117E-01
x= 3.665 014 8792E-01 Yout= 2.867 195 7181E-01
x= 3.691 717 7591E-01 Yout= 2.901 383 6040E-01
x= 3.718 420 6391E-01 Yout= 2.935 601 7371E-01
x= 3.745 123 5190E-01 Yout= 2.969 847 5155E-01
x= 3.771 826 3990E-01 Yout= 3.4 118 3685E-01
x= 3.798 529 2789E-01 Yout= 3.38 411 7556E-01
x= 3.825 232 1589E-01 Yout= 3.72 725 1663E-01
x= 3.851 935 0388E-01 Yout= 3.107 056 1206E-01
x= 3.878 637 9188E-01 Yout= 3.141 402 1682E-01
x= 3.905 340 7987E-01 Yout= 3.175 760 8884E-01
x= 3.932 043 6787E-01 Yout= 3.210 129 8900E-01
x= 3.958 746 5586E-01 Yout= 3.244 506 8111E-01
x= 3.985 449 4386E-01 Yout= 3.278 889 3189E-01
x= 4.12 152 3185E-01 Yout= 3.313 275 1094E-01
x= 4.38 855 1985E-01 Yout= 3.347 661 9073E-01
x= 4.65 558 0784E-01 Yout= 3.382 047 4657E-01
x= 4.92 260 9584E-01 Yout= 3.416 429 5660E-01
x= 4.118 963 8383E-01 Yout= 3.450 806 0178E-01
x= 4.145 666 7183E-01 Yout= 3.485 174 6582E-01
x= 4.172 369 5982E-01 Yout= 3.519 533 3524E-01
x= 4.199 072 4782E-01 Yout= 3.553 879 9927E-01
x= 4.225 775 3581E-01 Yout= 3.588 212 4987E-01
x= 4.252 478 2381E-01 Yout= 3.622 528 8174E-01
x= 4.279 181 1180E-01 Yout= 3.656 826 9222E-01
x= 4.305 883 9980E-01 Yout= 3.691 104 8135E-01
x= 4.332 586 8779E-01 Yout= 3.725 360 5178E-01
x= 4.359 289 7579E-01 Yout= 3.759 592 0882E-01
x= 4.385 992 6378E-01 Yout= 3.793 797 6036E-01
x= 4.412 695 5178E-01 Yout= 3.827 975 1690E-01
x= 4.439 398 3977E-01 Yout= 3.862 122 9147E-01
x= 4.466 101 2777E-01 Yout= 3.896 238 9970E-01
x= 4.492 804 1576E-01 Yout= 3.930 321 5969E-01
x= 4.519 507 0376E-01 Yout= 3.964 368 9207E-01
x= 4.546 209 9175E-01 Yout= 3.998 379 1998E-01
x= 4.572 912 7975E-01 Yout= 4.32 350 6898E-01
x= 4.599 615 6774E-01 Yout= 4.66 281 6712E-01
x= 4.626 318 5574E-01 Yout= 4.100 170 4486E-01
x= 4.653 021 4373E-01 Yout= 4.134 015 3507E-01
x= 4.679 724 3173E-01 Yout= 4.167 814 7301E-01
x= 4.706 427 1972E-01 Yout= 4.201 566 9632E-01
x= 4.733 130 0772E-01 Yout= 4.235 270 4498E-01
x= 4.759 832 9571E-01 Yout= 4.268 923 6132E-01
x= 4.786 535 8371E-01 Yout= 4.302 524 8995E-01
x= 4.813 238 7170E-01 Yout= 4.336 072 7782E-01
x= 4.839 941 5970E-01 Yout= 4.369 565 7412E-01
x= 4.866 644 4769E-01 Yout= 4.403 002 3032E-01
x= 4.893 347 3569E-01 Yout= 4.436 381 0012E-01
x= 4.920 050 2368E-01 Yout= 4.469 700 3942E-01
x= 4.946 753 1168E-01 Yout= 4.502 959 0636E-01
x= 4.973 455 9967E-01 Yout= 4.536 155 6125E-01
x= 5.158 8767E-01 Yout= 4.569 288 6655E-01
x= 5.26 861 7594E-01 Yout= 4.602 356 8688E-01
x= 5.53 564 6421E-01 Yout= 4.635 358 8899E-01
x= 5.80 267 5247E-01 Yout= 4.668 293 4174E-01
x= 5.106 970 4074E-01 Yout= 4.701 159 1606E-01
x= 5.133 673 2901E-01 Yout= 4.733 954 8498E-01
x= 5.160 376 1728E-01 Yout= 4.766 679 2359E-01
x= 5.187 079 0555E-01 Yout= 4.799 331 0901E-01
x= 5.213 781 9381E-01 Yout= 4.831 909 2037E-01
x= 5.240 484 8208E-01 Yout= 4.864 412 3882E-01
x= 5.267 187 7035E-01 Yout= 4.896 839 4750E-01
x= 5.293 890 5862E-01 Yout= 4.929 189 3148E-01
x= 5.320 593 4688E-01 Yout= 4.961 460 7784E-01
x= 5.347 296 3515E-01 Yout= 4.993 652 7555E-01
x= 5.373 999 2342E-01 Yout= 5.25 764 1552E-01
x= 5.400 702 1169E-01 Yout= 5.57 793 9055E-01
x= 5.427 404 9996E-01 Yout= 5.89 740 9532E-01
x= 5.454 107 8822E-01 Yout= 5.121 604 2640E-01
x= 5.480 810 7649E-01 Yout= 5.153 382 8219E-01
x= 5.507 513 6476E-01 Yout= 5.185 075 6291E-01
x= 5.534 216 5303E-01 Yout= 5.216 681 7062E-01
x= 5.560 919 4130E-01 Yout= 5.248 200 0916E-01
x= 5.587 622 2956E-01 Yout= 5.279 629 8419E-01
x= 5.614 325 1783E-01 Yout= 5.310 970 0308E-01
x= 5.641 028 0610E-01 Yout= 5.342 219 7501E-01
x= 5.667 730 9437E-01 Yout= 5.373 378 1085E-01
x= 5.694 433 8263E-01 Yout= 5.404 444 2320E-01
x= 5.721 136 7090E-01 Yout= 5.435 417 2639E-01
x= 5.753 447 1959E-01 Yout= 5.472 768 9747E-01
x= 5.806 852 9558E-01 Yout= 5.534 201 9858E-01
x= 5.860 258 7157E-01 Yout= 5.595 248 9499E-01
x= 5.913 664 4756E-01 Yout= 5.655 903 7646E-01
x= 5.967 070 2355E-01 Yout= 5.716 160 5819E-01
x= 6.20 475 9954E-01 Yout= 5.776 013 8035E-01
x= 6.73 881 7553E-01 Yout= 5.835 458 0758E-01
x= 6.127 287 5152E-01 Yout= 5.894 488 2852E-01
x= 6.180 693 2751E-01 Yout= 5.953 099 5531E-01
x= 6.234 099 0350E-01 Yout= 6.11 287 2320E-01
x= 6.287 504 7949E-01 Yout= 6.69 046 9001E-01
x= 6.340 910 5548E-01 Yout= 6.126 374 3575E-01
x= 6.394 316 3147E-01 Yout= 6.183 265 6208E-01
x= 6.447 722 0746E-01 Yout= 6.239 716 9201E-01
x= 6.505 578 3150E-01 Yout= 6.300 371 9014E-01
x= 6.612 389 8402E-01 Yout= 6.410 968 8634E-01
x= 6.719 201 3655E-01 Yout= 6.519 753 0493E-01
x= 6.826 012 8907E-01 Yout= 6.626 703 1539E-01
x= 6.932 824 4160E-01 Yout= 6.731 800 8003E-01
x= 7.39 635 9412E-01 Yout= 6.835 030 4138E-01
x= 7.146 447 4665E-01 Yout= 6.936 379 0985E-01
x= 7.253 258 9918E-01 Yout= 7.35 836 5183E-01
x= 7.360 070 5170E-01 Yout= 7.133 394 7804E-01
x= 7.466 882 0423E-01 Yout= 7.229 048 3225E-01
x= 7.573 693 5675E-01 Yout= 7.322 793 8031E-01
x= 7.680 505 0928E-01 Yout= 7.414 629 9944E-01
x= 7.787 316 6180E-01 Yout= 7.504 557 6798E-01
x= 7.894 128 1433E-01 Yout= 7.592 579 5525E-01
x= 8.939 6685E-01 Yout= 7.678 700 1187E-01
x= 8.107 751 1938E-01 Yout= 7.762 925 6032E-01
x= 8.214 562 7191E-01 Yout= 7.845 263 8577E-01
x= 8.321 374 2443E-01 Yout= 7.925 724 2725E-01
x= 8.428 185 7696E-01 Yout= 8.4 317 6903E-01
x= 8.534 997 2948E-01 Yout= 8.81 056 3235E-01
x= 8.641 808 8201E-01 Yout= 8.155 953 6737E-01
x= 8.748 620 3453E-01 Yout= 8.229 024 4544E-01
x= 8.855 431 8706E-01 Yout= 8.300 284 5154E-01
x= 8.962 243 3959E-01 Yout= 8.369 750 7708E-01
x= 9.69 054 9211E-01 Yout= 8.437 441 1285E-01
x= 9.175 866 4464E-01 Yout= 8.503 374 4237E-01
x= 9.282 677 9716E-01 Yout= 8.567 570 3523E-01
x= 9.389 489 4969E-01 Yout= 8.630 049 4094E-01
x= 9.496 301 0221E-01 Yout= 8.690 832 8281E-01
x= 9.603 112 5474E-01 Yout= 8.749 942 5215E-01
x= 9.709 924 0727E-01 Yout= 8.807 401 0265E-01
x= 9.816 735 5979E-01 Yout= 8.863 231 4501E-01
x= 9.923 547 1232E-01 Yout= 8.917 457 4174E-01
x= 1.3 035 8647E+00 Yout= 8.970 103 0218E-01
x= 1.13 717 0167E+00 Yout= 9.21 192 7773E-01
x= 1.24 398 1686E+00 Yout= 9.70 751 5724E-01
x= 1.35 079 3206E+00 Yout= 9.118 804 6268E-01
x= 1.45 760 4726E+00 Yout= 9.165 377 4482E-01
x= 1.56 441 6246E+00 Yout= 9.210 495 7930E-01
x= 1.67 122 7766E+00 Yout= 9.254 185 6268E-01
x= 1.77 803 9285E+00 Yout= 9.296 473 0882E-01
x= 1.88 485 0805E+00 Yout= 9.337 384 4534E-01
x= 1.99 166 2325E+00 Yout= 9.376 946 1021E-01
x= 1.109 847 3845E+00 Yout= 9.415 184 4858E-01
x= 1.120 528 5365E+00 Yout= 9.452 126 0976E-01
x= 1.131 209 6884E+00 Yout= 9.487 797 4421E-01
x= 1.141 890 8404E+00 Yout= 9.522 225 0086E-01
x= 1.152 571 9924E+00 Yout= 9.555 435 2448E-01
x= 1.163 253 1444E+00 Yout= 9.587 454 5310E-01
x= 1.173 934 2964E+00 Yout= 9.618 309 1578E-01
x= 1.184 615 4483E+00 Yout= 9.648 025 3025E-01
x= 1.195 296 6003E+00 Yout= 9.676 629 0085E-01
x= 1.205 977 7523E+00 Yout= 9.704 146 1654E-01
x= 1.216 658 9043E+00 Yout= 9.730 602 4901E-01
x= 1.227 340 0563E+00 Yout= 9.756 023 5089E-01
x= 1.238 021 2082E+00 Yout= 9.780 434 5414E-01
x= 1.248 702 3602E+00 Yout= 9.803 860 6841E-01
x= 1.259 383 5122E+00 Yout= 9.826 326 7963E-01
x= 1.270 064 6642E+00 Yout= 9.847 857 4869E-01
x= 1.280 745 8162E+00 Yout= 9.868 477 1006E-01
x= 1.291 426 9681E+00 Yout= 9.888 209 7070E-01
x= 1.302 108 1201E+00 Yout= 9.907 079 0891E-01
x= 1.312 789 2721E+00 Yout= 9.925 108 7337E-01
x= 1.323 470 4241E+00 Yout= 9.942 321 8217E-01
x= 1.334 151 5761E+00 Yout= 9.958 741 2198E-01
x= 1.344 832 7280E+00 Yout= 9.974 389 4725E-01
x= 1.355 513 8800E+00 Yout= 9.989 288 7954E-01
x= 1.366 195 0320E+00 Yout= 1.346 1068E+00
x= 1.376 876 1840E+00 Yout= 1.1 692 7830E+00
x= 1.387 557 3360E+00 Yout= 1.2 971 0272E+00
x= 1.398 238 4879E+00 Yout= 1.4 182 9237E+00
x= 1.408 919 6399E+00 Yout= 1.5 330 5211E+00
x= 1.419 600 7919E+00 Yout= 1.6 415 8323E+00
x= 1.430 281 9439E+00 Yout= 1.7 440 8340E+00
x= 1.440 963 0959E+00 Yout= 1.8 407 4668E+00
x= 1.451 644 2478E+00 Yout= 1.9 317 6345E+00
x= 1.462 325 3998E+00 Yout= 1.10 173 2045E+00
x= 1.473 006 5518E+00 Yout= 1.10 976 0072E+00
x= 1.483 687 7038E+00 Yout= 1.11 727 8367E+00
x= 1.494 368 8558E+00 Yout= 1.12 430 4497E+00
x= 1.505 050 0077E+00 Yout= 1.13 085 5666E+00
x= 1.515 731 1597E+00 Yout= 1.13 694 8708E+00
x= 1.526 412 3117E+00 Yout= 1.14 260 0093E+00
x= 1.537 093 4637E+00 Yout= 1.14 782 5927E+00
x= 1.547 774 6157E+00 Yout= 1.15 264 1949E+00
x= 1.558 455 7676E+00 Yout= 1.15 706 3540E+00
x= 1.569 136 9196E+00 Yout= 1.16 110 5723E+00
x= 1.579 818 0716E+00 Yout= 1.16 478 3161E+00
x= 1.590 499 2236E+00 Yout= 1.16 811 0164E+00
x= 1.601 180 3756E+00 Yout= 1.17 110 0695E+00
x= 1.611 861 5275E+00 Yout= 1.17 376 8363E+00
x= 1.633 223 8326E+00 Yout= 1.17 818 8001E+00
x= 1.654 586 1376E+00 Yout= 1.18 147 0965E+00
x= 1.675 948 4427E+00 Yout= 1.18 371 3994E+00
x= 1.697 310 7477E+00 Yout= 1.18 500 8623E+00
x= 1.718 673 0528E+00 Yout= 1.18 544 1320E+00
x= 1.740 035 3579E+00 Yout= 1.18 509 3634E+00
x= 1.761 397 6629E+00 Yout= 1.18 404 2343E+00
x= 1.782 759 9680E+00 Yout= 1.18 235 9608E+00
x= 1.804 122 2730E+00 Yout= 1.18 011 3136E+00
x= 1.825 484 5781E+00 Yout= 1.17 736 6337E+00
x= 1.846 846 8831E+00 Yout= 1.17 417 8489E+00
x= 1.884 017 2931E+00 Yout= 1.16 774 2641E+00
x= 1.926 741 9021E+00 Yout= 1.15 924 1156E+00
x= 1.969 466 5111E+00 Yout= 1.14 989 1077E+00
x= 2.12 191 1204E+00 Yout= 1.13 999 0070E+00
x= 2.76 705 2806E+00 Yout= 1.12 454 0729E+00
x= 2.162 154 4986E+00 Yout= 1.10 416 1703E+00
x= 2.420 495 9702E+00 Yout= 1.5 171 0198E+00
x= 3.104 089 7210E+00 Yout= 9.998 576 9550E-01
6. Графік перехідного процесу
Рунге-Кутта 4 порядку
Рунге-Кутта-Мерсона
Список літератури
1. С. М. Єжов, Методи обчислень, К. ВПЦ «Київський університет», 2000.
2. А. А. Самарский, Введение в численные методы, М., Наука, 1987.
3. Е. А. Волков, Численные методы, М., Наука, 1987.
4. Р. В. Хемминг, Численные методы, М., Наука, 1972.
5. Дж. Ортега, У. Пул, Введение в численные методы решения дифференциальных уравнений, М., Наука, 1986.
6. Дж. Форсайт, М. Малькольм, К. Моулер, Машинные методы вычислений, М., ``Мир'', 1980.
7. Д. Каханер, К. Моулер, С. Нэш, Численные методы и программное обеспечение, М., ``Мир'', 1998.
8.Лук'яненко С.О. Чиселові методи в інформатиці: Навч. посіб. — К.: НТУУ «КПІ», 2007. — 140с.
9.Фельдман Л. П., Петренко А.І., Дмітрієва О.А. Чисельні методи в інформатиці. — К.: Видавнича група BNV, 2006.-480с.
10. Лук’яненко С.О. «Основи обчислювальних методів розв’язування диференціальних рівнянь», Навч. Посібник. — К.: ІСДО, 1998 — 212с.