Реферат: Компьютерное моделирование движения тел

IF h < 0 THEN 200

GOTO 100

Где к - коэффициент сопротивления воздуха, который подбирается экспериментальным путем и может быть изменен в соответствии с условиями задачи, причем для вертикальной и горизонтальной составляющей может быть введен свой коэффициент, а выражение (4.3) обеспечивает пересчет скорости используя корень квадратного выражения.

Кроме этого, для большей информационности, выведем на экран значение скорости, высоты, дальности и времени:

LOCATE 1, 6: PRINT "V m/s"

LOCATE 2, 6: PRINT "H m"

LOCATE 2, 24: PRINT "S m"

LOCATE 1, 24: PRINT "T s"

LOCATE 1, 12: PRINT v

LOCATE 1, 28: PRINT t

LOCATE 2, 12: PRINT h

LOCATE 2, 28: PRINT s

Причем, первую половину необходимо разместить в программе до 80-й строки (кадра), а вторую в конце тела цикла.

Полностью программа представлена в приложении 3. Следует обратить внимание на то, что значения коэффициентов k и q в приложении 3 изменены для большей реалистичности полета ядра, так как часть ресурсов компьютера задействована для отображения текущих параметров, а выводимые значения высоты h и дальности s делятся на 5,5 и 2 соответственно для соблюдения масштаба.

К случайным событиям можно отнести и роль оператора на ход решения задачи. В данном случае изменим угол полета ядра уже в ходе решения задачи.

Что бы иметь возможность в процессе работы программы изменять угол стрельбы задействуем еще три функциональные клавиши F2, F3, F5, где F2 будет отвечать за подъем на один градус в верх, F3 – за один градус вниз, F5 – позволит прервать выполнение программы в любое время. Для этого в программе после выполнения подпрограммы выполнения фейерверка вставим блок:

ERASE x,y

210 KEY(2) ON

ON KEY(2) GOSUB 220

KEY(3) ON

ON KEY(3) GOSUB 230

KEY(1) ON

ON KEY(1) GOSUB 72

GOTO 210

GOTO 210

220 a = a + 1

LOCATE 1, 46: PRINT a

GOTO 210

230 a = a - 1