Multiple curves in a single graph

In a program i am plotting graph between two variables.

By changing Value of ''AOM2'' to 0,10,25,50 ,I will get different plots. I want to draw all those curves in a single figure.

I already tried many things like ''hold on'' and other commands also but nothing is working.

if true
    S1=[];
ADD=[];
AD=linspace(0.00001,1,50);
for AD1=1:size(AD,2);
    AD2=AD(1,AD1);
    ADD=[ADD:AD2];

       AOM2=25;%(ohm)

       AOG=5;                 %omega
       AOMG=AOM2./AOG;   %OHM UPON OMEGA
       AOMN=1+ (AOMG).^2 ; %1+OHM SQUARE UPON OMEGA SQUARE
       C4S=2*1i.*AOMG;
       ATT=27;%(T0=300K)
      ALM=3.17;% LAMBDA
      AMU=1.639;%MU
      RO=1.74; %RHO
      AT2=.005; %TAU NOT
      ACE=23.61; %C_E
      AK=0.492; %K
       ALF1=0.05;%ALPHA_T
      GAMA=1.13849;
      AD0=1.475; %d_0
      AXI=2; %axie
      AZT0=3.688;%ZETA NOT
      ZTA=.01;%ZETA 
      AJ=1.753;%J
      AT3=.05;%TAU
       AKS=1;%(small k wave number)
      AB=0.05;%a
      OGB=0.05;%b
       AT=0.05;%TAU 
      ALF2=0.05;%ALPHA_C=0.05
     ET=1i;

      ATMS=AT2+(ET./AOG);% (tau power m)
      ATNS=AT3+(ET./AOG);%(tau power n)
      C1S=(ALM+2.*AMU)/RO;
      C2S=AMU./RO;
      C3S=(AZT0/(RO*AJ))+AD0/(RO*AJ*AKS);%(delta= )
      %C3S=(AZT0/(RO*AJ))+D0/(RO*AJ*AKS);
      %ATMS=AT2*ONE+(ET/OG);
      %ATNS=AT3*ONE+(ET/OG);
      ABT1=(3.*ALM+2.*AMU)*ALF1;%beta 1
      ABT2=(3.*ALM+2.*AMU)*ALF2;%beta 2
      ABN=(ABT1*ATT)/(RO*ACE);%
      ABT1B=ABT1/RO;%beta 1 bar
      ABT2B=ABT2/RO; % beta 2 bar
      AD1=AK./(RO*ACE*ATMS); %d_1
      AD2=AD2*AB./ATNS;  %d_2
      AD3=AD2*ZTA./ATNS;  %d_3
      AD4=GAMA./RO; %d_4
      AD5=AD2*OGB./ATNS; %d_5
      AD6=AD2*ABT2./ATNS; %%d_6
      GAMAB=GAMA./(RO*AJ*AKS);%gamma bar
      AXIB=AXI./(RO*AJ*AKS); % exie bar 
      AZTB=ZTA./(RO*AJ*AKS); % zeta bar
      ABN=(ABT1*ATT)/(RO*ACE);% epsilon
      ABN0=(OGB*ATT)/(RO*ACE); % epsilon not 
      ABN1=(AXI*ATT)/(RO*ACE);%epsilon star

       B0=(-AOMN^2-C4S^2);
      B1=((C1S)*(AOMN) + (AOMN)*(C2S) + (AD1)*(AOMN)^2 + (AD1)*(C4S)^2 + (AOMN)^2*(AD2) + (C4S)^2*(AD2) + (AOMN)^2*(C3S) + (C4S)^2*(C3S) + (C4S)^2*(ABN1)*AXIB + (AOMN)*(ABT1B)*(ABN) + (ABN0)*(AOMN)^2*(AD5) + (ABN0)*(C4S)^2*(AD5) + (AOMN)^2*(ABN1)*AXIB);
      B2=(- (C1S)*(C2S) - (AOMN)^2*(AD2)*(C3S) + (AOMN)^2*(AD3)*(AZTB) - (C4S)^2*(AD2)*(C3S) + (C4S)^2*(AD3)*(AZTB) - (AD1)*(C1S)*(AOMN) - (AD1)*(AOMN)*(C2S) - (C1S)*(AOMN)*(AD2) + (AOMN)*(ABT2B)*(AD6) - (ABT1B)*(C2S)*(ABN) - (AOMN)*(C2S)*(AD2) - (C1S)*(AOMN)*(C3S) + (AOMN)*(AD4)*(GAMAB) - (AOMN)*(C2S)*(C3S) - (AD1)*(AOMN)^2*(AD2) - (AD1)*(C4S)^2*(AD2) - (AD1)*(AOMN)^2*(C3S) - (AD1)*(C4S)^2*(C3S) - (C1S)*(AOMN)*(ABN1)*(AXIB) - (AOMN)*(ABT1B)*(ABN)*(C3S) - (AOMN)*(ABT1B)*(ABN1)*(GAMAB) - (AOMN)*(AD4)*(ABN)*(AXIB) - (AOMN)*(C2S)*(ABN1)*(AXIB) - (ABN0)*(AOMN)^2*(AD5)*(C3S) - (ABN0)*(AOMN)^2*(AD3)*(AXIB) - (ABN0)*(C4S)^2*(AD5)*(C3S) - (ABN0)*(C4S)^2*(AD3)*(AXIB) - (AOMN)^2*(ABN1)*(AD5)*(AZTB) - (AOMN)^2*(ABN1)*(AD2)*(AXIB) - (C4S)^2*(ABN1)*(AD5)*(AZTB) - (C4S)^2*(ABN1)*(AD2)*(AXIB) - (ABN0)*(C1S)*(AOMN)*(AD5) - (ABN0)*(AOMN)*(ABT1B)*(AD6) - (ABN0)*(AOMN)*(C2S)*(AD5) - (AOMN)*(ABT1B)*(ABN)*(AD2) - (AOMN)*(ABT2B)*(ABN)*(AD5));
      B3= ((AD1)*(C1S)*(C2S) + (C1S)*(C2S)*(AD2) - (ABT2B)*(C2S)*(AD6) + (C1S)*(C2S)*(C3S) - (AD4)*(C2S)*(GAMAB) + (C1S)*(AOMN)*(AD2)*(C3S) + (C1S)*(C2S)*(ABN1)*(AXIB) - (AOMN)*(ABT2B)*(AD6)*(C3S) - (AOMN)*(AD4)*(AD6)*(AZTB) - (AOMN)*(AD4)*(AD2)*(GAMAB) + (ABT1B)*(C2S)*(ABN)*(C3S) + (ABT1B)*(C2S)*(ABN1)*(GAMAB) + (AD4)*(C2S)*(ABN)*(AXIB) - (C1S)*(AOMN)*(AD3)*(AZTB) - (AOMN)*(ABT2B)*(AD3)*(GAMAB) + (AOMN)*(C2S)*(AD2)*(C3S) - (AOMN)*(C2S)*(AD3)*(AZTB) + (AD1)*(AOMN)^2*(AD2)*(C3S) - (AD1)*(AOMN)^2*(AD3)*(AZTB) + (AD1)*(C4S)^2*(AD2)*(C3S) - (AD1)*(C4S)^2*(AD3)*(AZTB) + (AD1)*(C1S)*(AOMN)*(AD2) - (AD1)*(AOMN)*(ABT2B)*(AD6) + (ABN0)*(C1S)*(C2S)*(AD5) + (ABN0)*(ABT1B)*(C2S)*(AD6) + (AD1)*(AOMN)*(C2S)*(AD2) + (AD1)*(C1S)*(AOMN)*(C3S) - (AD1)*(AOMN)*(AD4)*(GAMAB) + (AD1)*(AOMN)*(C2S)*(C3S) + (ABT1B)*(C2S)*(ABN)*(AD2) + (ABT2B)*(C2S)*(ABN)*(AD5) + (ABN0)*(C1S)*(AOMN)*(AD5)*(C3S) + (ABN0)*(AOMN)*(ABT1B)*(AD6)*(C3S) + (ABN0)*(AOMN)*(AD4)*(AD6)*(AXIB) - (ABN0)*(AOMN)*(AD4)*(AD5)*(GAMAB) + (ABN0)*(C1S)*(AOMN)*(AD3)*(AXIB) + (ABN0)*(AOMN)*(ABT1B)*(AD3)*(GAMAB) + (ABN0)*(AOMN)*(C2S)*(AD5)*(C3S) + (ABN0)*(AOMN)*(C2S)*(AD3)*(AXIB) + (C1S)*(AOMN)*(ABN1)*(AD5)*(AZTB) + (C1S)*(AOMN)*(ABN1)*(AD2)*(AXIB) + (AOMN)*(ABT1B)*(ABN)*(AD2)*(C3S) + (AOMN)*(ABT1B)*(ABN1)*(AD6)*(AZTB) + (AOMN)*(ABT1B)*(ABN1)*(AD2)*(GAMAB) + (AOMN)*(ABT2B)*(ABN)*(AD5)*(C3S) - (AOMN)*(ABT2B)*(ABN1)*(AD6)*(AXIB) + (AOMN)*(ABT2B)*(ABN1)*(AD5)*(GAMAB) + (AOMN)*(AD4)*(ABN)*(AD5)*(AZTB) + (AOMN)*(AD4)*(ABN)*(AD2)*(AXIB) - (AOMN)*(ABT1B)*(ABN)*(AD3)*(AZTB) + (AOMN)*(ABT2B)*(ABN)*(AD3)*(AXIB) + (AOMN)*(C2S)*(ABN1)*(AD5)*(AZTB) + (AOMN)*(C2S)*(ABN1)*(AD2)*(AXIB));
      B4  = (- (C1S)*(C3S)*(AD2)*(C3S) + (ABT2B)*(C3S)*(AD6)*(C3S) + (AD4)*(C3S)*(AD6)*(AZTB) + (AD4)*(C3S)*(AD2)*(GAMAB) + (C1S)*(C3S)*(AD3)*(AZTB) + (ABT2B)*(C3S)*(AD3)*(GAMAB) - (AD1)*(C1S)*(C3S)*(AD2) + (AD1)*(ABT2B)*(C3S)*(AD6) - (AD1)*(C1S)*(C3S)*(C3S) + (AD1)*(AD4)*(C3S)*(GAMAB) - (AD1)*(C1S)*(AOMN)*(AD2)*(C3S) + (AD1)*(AOMN)*(ABT2B)*(AD6)*(C3S) + (AD1)*(AOMN)*(AD4)*(AD6)*(AZTB) + (AD1)*(AOMN)*(AD4)*(AD2)*(GAMAB) + (AD1)*(C1S)*(AOMN)*(AD3)*(AZTB) + (AD1)*(AOMN)*(ABT2B)*(AD3)*(GAMAB) - (ABN0)*(C1S)*(C3S)*(AD5)*(C3S) - (ABN0)*(ABT1B)*(C3S)*(AD6)*(C3S) - (ABN0)*(AD4)*(C3S)*(AD6)*(AXIB) + (ABN0)*(AD4)*(C3S)*(AD5)*(GAMAB) - (AD1)*(AOMN)*(C3S)*(AD2)*(C3S) - (ABN0)*(C1S)*(C3S)*(AD3)*(AXIB) - (ABN0)*(ABT1B)*(C3S)*(AD3)*(GAMAB) + (AD1)*(AOMN)*(C3S)*(AD3)*(AZTB) - (C1S)*(C3S)*(ABN1)*(AD5)*(AZTB) - (C1S)*(C3S)*(ABN1)*(AD2)*(AXIB) - (ABT1B)*(C3S)*(ABN)*(AD2)*(C3S) - (ABT1B)*(C3S)*(ABN1)*(AD6)*(AZTB) - (ABT1B)*(C3S)*(ABN1)*(AD2)*(GAMAB) - (ABT2B)*(C3S)*(ABN)*(AD5)*(C3S) + (ABT2B)*(C3S)*(ABN1)*(AD6)*(AXIB) - (ABT2B)*(C3S)*(ABN1)*(AD5)*(GAMAB) - (AD4)*(C3S)*(ABN)*(AD5)*(AZTB) - (AD4)*(C3S)*(ABN)*(AD2)*(AXIB) + (ABT1B)*(C3S)*(ABN)*(AD3)*(AZTB) - (ABT2B)*(C3S)*(ABN)*(AD3)*(AXIB));
      B5=((AD1)*(C1S)*(C3S)*(AD2)*(C3S) - (AD1)*(ABT2B)*(C3S)*(AD6)*(C3S) - (AD1)*(AD4)*(C3S)*(AD6)*(AZTB) - (AD1)*(AD4)*(C3S)*(AD2)*(GAMAB) - (AD1)*(C1S)*(C3S)*(AD3)*(AZTB) - (AD1)*(ABT2B)*(C3S)*(AD3)*(GAMAB));

      Y=[B0 B1 B2 B3 B4 B5];
      R=roots(Y);
  S=real(sqrt(R));
S1=[S1 S];

end
figure(1)
plot(AD,S1(5,:));
end

Kindly Help me with this.

Answer 1

[KSSV]

USe :

plot(AD,S1);

instead of:

plot(AD,S1(5,:));

[Himanshu SINGLA]

But it is not giving different curves for different values of Ohm.

[KSSV]

Then you need to check your coding.....

[Himanshu SINGLA]

ok, this is right, It is showing the graph of all the five roots.

But i am asking you that if i put AOM2 different values say 0,10,50,100 then curves will differ. so corresponding to one root(for different values of AOM2) I want curves in one graph.

[KSSV]

O = [0,10,50,100] ;
figure
hold on
for i = 1:length(O)
    S1=[];
    ADD=[];
    AD=linspace(0.00001,1,50);
    for AD1=1:size(AD,2);
        AD2=AD(1,AD1);
        ADD=[ADD:AD2];
        AOM2=O(i);%(ohm)
        AOG=5;                 %omega
        AOMG=AOM2./AOG;   %OHM UPON OMEGA
        AOMN=1+ (AOMG).^2 ; %1+OHM SQUARE UPON OMEGA SQUARE
        C4S=2*1i.*AOMG;
        ATT=27;%(T0=300K)
        ALM=3.17;% LAMBDA
        AMU=1.639;%MU
        RO=1.74; %RHO
        AT2=.005; %TAU NOT
        ACE=23.61; %C_E
        AK=0.492; %K
        ALF1=0.05;%ALPHA_T
        GAMA=1.13849;
        AD0=1.475; %d_0
        AXI=2; %axie
        AZT0=3.688;%ZETA NOT
        ZTA=.01;%ZETA
        AJ=1.753;%J
        AT3=.05;%TAU
        AKS=1;%(small k wave number)
        AB=0.05;%a
        OGB=0.05;%b
        AT=0.05;%TAU
        ALF2=0.05;%ALPHA_C=0.05
        ET=1i;
        ATMS=AT2+(ET./AOG);% (tau power m)
        ATNS=AT3+(ET./AOG);%(tau power n)
        C1S=(ALM+2.*AMU)/RO;
        C2S=AMU./RO;
        C3S=(AZT0/(RO*AJ))+AD0/(RO*AJ*AKS);%(delta= )
        %C3S=(AZT0/(RO*AJ))+D0/(RO*AJ*AKS);
        %ATMS=AT2*ONE+(ET/OG);
        %ATNS=AT3*ONE+(ET/OG);
        ABT1=(3.*ALM+2.*AMU)*ALF1;%beta 1
        ABT2=(3.*ALM+2.*AMU)*ALF2;%beta 2
        ABN=(ABT1*ATT)/(RO*ACE);%
        ABT1B=ABT1/RO;%beta 1 bar
        ABT2B=ABT2/RO; % beta 2 bar
        AD1=AK./(RO*ACE*ATMS); %d_1
        AD2=AD2*AB./ATNS;  %d_2
        AD3=AD2*ZTA./ATNS;  %d_3
        AD4=GAMA./RO; %d_4
        AD5=AD2*OGB./ATNS; %d_5
        AD6=AD2*ABT2./ATNS; %%d_6
        GAMAB=GAMA./(RO*AJ*AKS);%gamma bar
        AXIB=AXI./(RO*AJ*AKS); % exie bar
        AZTB=ZTA./(RO*AJ*AKS); % zeta bar
        ABN=(ABT1*ATT)/(RO*ACE);% epsilon
        ABN0=(OGB*ATT)/(RO*ACE); % epsilon not
        ABN1=(AXI*ATT)/(RO*ACE);%epsilon star
        B0=(-AOMN^2-C4S^2);
        B1=((C1S)*(AOMN) + (AOMN)*(C2S) + (AD1)*(AOMN)^2 + (AD1)*(C4S)^2 + (AOMN)^2*(AD2) + (C4S)^2*(AD2) + (AOMN)^2*(C3S) + (C4S)^2*(C3S) + (C4S)^2*(ABN1)*AXIB + (AOMN)*(ABT1B)*(ABN) + (ABN0)*(AOMN)^2*(AD5) + (ABN0)*(C4S)^2*(AD5) + (AOMN)^2*(ABN1)*AXIB);
        B2=(- (C1S)*(C2S) - (AOMN)^2*(AD2)*(C3S) + (AOMN)^2*(AD3)*(AZTB) - (C4S)^2*(AD2)*(C3S) + (C4S)^2*(AD3)*(AZTB) - (AD1)*(C1S)*(AOMN) - (AD1)*(AOMN)*(C2S) - (C1S)*(AOMN)*(AD2) + (AOMN)*(ABT2B)*(AD6) - (ABT1B)*(C2S)*(ABN) - (AOMN)*(C2S)*(AD2) - (C1S)*(AOMN)*(C3S) + (AOMN)*(AD4)*(GAMAB) - (AOMN)*(C2S)*(C3S) - (AD1)*(AOMN)^2*(AD2) - (AD1)*(C4S)^2*(AD2) - (AD1)*(AOMN)^2*(C3S) - (AD1)*(C4S)^2*(C3S) - (C1S)*(AOMN)*(ABN1)*(AXIB) - (AOMN)*(ABT1B)*(ABN)*(C3S) - (AOMN)*(ABT1B)*(ABN1)*(GAMAB) - (AOMN)*(AD4)*(ABN)*(AXIB) - (AOMN)*(C2S)*(ABN1)*(AXIB) - (ABN0)*(AOMN)^2*(AD5)*(C3S) - (ABN0)*(AOMN)^2*(AD3)*(AXIB) - (ABN0)*(C4S)^2*(AD5)*(C3S) - (ABN0)*(C4S)^2*(AD3)*(AXIB) - (AOMN)^2*(ABN1)*(AD5)*(AZTB) - (AOMN)^2*(ABN1)*(AD2)*(AXIB) - (C4S)^2*(ABN1)*(AD5)*(AZTB) - (C4S)^2*(ABN1)*(AD2)*(AXIB) - (ABN0)*(C1S)*(AOMN)*(AD5) - (ABN0)*(AOMN)*(ABT1B)*(AD6) - (ABN0)*(AOMN)*(C2S)*(AD5) - (AOMN)*(ABT1B)*(ABN)*(AD2) - (AOMN)*(ABT2B)*(ABN)*(AD5));
        B3= ((AD1)*(C1S)*(C2S) + (C1S)*(C2S)*(AD2) - (ABT2B)*(C2S)*(AD6) + (C1S)*(C2S)*(C3S) - (AD4)*(C2S)*(GAMAB) + (C1S)*(AOMN)*(AD2)*(C3S) + (C1S)*(C2S)*(ABN1)*(AXIB) - (AOMN)*(ABT2B)*(AD6)*(C3S) - (AOMN)*(AD4)*(AD6)*(AZTB) - (AOMN)*(AD4)*(AD2)*(GAMAB) + (ABT1B)*(C2S)*(ABN)*(C3S) + (ABT1B)*(C2S)*(ABN1)*(GAMAB) + (AD4)*(C2S)*(ABN)*(AXIB) - (C1S)*(AOMN)*(AD3)*(AZTB) - (AOMN)*(ABT2B)*(AD3)*(GAMAB) + (AOMN)*(C2S)*(AD2)*(C3S) - (AOMN)*(C2S)*(AD3)*(AZTB) + (AD1)*(AOMN)^2*(AD2)*(C3S) - (AD1)*(AOMN)^2*(AD3)*(AZTB) + (AD1)*(C4S)^2*(AD2)*(C3S) - (AD1)*(C4S)^2*(AD3)*(AZTB) + (AD1)*(C1S)*(AOMN)*(AD2) - (AD1)*(AOMN)*(ABT2B)*(AD6) + (ABN0)*(C1S)*(C2S)*(AD5) + (ABN0)*(ABT1B)*(C2S)*(AD6) + (AD1)*(AOMN)*(C2S)*(AD2) + (AD1)*(C1S)*(AOMN)*(C3S) - (AD1)*(AOMN)*(AD4)*(GAMAB) + (AD1)*(AOMN)*(C2S)*(C3S) + (ABT1B)*(C2S)*(ABN)*(AD2) + (ABT2B)*(C2S)*(ABN)*(AD5) + (ABN0)*(C1S)*(AOMN)*(AD5)*(C3S) + (ABN0)*(AOMN)*(ABT1B)*(AD6)*(C3S) + (ABN0)*(AOMN)*(AD4)*(AD6)*(AXIB) - (ABN0)*(AOMN)*(AD4)*(AD5)*(GAMAB) + (ABN0)*(C1S)*(AOMN)*(AD3)*(AXIB) + (ABN0)*(AOMN)*(ABT1B)*(AD3)*(GAMAB) + (ABN0)*(AOMN)*(C2S)*(AD5)*(C3S) + (ABN0)*(AOMN)*(C2S)*(AD3)*(AXIB) + (C1S)*(AOMN)*(ABN1)*(AD5)*(AZTB) + (C1S)*(AOMN)*(ABN1)*(AD2)*(AXIB) + (AOMN)*(ABT1B)*(ABN)*(AD2)*(C3S) + (AOMN)*(ABT1B)*(ABN1)*(AD6)*(AZTB) + (AOMN)*(ABT1B)*(ABN1)*(AD2)*(GAMAB) + (AOMN)*(ABT2B)*(ABN)*(AD5)*(C3S) - (AOMN)*(ABT2B)*(ABN1)*(AD6)*(AXIB) + (AOMN)*(ABT2B)*(ABN1)*(AD5)*(GAMAB) + (AOMN)*(AD4)*(ABN)*(AD5)*(AZTB) + (AOMN)*(AD4)*(ABN)*(AD2)*(AXIB) - (AOMN)*(ABT1B)*(ABN)*(AD3)*(AZTB) + (AOMN)*(ABT2B)*(ABN)*(AD3)*(AXIB) + (AOMN)*(C2S)*(ABN1)*(AD5)*(AZTB) + (AOMN)*(C2S)*(ABN1)*(AD2)*(AXIB));
        B4  = (- (C1S)*(C3S)*(AD2)*(C3S) + (ABT2B)*(C3S)*(AD6)*(C3S) + (AD4)*(C3S)*(AD6)*(AZTB) + (AD4)*(C3S)*(AD2)*(GAMAB) + (C1S)*(C3S)*(AD3)*(AZTB) + (ABT2B)*(C3S)*(AD3)*(GAMAB) - (AD1)*(C1S)*(C3S)*(AD2) + (AD1)*(ABT2B)*(C3S)*(AD6) - (AD1)*(C1S)*(C3S)*(C3S) + (AD1)*(AD4)*(C3S)*(GAMAB) - (AD1)*(C1S)*(AOMN)*(AD2)*(C3S) + (AD1)*(AOMN)*(ABT2B)*(AD6)*(C3S) + (AD1)*(AOMN)*(AD4)*(AD6)*(AZTB) + (AD1)*(AOMN)*(AD4)*(AD2)*(GAMAB) + (AD1)*(C1S)*(AOMN)*(AD3)*(AZTB) + (AD1)*(AOMN)*(ABT2B)*(AD3)*(GAMAB) - (ABN0)*(C1S)*(C3S)*(AD5)*(C3S) - (ABN0)*(ABT1B)*(C3S)*(AD6)*(C3S) - (ABN0)*(AD4)*(C3S)*(AD6)*(AXIB) + (ABN0)*(AD4)*(C3S)*(AD5)*(GAMAB) - (AD1)*(AOMN)*(C3S)*(AD2)*(C3S) - (ABN0)*(C1S)*(C3S)*(AD3)*(AXIB) - (ABN0)*(ABT1B)*(C3S)*(AD3)*(GAMAB) + (AD1)*(AOMN)*(C3S)*(AD3)*(AZTB) - (C1S)*(C3S)*(ABN1)*(AD5)*(AZTB) - (C1S)*(C3S)*(ABN1)*(AD2)*(AXIB) - (ABT1B)*(C3S)*(ABN)*(AD2)*(C3S) - (ABT1B)*(C3S)*(ABN1)*(AD6)*(AZTB) - (ABT1B)*(C3S)*(ABN1)*(AD2)*(GAMAB) - (ABT2B)*(C3S)*(ABN)*(AD5)*(C3S) + (ABT2B)*(C3S)*(ABN1)*(AD6)*(AXIB) - (ABT2B)*(C3S)*(ABN1)*(AD5)*(GAMAB) - (AD4)*(C3S)*(ABN)*(AD5)*(AZTB) - (AD4)*(C3S)*(ABN)*(AD2)*(AXIB) + (ABT1B)*(C3S)*(ABN)*(AD3)*(AZTB) - (ABT2B)*(C3S)*(ABN)*(AD3)*(AXIB));
        B5=((AD1)*(C1S)*(C3S)*(AD2)*(C3S) - (AD1)*(ABT2B)*(C3S)*(AD6)*(C3S) - (AD1)*(AD4)*(C3S)*(AD6)*(AZTB) - (AD1)*(AD4)*(C3S)*(AD2)*(GAMAB) - (AD1)*(C1S)*(C3S)*(AD3)*(AZTB) - (AD1)*(ABT2B)*(C3S)*(AD3)*(GAMAB));
        Y=[B0 B1 B2 B3 B4 B5];
        R=roots(Y);
        S=real(sqrt(R));
        S1=[S1 S];
    end

      plot(AD,S1(5,:));
  end
  legend

[Himanshu SINGLA]

For legend I use

if true
  legend('\Omega=0','\Omega=10','\Omega=50')
end

But it will show all the line of same colour(black) on print.

Isn't their option to draw dotted, dash dotted etc. line

[KSSV]

YOu can specify markers..read about plot.

[Himanshu SINGLA]

Sir, I tried this

if true
  plot(AD,S1(1,:),AD,S1(1,:),'+--',AD,S1(1,:),AD,S1(1,:),'*--')
end

but in figure, all the curves show the same marker (last marker style star).

I am plotting graph between AD and first root of 5 degree polynomial by taking different values of AOM=0,10,50,100.

I also tried following

if true
  plot(AD,S1(1,:),
%hold on
  plot(AD,S1(1,:),'+--')
% hold on
plot(AD,S1(1,:),'--')
% hold on
plot(AD,S1(1,:),'*--');
% hold on  end

but in this case also figure shows same marker style in all curves.