function dAdt = induced(t,y) % Rate of change of GFP with inducer
global beta km n d2 inducer %declare parameters as global
dAdt = (beta*inducer^n)/(km^n+inducer^n) - d2*y; %equation