бчислення правих частин}
begin [1]: = x [3]; [2]: = (fs-k23 * x [2])/k22; [3]: = (x [2]-x [3] * k2- x [1] * k3)/k1;; eu (var t, t0, h: real; var x0: mas; var x: mas);
{процедура, що виконує дії за методом Ейлера}
var j: integer; (t, x0, f); j: = 1 to 3 do [j]: = x0 [j] + h * f [j];: = t0 + h;
end;
begin
writeln ('введіть шлях до файлу з вихідними даними');
readln (l); (d, l); (d); (d, k1, k2, k3, k22, k23, fs); (d);
writeln ('введіть шлях до файлу виводу');
readln (l); (d, l); (d);: = 0.0; tk: = 5.0;: = 0.01;
{k1: = 1.0; k2: = 1.0; k3: = 1.0; k22: = 1.0; k23: = 1.0; fs: = 1.0;}: = 0; [1]: = 0.0; x0 [2] : = 0.0; x0 [3]: = 0.0;: = t0; ('t =', t: 2:2, '; x1 =', x0 [1]: 5:6, '; x2 =', x0 [3]: 5:6, '; x3 =', x0 [2]: 5:6, ';'); (d, 't =', t: 2:2, '; x1 =', x0 [ 1]: 5:6, '; x2 =', x0 [3]: 5:6, '; x3 =', x0 [2]: 5:6, ';'); t <= tk do (t, t0, h, x0, x);: = n +1; n = 10 then (d, 't =', t: 2:2, '; x1 =', x [1]: 5:6, '; x2 = ', x [3]: 5:6,'; x3 = ', x [2]: 5:6,'; '); (' t = ', t: 2:2,'; x1 = ' , x [1]: 5:6, '; x2 =', x [3]: 5:6, '; x3 =', x [2]: 5:6, ';');: = 0;; j: = 1 to 3 do0 [j]: = x [j];
t0: = t;
end;
close (d);
end.
Файл з вихідними даними
В
Результат роботи програми
tx 1 x 2 х 2 1 x 2 х 2 1 x 2 х 2 Графік залежності x (t)
В
Список використаної літератури
1. А.М. Епанешников, В.А. Епанешников. Програмування в середовищі TURBO PASCAL 7.0. - 3-е вид., Стер. - М.: В«ДІАЛОГ-МІФІВ», 2004 - 367с. p align="justify">. Вірт Н. Алгоритми і структури даних з прикладами на Паскалі. - 2-е вид. - СПб. : Невський діалект, 2007. - 352с. p align="justify">. Марченко А.І. Програмування в середовищі Turbo Pascal 7.0/А.І. Марченко, Л.А. Марченко. - 8-е вид. - Київ: Корона Принт, 2004. - 464с. p align="justify">. Фаронов, В.В. Turbo Pascal: навч. посібник для вузів/В.В. Фаронов. - М. [та ін]: Пітер, 2007. - 367с
. Шафрін Ю. А. Інформаційні технології: У 2 ч. - М.: Лабораторія Базових Знань, 2001. br/>