Майер роберт Валерьевич основы компьютерного моделирования: Учебное пособие. Глазов: ггпи, 2005. 25 с

3. Алгоритм.
Дифференциальное уравнение второго порядка может быть 
решено методом конечных разностей Эйлера. Он состоит в том, что 
бесконечно малые приращения функции 
и 
ее первых двух 
производных 
заменяются конечными разностями, а 
бесконечно 
малые 
приращения 
-- 
малыми, 
но 
конечными 
приращениями 
Сначала, исходя из параметров системы 
, ее 

---

координаты и скорости 
в момент времени , расчитывается ее 
координата и скорость в следующий момент 
Затем это состояние рассматривается как исходное и процедура расчета 
повторяется для момента времени 
и так далее. Одновременно с 
вычислениями строятся графики 
Рассмотрим алгоритм численного решения уравнения. 
1. Задают параметры физической системы 
зависимость внешнего 
воздействия от времени 
а также начальные условия 
и шаг 
по времени 
2. Начало 
цикла 
по Дают 
приращение 
по 
времени: 
переменной присваивают значение 
3. Определяют ускорение, скорость и координату тела в момент 
4. Результаты 
вычислений 
выводят 
на 
экран в числовом виде либо строят соответствующие точки на координатной 
плоскости. 
5. Возвращение к операции 2. Если цикл по закончился, -- выход из цикла. 

---

4. Компьютерная программа.
При запуске программа рисует графики 
зависимостей 
координаты 
проекции 
скорости 
и 
проекции ускорения 
. 
Некоторые строчки программы заключены в скобки "(*" и "*)". Убрав скобки 
и активизировав соответствующие операторы, можно промоделировать 
различные явления. 
program PROGRAMMA2; 
uses dos, crt, graph; 
Const Fm=10;w=5;m=2;r=0;k=0; 
Mx=20; Mv=40; Ma=8; Mf=2; Mt=100; 
dt=0.00006; 
Var x,v,a,F,t : Real; 
j,xx,vv,aa,FF,tt,Gd,Gm : Integer; 
BEGIN 
Gd:= Detect; 
InitGraph(Gd, Gm, 'c:\bp\bgi'); 
if GraphResult<>grOk then Halt(1); 
t:=0; v:=0; x:=-3; 
line(30,300,650,300); 
line(31,500,31,10); 
OutTextXY(50,20,'X, V, A'); 
Repeat 
begin {Задание функции F=F(t)} 
t:=t+dt; (* F:=Fm*sin(w*t); *) 
(*If sin(w*t)<0 then F:=0; 
If sin(w*t)>0 then F:=Fm;*) 
F:=0; If t<1 then F:= Fm; 
If t>3 then F:=-Fm; 
a:=(F-r*v-k*x)/m; x:=x+v*dt; v:=v+a*dt; tt:=round(t*Mt); 
xx:=round(x*Mx); vv:=round(v*Mv); aa:=round(a*Ma); FF:=round(F*Mf); 
circle(30+tt,300-xx,1); circle(30+tt,300-vv,1); circle(30+tt,300-aa,2); 
end; 
until KeyPressed; 
CloseGraph; 
END.