% Fun Population dynamics

function y = Fun_PPdynamics(t,x,a,b,c,d)
y(1,1)= x(1,1)*(a-b*x(2,1) );  
y(2,1) =x(2,1)*(c*x(1,1)-d );  


%= () 
