c = [a b];
dxdt = @(t, c)[0.03 - 0.02*c(1) - c(1)*c(2);c(1)*c(2) - (0.11*c(2))/(c(2) + 0.08) + 0.0005];

[jared,baller] = ode45 (dxdt,[0 200],[0 1]);

hold all;

plot(jared,baller(:,1),'r');
plot(jared,baller(:,2),'g');