Multiple curves in a single graph

조회 수: 2(최근 30일)
Himanshu SINGLA
Himanshu SINGLA 2018년 5월 22일
답변: VBBV 2021년 10월 25일
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.

답변(2개)

KSSV
KSSV 2018년 5월 22일
USe :
plot(AD,S1);
instead of:
plot(AD,S1(5,:));
  댓글 수: 7
Himanshu SINGLA
Himanshu SINGLA 2018년 5월 30일
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.

댓글을 달려면 로그인하십시오.


VBBV
VBBV 2021년 10월 25일
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
plot(AD,S1(1,:),AD,S1(2,:),'+--',AD,S1(3,:),AD,S1(4,:),'*--', AD,S1(5,:),'s-')
legend
Which version of Matlab are you using ?

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by