REM  VIVERO04.BAS     

SCREEN 12, 0: COLOR 8: REM tela 640 X 480 texto cinza (color 10)
LINE (20, 0)-(339, 280), 3, B:  REM LINE (0, 0)-(319, 180), 3, B
LINE (20, 100)-(339, 100), 3:   REM LINE (0, 50)-(319, 50), 3
LINE (20, 220)-(339, 220), 3:   REM LINE (0, 120)-(319, 120), 3

J = .004: REM J = 0.005
Q1 = .05
Q2 = 100
Q3 = 100: REM 200
F = 0
LIM1 = 100: REM 1900
LIM2 = 360: REM 2335

K0 = 1: REM 113:
K1 = .2: REM .268:
K2 = .005:
K3 = .000003: REM .0000016:
K4 = .03: REM .03:
K5 = .000268: REM .000268:
K6 = .00165: REM 0.0016
K7 = .003: REM .005:
K8 = 1.34E-06:
K9 = .0012#
K10 = .001432#
K11 = .00002

N0 = 240:   REM 20
Q0 = .045: REM .02
P0 = .2: REM .02
I0 = .01
ANOS = .5
DIAS = 365 * ANOS
T0 = 320 / DIAS
DT = .1

150
I = 3200 + 800 * SIN(T * .017)
R = I / (1 + K0 * Q1)
D1 = J + (K8 * Q2) + (K3 * Q3) - (K4 * Q1) - (K8 * Q1 * R): REM NUTRIENTES
D2 = F + (K1 * Q1 * R) - (K2 * Q2) - (K7 * Q2 * Q3): REM BIOMASSA
REM PEIXES
IF (Q3 < LIM1) THEN
                    D3 = (K11 * Q2 * Q3) - (K6 * Q3)
ELSEIF (Q3 < LIM2) THEN
                    D3 = (K11 * Q2 * Q3) - (K6 * Q3) - (K9 * Q3)
ELSE
                    D3 = (K11 * Q2) - (K6 * Q3) - (K10 * Q3 * Q3)
END IF
Q1 = Q1 + D1 * DT
IF Q1 < .0000001 THEN Q1 = .0000001
Q2 = Q2 + D2 * DT
Q3 = Q3 + D3 * DT
T = T + DT

REM Y = INT(T / 365): W = 12 * Y: M = INT(T / 30.4) - W
PSET (20 + T0 * T, 100 - I * I0), 14: REM I ENERGIA SOLAR
PSET (20 + T0 * T, 220 - Q2 * Q0), 2: REM Q2 BIOMASSA
PSET (20 + T0 * T, 220 - Q3 * P0), 6: REM Q3 PEIXE
PSET (20 + T0 * T, 280 - Q1 * N0), 5: REM Q1 NUTRIENTES
LOCATE 24, 1: PRINT Q1; Q2; Q3:

IF T * T0 < 319 GOTO 150
T = 0: Q1 = .05: Q2 = 100: Q3 = 100:
F = F + 200!
IF F < 600! GOTO 150
END