Курсовая работа: Моделирование физических процессов

c:=c+g/6;

end;

readln;

gd:=detect;

init('c:\tp70\bgi');

setbkcolor(15);

setcolor(3);

line(0,240,640,240);

line(320,0,320,480);

for i:=1 to n do

begin

if i+1<=n then

begin

setlinestyle(0,0,3);

setcolor(7);

line(320+Xf[i],Yf[i]-160,320+Xf[i+1],Yf[i+1]-160);

putpixel(460,240,15);

putpixel(320,40,15);

putpixel(Xf[i]+320,Yf[i]-160,0);

end;

end;

readln;

closeGraph;

end.

Таблица измерений в Pascal, Mathcad:

t	X1	X2	Xm
2	2,077
2,1	2,16
2,2	2,247
2,3	2,341
2,4	2,44
2,5	2,545
2,6	2,656
2,7	2,773
2,8	2,896

X1 – метод Эйлера модифицированный, X2 – метод Рунге – Кутта, Xm – решение в Mathcad

Решение в Mathcad.

Вывод