'STACIONARNA STANJA ENZIMSKE REAKCIJE
'METODA NEWTON-RAPHSONA
CLS
10 REM
PRINT "Stacionarna stanja enzimske reakcije "
PRINT "Metoda Newton-Raphsona"
INPUT "Upi{ite naziv enzima: enzim= "; enzim$
INPUT "Upi{ite naziv supstrata: supstrat="; supstrat$
INPUT "Upi{ite vrijednost temperature: T(C)="; temp
INPUT "Upi{ite vrijednost pH ="; pH
INPUT "Upi{ite koncentraciju enzima Eo (mol/l) ="; Eo
INPUT "Upi{ite koncentraciju supstrata u ulaznom toku So(mol/l)="; So
INPUT "Upi{ite vrijednost kineti~ke konstante k2 (1/s)="; k2
INPUT "Upi{ite vrijednost konstante 1/Km (l/mol) ="; kminv
INPUT "Upi{ite brzinu razrje|enja D(1/s)="; D
CLS
PRINT " ZADANI PODACI"
PRINT "---------------"
PRINT "enzim="; enzim$
PRINT "supstrat="; supstrat$
PRINT "temperatura (C)="; temp
PRINT "pH ="; pH
PRINT "koncentracija enzima Eo(mol/l)="; Eo
PRINT "koncentracija supstrata u ulaznom toku So (mol/l)="; So
PRINT "kineti~ka konstanta k2(1/s)="; k2
PRINT "recipro~na Michaelis-Mentenova konstanta 1/Km(l/mol)="; kminv
PRINT "brzina razre|enja D(1/s)="; D
PRINT
INPUT "Da li su podaci dobro upisani (da/ne) ="; r$
IF r$ = "ne" THEN PRINT : PRINT "*******  PONOVITE UPIS  ***** ": GOTO 10
km = 1 / kminv
vmax = k2 * Eo
CLS
PRINT "Michaelis-Mentenova konstanta  Km (mol/l)="; km
PRINT "Maksimalna brzina enzimske reakcije Vmax (1/s )="; vmax
INPUT "ZADAJTE APSOLUTNU VRIJEDNOST DOZVOLJENE POGRE[KE eps(mol/l)="; eps
INPUT "Zadajte pribli`nu vrijednost stacionarnog stanja "; s
p = km + vmax / D - So
segz = -p / 2 + SQR(p * p / 4 + So * km)
k = 1
PRINT
PRINT "-----------------------------------------"
PRINT "po~etak NEWTON-RAPHSONOVE iteracije"
PRINT "-----------------------------------------"
20 e = s - segz
PRINT
PRINT "Redni broj iteracije k="; k
PRINT "Koncentracija supstrata s(mol/l)="; s
PRINT "Pogre{ka e="; e
INPUT "Nova iteracija (*/*)"; r$
IF ABS(e) < eps THEN 30
k = k + 1
f = D * (So - s) - vmax / (km + s)
df = -D - vmax * km / (km + s) ^ 2
s = s - f / df
GOTO 20
30 PRINT "------------------"
CLS
PRINT "  R J E [ E N J E "
PRINT "*******************"
es = Eo * s / (km + s)
ee = Eo - es
p = So - s
PRINT "enzim="; enzim$
PRINT "supstrat="; supstrat$
PRINT "temperatura (C)="; temp
PRINT "pH ="; pH
PRINT "koncentracija enzima Eo(mol/l)="; Eo
PRINT "koncentracija supstrata u ulaznom toku So="; So
PRINT "maksimalna brzina vmax(1/s)="; vmax
PRINT "Michaelis-Mentenova konstanta Km(mol/l)="; km
PRINT "brzina razre|enja D(1/s)="; D
PRINT "---------------------------------------"
PRINT "STACIONARNA STANJA "
PRINT "KONCENTRACIJA SUPSTRATA  s(mol/l)="; s
PRINT "KONCENTRACIJA ME\UKOMPLEKSA ENZIM-SUPSTRAT ES(mol/l)="; es
PRINT "KONCENTRACIJA SLOBODNOG ENZIMA E(mol/l) ="; ee
PRINT "KONCENTRACIJA PRODUKTA p(mol/l) ="; p
PRINT "POGRE[KA U PRORA^UNU KONCENTRACIJE SUPSTRATA e(mol/l)="; s - segz
END