clear
%y-x
for i=0.1:0.2:0.9
    for j=0.1:0.2:0.9
        [T,Y]=ode45('differential1',[0 20],[i j]);
        figure(1)
        grid on
        plot(T,Y(:,1),'k*-','lineWidth',1);
        hold on
        plot(T,Y(:,2),'k--','lineWidth',1);
        hold on
    end
end
axis([0 20 -0.1 1.1]);
set(gca,'XTick',[0:2:20], 'YTick',[-0.1:0.1:1.1]);
xlabel('evolutionary steps');
ylabel('population');
legend('OSN service providers','Online platforms')


clear
%y-x
for i=0.1:0.2:0.9
    for j=0.1:0.2:0.9
        [T,Y]=ode45('differential2',[0 20],[i j]);
        figure(1)
        grid on
        plot(T,Y(:,1),'k*-','lineWidth',1);
        hold on
        plot(T,Y(:,2),'k--','lineWidth',1);
        hold on
    end
end
axis([0 20 -0.1 1.1]);
set(gca,'XTick',[0:2:20], 'YTick',[-0.1:0.1:1.1]);
xlabel('evolutionary steps');
ylabel('population');
legend('OSN service providers','Online platforms')

clear
%y-x
for i=0.1:0.2:0.9
    for j=0.1:0.2:0.9
        [T,Y]=ode45('differential3',[0 20],[i j]);
        figure(1)
        grid on
       plot(T,Y(:,1),'k*-','lineWidth',1);
        hold on
        plot(T,Y(:,2),'k--','lineWidth',1);
        hold on
    end
end
axis([0 20 -0.1 1.1]);
set(gca,'XTick',[0:2:20], 'YTick',[-0.1:0.1:1.1]);
xlabel('evolutionary steps');
ylabel('population');
legend('OSN service providers','Online platforms')

clear
%y-x
for i=0.1:0.2:0.9
    for j=0.1:0.2:0.9
        [T,Y]=ode45('differential4',[0 20],[i j]);
        figure(1)
        grid on
        plot(T,Y(:,1),'k*-','lineWidth',1);
        hold on
        plot(T,Y(:,2),'k--','lineWidth',1);
        hold on
    end
end
axis([0 20 -0.1 1.1]);
set(gca,'XTick',[0:2:20], 'YTick',[-0.1:0.1:1.1]);
xlabel('evolutionary steps');
ylabel('population');
legend('OSN service providers','Online platforms')

clear
%y-x
for i=0.1:0.2:0.9
    for j=0.1:0.2:0.9
        [T,Y]=ode45('differential6',[0 20],[i j]);
        figure(1)
        grid on
        plot(T,Y(:,1),'k*-','lineWidth',1);
        hold on
        plot(T,Y(:,2),'k-','lineWidth',1);
        hold on
    end
end
axis([0 20 -0.1 1.1]);
set(gca,'XTick',[0:2:20], 'YTick',[-0.1:0.1:1.1]);
xlabel('evolutionary steps');
ylabel('population');
legend('OSN service providers','Online platforms')

clear
%y-x
for i=0.1:0.2:0.9
    for j=0.1:0.2:0.9
        [T,Y]=ode45('freeriding',[0 20],[i j]);
        figure(1)
        grid on
       plot(T,Y(:,1),'k*-','lineWidth',1);
        hold on
        plot(T,Y(:,2),'k--','lineWidth',1);
        hold on
    end
end
axis([0 20 -0.1 1.1]);
set(gca,'XTick',[0:2:20], 'YTick',[-0.1:0.1:1.1]);
xlabel('evolutionary steps');
ylabel('population');
legend('OSN service providers','Online platforms')