RHO= (1-N)*RHO1 + N*RHO2;
BETA1= BETA1D*ET*T1S./OMG*RHO1;
BETA2= BETA1D*ET*T1S./OMG*RHO1;
KD= KB./(1-ET*OMG*T0*OMEGA)*(ET*OMG);
TAUM = T0S./(1-ET*OMG*T0*OMG)*(ET*OMG);
SD= SB./(1-ET*OMG*T0*OMEGA);
A1=(C11*KD*OMGS - C12D*KD*OMGS - A*BETA1*C12D*OMGS + B*BETA2*C11*OMGS - C11*C12D*OMGS*TAUM + C11*KD*OMGS*XI2B - C12D*KD*OMGS*XI1B);
A2=(KD*OMGS^2 + A*BETA1*OMGS^2 + B*BETA2*OMGS^2 + C11*OMGS^2*TAUM - C12D*OMGS^2*TAUM + KD*OMGS^2*XI1B + KD*OMGS^2*XI2B + A*BETA1*OMGS^2*XI2B - A*BETA2*OMGS^2*XI1B - B*BETA1*OMGS^2*XI2B + B*BETA2*OMGS^2*XI1B + C11*OMGS^2*TAUM*XI2B - C12D*OMGS^2*TAUM*XI1B);
A3= OMGS^3*TAUM + OMGS^3*TAUM*XI1B + OMGS^3*TAUM*XI2B;
r = solve(poly2sym(Y), 'MaxDegree', 3);
B1= C21*XI2B+C21-C22D*XI1B-C22D;
R = solve(poly2sym(Z), 'MaxDegree', 3);
ZETA1=(C11*K1S+(XI1B+1)*OMGS)*BETA2*OMGS+XI2B*BETA1*OMGS*OMGS./XI1B*BETA2*OMGS*OMGS+(-C12D*K1S+ (XI2B+1)*OMGS)*BETA1*OMGS;
ZETA2=(C11*K2S+(XI1B+1)*OMGS)*BETA2*OMGS+XI2B*BETA1*OMGS*OMGS./XI1B*BETA2*OMGS*OMGS+(-C12D*K2S+ (XI2B+1)*OMGS)*BETA1*OMGS;
ZETA3=(C11*K3S+(XI1B+1)*OMGS)*BETA2*OMGS+XI2B*BETA1*OMGS*OMGS./XI1B*BETA2*OMGS*OMGS+(-C12D*K3S+ (XI2B+1)*OMGS)*BETA1*OMGS;
ETA1=XI1B*OMGS*ZETA1-(C11*K1S+(XI1B+1)*OMGS);
ETA2=XI1B*OMGS*ZETA2-(C11*K2S+(XI1B+1)*OMGS);
ETA3=XI1B*OMGS*ZETA3-(C11*K3S+(XI1B+1)*OMGS);
XI4= C21*K4S+ (XI1B+1)*OMGS;
XI5= C21*K5S+ (XI1B+1)*OMGS;
A11= ALM1*K1S+2*MU1*M1S+BETA0*T1S*ETA1;
A12= ALM1*K2S+2*MU1*M2S+BETA0*T1S*ETA2;
A13= ALM1*K3S+2*MU1*M3S+BETA0*T1S*ETA3;
A41= (ALM2*K1S+2*MU2*M1S)*ZETA1;
A42= (ALM2*K2S+2*MU2*M2S)*ZETA2;
A43= (ALM2*K3S+2*MU2*M3S)*ZETA3;
M=[A11 A12 A13 A14 A15; A21 A22 A23 A24 A25; A31 A32 A33 A34 A35; A41 A42 A43 A44 A45; A51 A52 A53 A54 A55];
RRR = solve(D==0, K, 'MaxDegree', 4);
RRd
RRd =
-0.488289840687524 + 0.482618182235875i
-0.000000377515847 - 0.000001826506960i
-0.000001826150043 + 0.000000377239853i
-8.393573674339343 + 2.382886985058038i
0.488289840687524 - 0.482618182235875i
0.000000377515847 + 0.000001826506960i
0.000001826150043 - 0.000000377239853i
8.393573674339343 - 2.382886985058038i
●
Best Answer