function dAdt = GFP(t,y) % Rate of change of GFP
global k1 d1 %declare k1 and d1 as global
dAdt = k1 - d1*y; %equation