%Vector c with a and b
c = [a b];

%Puts function into a vector format
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];

%Solves using ODE45
[jared,baller] = ode45 (dxdt,[0 200],[0 1]);

hold all;

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

%Solves using RK4
[why, amidoingthis] = RK4(dxdt,[0 200],[0 1],200);

%Plots RK4
plot(why, amidoingthis);

%Legen....wait for it...DARY
legend('Ode45 A','Ode45 B','RK4 A','RK4 B');