[t,y] = ode45(@predprey,[0 25],[20;20]);
plot(t,y(:,1),'b-o',t,y(:,2),'r--o')
title('Populations of foxes (red) and bunnies(blue)')
xlabel('time')
ylabel('population')
