RANDOMIZE TIMER

50 RANDOMIZE TIMER

mt = mt + 1

IF mt > 12 THEN s = 0 AND ws = 0

ds = 0

IF mt > 12 THEN y = y + 1

IF y > f THEN END

'ПРЕРЫВАНИЕ ПРОГРАММЫ

IF mt > 12 THEN PRINT "++++++++++++++++++++++++++++++++++++++++++++++++++++++"

IF mt > 12 THEN INPUT "Если введёте 999 то программа остановится-"; i

IF i = 999 THEN GOTO 3000

IF mt > 12 THEN mt = 1

IF mt = 1 THEN INPUT "Введите цену килограмма пшеницы в $-"; pw

'МЕСЯЦ СБОРА УРОЖАЯ ИЛИ НЕТ

dw = 0

IF mt = 8 THEN dw = dw + h ELSE dw = dw + 0

'ЕСТЕСТВЕННАЯ СМЕРТНОСТЬ МЫШЕЙ

dm = 0

dm = dm - 1 * INT(m / 12)

'ОПРЕДЕЛЕНИЕ КОЛ-ВА ПШЕНИЦЫ ЗА ВЫЧЕТОМ СЪЕД. МЫШАМИ

dw = dw - 2 * INT(m)

IF w <= 0 THEN END

'СКОЛЬКО ПШЕНИЦЫ ЗАБИРАЕТСЯ В МЕСЯЦ

IF w > t THEN dt = t ELSE dt = w

dw = dw - dt

'КОЛИЧЕСТВО МЫШЕЙ,СЪЕДЕННЫХ ЗА МЕСЯЦ

SELECT CASE INT(m)

CASE IS > 60 * c

dm = dm - 40 * c

CASE IS <= 60 * c

dm = dm - 15 * c

CASE 0

dm = dm + 0

END SELECT

'ПО КОЛ-ВУ ЗЕРНА ОПРЕДЕЛЯЕТСЯ ПРИРОСТ МЫШЕЙ

IF m <= w / 2 THEN dm = dm + INT(1.5 * m) ELSE dm = dm + INT(m / 10)

'ОПРЕДЕЛЯЕТСЯ НОВОЕ КОЛ-ВО МЫШЕЙ

m = m + dm

'ЕСТЕСТВЕНАЯ СМЕРТНОСТЬ КОШЕК

dc = 0

IF c > 120 THEN dc = dc - СINT(c / 120) ELSE dc = dc - CINT(RND(1))

IF c < 0 THEN c = 0

'ПРИРОСТ КОЛИЧЕСТВА КОШЕК

SELECT CASE mt

CASE 3, 9

GOSUB 1000

END SELECT

'ПРОДАЖА КОШЕК ЗА НЕНУЖНОСТЬЮ

IF m > mn THEN GOTO 80

60 dc = dc - c + mnc

70 ds = ds + pc * (c - mnc): GOTO 85

80 dc = dc + INT(m \ 30)

ds = ds - INT(m \ 30) * pc

85 s = s + ds

c = c + dc

'В СЛУЧАЕ ГИБЕЛИ ВСЕХ МЫШЕЙ ПРИХОДИТ ГРУППА МЫШЕЙ С ПОЛЯ

IF INT(m) <= 0 THEN m = INT(RND(1) * 20 + 1)

'ОПРЕДЕЛЕНИЕ КОЛИЧЕСТВА ПШЕНИЦЫ

w = w + dw

IF INT(w) < 0 THEN w = 0

IF INT(w) <= 0 THEN GOTO 3000

'ИНФОРМАЦИЯ,ВЫВОДИМАЯ НА ЭКРАН

PRINT "y="; by + y; " ";

PRINT "mt="; mt; " ";

PRINT "w="; INT(w); " ";

PRINT "dw="; INT(dw); " ";

PRINT "m="; INT(m); " ";

PRINT "dm="; INT(dm); " ";

PRINT "c="; INT(c); " ";

PRINT "dc="; INT(dc); " ";

PRINT "ds="; ds

'СОБЫТИЯ ДЕКАБРЯ

IF mt = 12 THEN GOSUB 2500

GOTO 50

'ПРИРОСТ КОШЕК

1000 SELECT CASE INT(m)

CASE IS > 60 * c

dc = dc + 3 * c

CASE IS < 20

dc = dc + 0

CASE 0

dc = dc - .8 * c

CASE ELSE

dc = dc + СINT(RND(1)) + 3

END SELECT

RETURN

'СОБЫТИЯ ДЕКАБРЯ

2500 PRINT "*********************************************************"

'ОПРЕДЕЛЕНИЕ ГОДОВЫХ РАСХОДОВ НА КОШЕК

PRINT "Годовой доход от продажи кошек="; INT(s); "$ "

'РЕАЛИЗАЦИЯ ЗЕРНА С УЧЁТОМ ГОДОВЫХ РАСХОДОВ

INPUT "Сколько кг пшеницы продать"; ws

'РАСЧЁТ ПРИБЫЛИ ОТ РЕАЛИЗАЦИИ ЗЕРНА И ИЗЛИШКА КОШЕК

dw = dw - ws

w = w + dw

ds = ws * pw

s = s + ds

st = st + s

'ГОДОВОЙ БАЛАНС

PRINT "Годовой доход с продажи зерна="; pw * ws; "$ "

PRINT "Общий годовой доход="; INT(s); "$"

PRINT "Общая накопленная сумма="; INT(st); "$"

IF st < 0 THEN GOTO 3000

RETURN

3000 END

5.Тестирование программы.

В связи с тем , что программу протестировать очень сложно из – за большого количества взаимозависмых параметров , при те-

стировании , я использую реакцию программы на исключительные

значения. За год начала прогнозирования я беру 1999 г. Длитель-

ность прогнозирования определяю как равную двум годам. Для того , чтобы прекратить работу программы на прогнозировании второго года , необходимо соответствующим образом задать зна-

чения переменных BW , H и T. Например , чтобы запасы зерна ис-

сякли на втором году прогноза , необходимо , чтобы W равнялось

16000 кг , H также равнялось 16 тоннам , а ежемесячно забиралось

2 тонны пшеницы ( т.е. Т = 2000 кг ).

И действительно программа нашла что , при “содействии” мышей

примерно в 250 кг съеденной за год пшеницы и при затратах на кошек ещё в 20 кг проданной , запасы зерна закончились в марте

второго года прогнозирования.Разумеется , вред , наносимый мышами не учитывать невозможно , но и учесть трудно , т.к. при определении изначального и переопределении нулевого уровня численности мышиной экогруппы используется генератор псевдослучайных чисел и циклическая инициация генератора – RND и RANDOMIZE соответственно. Но приблизительно срок должен совпадать.Кроме того , зная данные за какой – либо ме- сяц , можно по формулам подсчитать значения соответствующих параметров в следующем месяце. Допустим , что MNC=3 ,М в первом месяце прогноза = 27 , а BW , H и T взяты из предыдущего примера. Предельно терпимое количество мышей я взял равным 5.

За месяц меняется только DM , а М остаётся постоянным.Поэтому

 в начале месяца DM = 0 , M = 27 . Далее рассчитывается естест- венная убыль мышей DM = DM – INT ( M/12 ) , т.е. DM = = 0 – INT ( 27/12 ) = 0 – 2 = -2. Потом считается естественный прирост DM = DM + INT ( 1.5*M ) = -2 + INT ( 40.5 ) = -2 + 40 =

= 38. В последнюю очередь считается количество уничтоженных

за месяц мышей : DM = DM – 15*C = 38 – 45 = - 7.

Именно эти данные при тестовом запуске выдала программа.

Таким образом , программа действительно годится для прогнози-

рования и изучения состояния системы взаимовлияющих матема-

тических параметров.