clear all
clc
timerange=[0,20];
initialValues=[100,90,80];
[t,s]=ode45(@testfunction,timerange,initialValues);
grid on
plot(s(:,1),s(:,2))
ylabel('rabbits')
xlabel('foxes')