Modelo Depredador-Presa por Juan José Moreno Balcázar

Caso Lince - Conejo

a=Input["Introduce el número de conejos"];
b=Input["Introduce el número de linces"];
n=Input["Introduce el número de años"];

sol=NDSolve[{y'[t]==0.6*y[t]-0.018*y[t]*z[t],z'[t]==-0.7*z[t]+0.012*y[t]*z[t], y[0]==a,z[0]==b},{y,z},{t,0,n}];

lista=Table[i*0.001,{i,0,Floor[(n/0.001)]}];
y1=Flatten[y[lista]/.sol];
p1=Transpose[{lista,y1}];
z1=Flatten[z[lista]/.sol];
p2=Transpose[{lista,z1}];
dib[1]=ListPlot[p1, PlotStyle->{RGBColor[0,0,1]},DisplayFunction->Identity];
dib[2]=ListPlot[p2, PlotStyle->{RGBColor[0,1,0]},DisplayFunction->Identity];
Show[dib[1],dib[2],DisplayFunction->$DisplayFunction]

La curva azul representa la población de conejos y la verde la de linces, ambas en función del tiempo.

puntos=Transpose[{y1,z1}];
ListPlot[puntos, PlotStyle->{RGBColor[1,0,0]}, AxesLabel->{Conejos, Linces}]