function [out1,out2]=competitive_system(P_0,Q_0,a,b,m,n,length)

out1=zeros(1,length+1);
out2=zeros(1,length+1);
out1(1)=P_0;
out2(1)=Q_0;

for i=1:length
    out1(i+1)=out1(i)+a*out1(i)-b*out1(i)*out2(i);
    out2(i+1)=out2(i)-n*out2(i)+m*out1(i)*out2(i);
end

    