global CGS DGS CGO DGO kfGLU krGLU kfNH3 krNH3 kfGLN krGLN kfALPK krALPK kcatGS kcatGO Con1 Con2;

clc;
clear all;
close all;

CGS=0.01; %%find a realistic value
DGS=0.01;  %%find a realistic value
CGO=0.01; %%find a realistic value
DGO=0.01;  %%find a realistic value


kfGLU=2; %Calculated from Km- set
krGLU=13; %Calculated in proportion to set KfGLU
kfNH3=15; %Set
krNH3=8; %Parametrized from Kf
%KmALPK = 0.01
kfGLN=2; %CHANGE
krGLN=13; %CHANGE
%KmGLN=0.1
kfALPK=1; %CHANGE
krALPK=8; %CHANGE

Con1= 4; %%find a realistic value
Con2= 40;%%find a realistic value. 

kcatGS=43;
kcatGO= 40;  %turnover from brenda azospirillum


[T,Y]= ode15s('MichaelisGSGO3', [0 100],[1 1 1 1 1 1 1 1]);
%plot(T,Y);
%legend('GS','GLU','GSGLU','GLN','GO','GOGLN','GLNGOALPK','GLUGSNH3');