clear
R(1)=100.0   %initial population of rabbits
F(1)=40.0    %initial population of foxs
BR_rabbit=0.5     %birth rate of rabbit 
DR_rabbit_INT=0.02   %death date of rabbit (prey) due to interaction
DR_fox=0.10     %death rate of fox
BR_fox_INT= 0.001 % birth rate of fox (predator) due to interaction

for i=1:50
 R(i+1)=R(i)+BR_rabbit*R(i)-DR_rabbit_INT*R(i)*F(i)
 F(i+1)=F(i)-DR_fox*F(i)+BR_fox_INT*F(i)*R(i)
end

hold off  %clear the previous plot setting
figure(1)  %make plot 1 for rabbit population evolution
plot(R)   
figure(2)  %make plot 2 in another window for fox population evolution
plot(F)
