матеріальних точок, записуючи їх в масиви xx [i], yy [i]. p align="justify">. У циклі перебирають всі матеріальні точки і визначають проекції прискорення, швидкості і координати для кожної з них у момент t + ? T :
axi (t + ? t) = Fxi (t + ? t)/mi, (t + ? t) = vix (t) + aix (t + ? t) ? t, (t + ? t) = xi (t) + vix (t + ? t) ? t.
З аналогічних формулами обчислюють проекції на вісь OY. Результати записують у масиви x [i], y [i], vx [i], vy [i]. p align="justify">. Перуть зображення матеріальних точок у попередній момент часу t, координати яких збережені в масивах xx [i], yy [i]. p align="justify">. На екрані будують точки в наступний момент t + ? T , або малюють графіки або виводять результат в числовому вигляді.
. Повернення до операції 2. Якщо цикл по t закінчився, - вихід з циклу. p align="justify"> 4. Комп'ютерна програма. Нижче приведений код програми, яка моделює рух 50 молекул газу в прямокутній посудині, що знаходиться в однорідному гравітаційному полі.
program PROGRAMMA4; dos, crt, graph; N = 20; dt = 0.01; m, Fx, Fy, x, y, vx, vy: array [1 .. N] of real;
Gd, Gm, i, j: integer;
ax, ay, F, l: real; Metka; Sila; {Обчислення діючих сил}
label Metka; i: = 1 to N do [i]: = 0; [i]: = 0;; i: = 1 to N doj: = 1 to N doj = iMetka;: = sqrt (sqr (x [i]-x [j]) + sqr (y [i]-y [j])); l <2l: = 2;: =-m [i] * m [j]/ sqr (l);
Fx [i]: = Fx [i] + F * (x [i]-x [j])/l-(random (20) +10); [i]: = Fy [i] + F * (y [i]-y [j])/l + m [i] * 10;
Metka:;; Nach_uslov;; {Завдання випадкових координат і швидкостей}
for i: = 1 to N do
m [i]: = 2; [i]: = random (80) +60; [i]: = random (80) +60; [i]: = random (30) + 15;
vx [i]: = random (30) +15;;; (Gd, Gm, 'c: bp bgi'); {Ініціалізація графіки} _uslov; Sila; i: = 1 to N do: = Fx [i]/m [i];
ay: = Fy [i]/m [i]; {Обчислення прискорень}
vx [i]: = vx [i] + ax * dt; [i]: = vy [i] + ay * dt; {швидкостей}
y [i]: = y [i] + vy [i] * dt; {координат}
if (y...
