clear
tspan = [0 10];% s
y0=[2 1 0];% uM
[t,y] = ode45('hw3_1',tspan,y0);
plot(t,y)
xlabel('TIME (s)')
ylabel('CONCENTRATION (µM)')
axis([0 tspan(2) 0 max(y0)])