10 REM IBM 20 REM CLIMAX 30 CLS 40 SCREEN 1, 0: COLOR 15, 0 50 LINE (0, 0)-(319, 180), 3, B 60 REM COEFFICIENTS 70 I = 1 80 S = 1 90 B = 100 100 D = 1 110 K0 = .027 120 K1 = .0002 130 K2 = 2 140 K3 = .00417 150 K4 = .000025 160 K5 = .01 170 K6 = .000175 180 K7 = .00004 190 K8 = .000005 200 K9 = .002 205 TX = 5400 210 T0 = 360 / TX 220 B0 = .003 230 G0 = 3 240 P0 = 3.2 250 D0 = 3 260 DT = 10 280 PSET (T * T0, 180 - D * D0), 1 290 PSET (T * T0, 180 - B * B0), 2 300 PSET (T * T0, 180 - PG * P0), 3 305 PSET (T * T0, 180 - DB * G0), 3 310 R = I / (1 + K0 * D + K1 * B) 320 DD = K8 * B * S - K9 * D - K7 * D * D 330 PG = K2 * R * D + K3 * R * B: REM gross production 340 DB = PG - K6 * B - K4 * B - K5 * D * D 350 B = B + DB * DT 360 D = D + DD * DT 370 T = T + DT 380 IF T * T0 < 319 GOTO 280 385 LOCATE 2, 2: PRINT "Climax.bas" 390 SLEEP: SLEEP 400 END