% script file for solving concentrationsODE as a differential equation
clear all;
global k1 k2 k3 k4 km d1 d2 n; %Declare all parameters global

k1 = 0.5; %requests input for k1
k2 = 0.5;
k3 = 1.0; %requests input for k1
k4 = 1.0;
km = 0.5;
d1 = 0.1;
d2 = 0.1;
n = 2;
initIPTG = input('initial concentration of IPTG ');
time = 1000;
[t,y] = ode45(@concentrationsODE_IPTG,[0 time],[initIPTG;10;0;0]);
plot(t,y(:,1),t,y(:,2),t,y(:,3),t,y(:,4))
legend('[IPTG]','[LacI]','[IPTG-LacI complex]','[GFP]','location','Best')
title('Inducible GFP expression ')
xlabel('time')
ylabel('amount')