Как построить 2 графика на одном рисунке?

x = 0:0.5:10;
y1 = cos(x) - 2;
y2 = sin(x);
plot(x,y1,'b');
hold on;
plot(x,y2,'r');
hold off;

x = 0:0.5:10;
y1 = cos(x) - 2;
y2 = sin(x);
subplot(1,2,1); plot(x,y1,'b');
subplot(1,2,2); plot(x,y2,'r');

S=50;
Talone=[1.5 2 2 2 2 2 2 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 3 3 3 3 3 3 3 3 3 3 3.5 3.5 3.5 3.5 3.5 3.5 3.5 3.5 3.5 4 4 4 4 4 4 4.5 4.5 4.5 4.5 5 5 5.5 5.5 6];
Talone=Talone*1000;
Tmin_fix=0.9*1000;
Tmin_min=0.6*1000;
Tmin_max=1.2*1000;
shagTmin=0.1*1000;
kolShagT=(Tmin_max-Tmin_min)/shagTmin+1;
C=((1/S)*sum(1./Talone))^(-1);
A=1500;
Ro=A/(2*C);
maxM=10;
minM=4;
mFix=5;
m_plot=zeros(1,(maxM-minM)+1);
p_plot=zeros(1,(maxM-minM)+1);
for countM=1 : 1 : (maxM-minM)+1;
m=countM+minM-1;
m_plot(countM)=m;
t=Talone./m;
Pr1=(Ro^m)*(1-Ro);
Em=Ro/(1-Ro);
Pr2=(m*Pr1)/Em;
T=((1-Ro)/S)*sum(Talone);
Tmin=Tmin_fix;
k=(Talone)./Tmin;
k1=floor(k);
a=zeros(1,length(k1));
a1=zeros(1,length(k1));
for j = 1 : 1: S;
if k1(j)>=m 
    a(j)=(m*(Ro)^(m-1))*(1-Ro)^2;
else
    continue;
end;
if j>1
    if k1(j)~=k1(j-1)
        a1(j)=sum(a);
    else
        a1(j)=a1(j-1);
    end;
else
    a1(1)=a(1);
end;
end;
Pr1Tmin=sum(a1)/S;
p_plot(countM)=Pr1Tmin;
end;
plot(m_plot,p_plot,'-r','LineWidth',2);
hold on
% pause;
t_plot=zeros(1,kolShagT);
p_plot=zeros(1,kolShagT);
for countT=1 : 1 : kolShagT;
m=mFix;
t=Talone./m;
Pr1=(Ro^m)*(1-Ro);
Em=Ro/(1-Ro);
Pr2=(m*Pr1)/Em;
T=((1-Ro)/S)*sum(Talone);
Tmin=Tmin_min+(countT-1)*shagTmin;
t_plot(countT)=Tmin;
k=(Talone)./Tmin;
k1=floor(k);
a=zeros(1,length(k1));
a1=zeros(1,length(k1));
for j = 1 : 1: S;
if k1(j)>=m 
    a(j)=(m*(Ro)^(m-1))*(1-Ro)^2;
else
    continue;
end;
if j>1
    if k1(j)~=k1(j-1)
        a1(j)=sum(a);
    else
        a1(j)=a1(j-1);
    end;
else
    a1(1)=a(1);
end;
end;
Pr1Tmin=sum(a1)/S;
p_plot(countT)=Pr1Tmin;
end;
plot(t_plot,p_plot,'LineWidth',2);
 
% случай с двойной несущей домножаю Talone на 2 в сторке 26 и 62  
% уменьшаю спорс на трафик А в два раза(строка) 10

clc;
clear all;
step=0.01;
x1=-1:step:0;
y1=(2+sin(x1).^2)/(1 + x1.^2);
x2=0:step:2;
y2=(4*cos(3*x2))/(1+exp(3*x2));
x=[x1 x2];
y=[y1 y2];
plot(x,y,'r')
hold on;
xlabel ('X')
ylabel ('Y')
grid on;
xlim([-10 10]);
ylim([0 10]);
title('График');
set(gca,'XTick',-pi:pi/1:pi);