clear all
global r1 rn1 r3 r2 rn2 

t0 = 0;
t1 = 10;
a0 = 2; 
b0 = 0.9;
c0 = 4.2;


r1 = 1;
rn1 = 1;
r3 = 4; 
r2 = 1; 
rn2 = 1;  

S0  = [a0;b0;c0];

[t,St] = ode45('GomesMorrisPresentationfn',[t0,t1],S0);
    
plot(t,St(:,1), t, St(:,2), t, St(:,3)); 
title('Nitrogen Metabolism Model');
xlabel('Time');
ylabel('Concentration');

legend('alpha-ketoglutarate', 'glutamate', 'glutamine');