function dxdt = ODE_resources(t,x);

global a alphdot m;

x1 = x(1)   %[volume]
x2 = x(2)   %[x = Resource/Km]

dx1dt = a*((x2^m) / (1 + (x2^m)))*x1;
dx2dt = -(((alphdot)*a)*((x2^m) / (1 + (x2^m))))* x1;

dxdt = [dx1dt; dx2dt];