% population dynamics (wolfs reindeer) 

a = 8;
b = 0.25;
c = 0.02;
d = 20;
timint = [0;2];
init   = [600;60];

[t,y] = ode45(@Fun_PPdynamics,timint,init,[],a,b,c,d);

plot(t,y,'-*')
