%script file for solving gfp.m as a differential equation
clear all;
global k1 d1; % declare k1 and d1 global

k1 = input('k1 '); %requests input for k1
d1 = input('d1 '); %requests input for d1
time = input('timepoints '); %requests number of time points
initGFP = input('initial amount of GFP '); %requests initial condition

[t,y] = ode45(@gfp,[0 time],initGFP);