CLS
PRINT "Laboratorij MRA"
PRINT "Prehrambeno biotehnolo{ki fakultet"
PRINT "----------------------------------"
PRINT "VJE@BA: MODEL STACIONARNIH I DINAMI^KIH STANJA pH"
PRINT : PRINT :
PRINT "1 ODRE\IVANJE KRIVULJE NEUTRALIZACIJE ( stacionarna karakteristika)"
PRINT "2 DINAMI^KE PROMJENE pH "
5 PRINT : PRINT "Odaberite 1 ili 2 ": INPUT r$
IF r$ = "1" GOTO 100
IF r$ = "2" GOTO 200
PRINT "Upi{ite 1 ili 2": GOTO 5
100 CLS : PRINT "STACIONARNA STANJA pH"
PRINT "ODRE\IVANJE KRIVULJE NEUTRALIZACIJE"
PRINT "-----------------------------------"
PRINT
PRINT "Krivulju neutralizacije odredite tako da izaberete koncentraciju "
PRINT "CH3COOH (mol/l) koja je konstantna, a mijenja se koncentraciju NaOH."
PRINT "Koncentracija NaOH se uvecava za konstani iznos, po~ev{i od 0, i"
PRINT "Iterativnim postupkom se izra~unavaju vrijednosti pH u stacionarnim stanjima"
PRINT : PRINT : PRINT "IZABERITE KONCENTRACIJU CH3COOH (mol/l)  "
INPUT y
IF y <= 0 THEN GOTO 100
xmax = 5 * y
dx = y / 10
x = 0
110 REM
IF x > xmax THEN CLS : PRINT "ZAVR[EN PRORA^UN KRIVULJE NEUTRALIZACIJE": END
CLS
PRINT "U PROCESU JE STALNA KONCENTRACIJA CH3COOH (mol/l)="; y
PRINT "KONCENTRACIJA NaOH (mol/l) ="; x
PRINT "--------------------------------------------"

Kv = 1E-14
Ka = .0000185
DEF fnF (pH) = 10 ^ (-2 * pH) + (x + Ka) * 10 ^ (-pH) + ((x - y) * Ka - Kv) - Kv * Ka * 10 ^ pH
A = 0: B = 16
IF fnF(A) * fnF(B) > 0 THEN PRINT "pH>16": END
i = 1
10 F1 = fnF(A)
20 xp = (A + B) / 2
30 F2 = fnF(xp)
PRINT "iteracija=", i, "pH="; xp
IF F2 = 0 THEN 40
IF i > 17 THEN 40
IF F1 * F2 < 0 THEN B = xp: i = i + 1: GOTO 20
A = xp: i = i + 1: GOTO 10
40 PRINT "------------------------"
PRINT "Stacionarna vrijednost pH u procesu ="; xp
IF xp >= 14 THEN END
PRINT
INPUT "NOVO STACIONARNO STANJE (y/n)"; r$: IF r$ = "n" THEN CLS : END
x = x + dx
GOTO 110
END
200 CLS
PRINT "DINAMI^KE PROMJENE pH U PROCESU"
PRINT "---------------------------------------"
PRINT "ZADANI SU SLIJEDE]I PODACI "
PRINT
PRINT "VOLUMEN                  V = 10   l "
PRINT "VOLUMNI PRITOK CH3COOH   Q = 0.75 l/h "
PRINT "VOLUMNI PRITOK NaOH      Q = 0.25 l/h "
PRINT "VOLUMNI PROTOK SIROVINE  Q = 1    l/h "
PRINT "BRZINA  RAZRJE\ENJA JE    D = 0.2  1/h "
D0 = .2: D1 = .075: D2 = .025
PRINT
201 PRINT "IZABERITE SLIJEDE]E PROCESNE PARAMETRE "
PRINT "--------------------------------------"
INPUT "PO^ETNO STANJE KONCENTRACIJE CH3COOH (mol/l) U PROCESU"; y0
INPUT "PO^ETNO STANJE KONCENTRACIJE NaOH    (mol/l) U PROCESU"; x0
INPUT "KONCENTRACIJA CH3COOH (mol/l) U PRITOKU JE"; ypp
INPUT "KONCENTRACIJA NaOH    (mol/l) U PRITOKU JE"; xpp
IF y0 < 0 OR x0 < 0 OR yp < 0 OR xp < 0 THEN CLS : PRINT "koncentracije moraju biti >0": GOTO 201
CLS
PRINT "PO^ETNO STANJE KONCENTRACIJE CH3COOH (mol/l) U PROCESU"; y0
PRINT "PO^ETNO STANJE KONCENTRACIJE NaOH    (mol/l) U PROCESU"; x0
PRINT "KONCENTRACIJA CH3COOH (mol/l) U PRITOKU JE"; ypp
PRINT "KONCENTRACIJA NaOH    (mol/l) U PRITOKU JE"; xpp
PRINT : INPUT "DA LI SU ISPRAVNO UPISANI PODACI (y/n)"; r$
IF r$ = "n" OR r$ = "N" GOTO 300


CLS
t = 0: dt = 1
300 PRINT "vrijeme t(h)                           ="; t
        y = y0 * EXP(-D0 * t) + D1 / D0 * ypp * (1 - EXP(-D0 * t))
x = x0 * EXP(-D0 * t) + D2 / D0 * xpp * (1 - EXP(-D0 * t))
Kv = 1E-14
Ka = .0000185
A = 0: B = 16
IF fnF(A) * fnF(B) > 0 THEN PRINT "pH>16": GOTO 540
i = 1
510 F1 = fnF(A)
520 xp = (A + B) / 2
530 F2 = fnF(xp)
IF F2 = 0 THEN 540
IF i > 20 THEN 540
IF F1 * F2 < 0 THEN B = xp: i = i + 1: GOTO 520
A = xp: i = i + 1: GOTO 510
540 REM
PRINT "pH u procesu                           ="; xp
PRINT "Koncentracija Na+ (mol/l)              ="; x
PRINT "Koncentracija CH3COOH i CH3COO- (mol/l)="; y
PRINT "------------------------"
INPUT "Novi podatak "; r$
IF t > 15 THEN END
t = t + dt
GOTO 300

END