.1 t [i]: = t [1] + (i-1) * h;
.2 v [i]: = v [i-1] + h * f (t [i-1]);
Алгоритм функції fun1
1. fun1: = 12/(2100 +0.56 * w-60)
Алгоритм функції fun2
1. fun2: = (2100-20 * w-60)/12
4. Схема алгоритму
Функція fun1
В
: real
Функція fun2
В
: real
Процедура EILER
В
: integer; tn, tk, vn: real; f: fun;
Var t, v: vect
5. Таблиця ідентифікаторів
Найменування фізичний смислідентіфікаторВремя t t Кутова швидкість двигуна ? w Постійні коефіцієнти a 1 , a 2 , b 1 , b 2 a1, a2, b1, b2 Момент опору на валу двигуна M C Mc Момент інерції двигуна JJ Крок часу D tdt Відносна помилка вимірювання ? epsilon Певне і встановлений час t опр , ? вуст topr, wopr Певна і встановлена ​​кутова швидкість ? опр wopr Число інтервалів розбиття nn
6. Текст програми
kyrsa4; crt; fun = function (y: real): real; = array [1 .. 100] of real; f1, f2: fun;, w: vect; f3, f4: text ;, n, m: integer;, wn, wk, Mc, a1, a2, b1, b2, dt, epsilon, j, topr, wopr: real;
{$ F +} fun1 (w: real): real;: = 12/(2100 +0.56 * w-60); fun2 (w: real): real;: = (2100-20 * w-60)/12;
{$ F-} EILER (n: integer; tn, tk, vn: real; f: fun; var t, v: vect); i: integer; h: real; [1] : = tn; t [n +1]: = tk; v [1]: = vn;: = (t [n +1] - t [1])/n; i: = 2 to n +1 do begin [i]: = t [1] + (...
