{2,3}
ffoutsuper[flag_,LRc_,cL_,cR_,alphaL_,alphaR_,LRC_,CL_,CR_,AlphaL_,AlphaR_]:=(solz=eco2[LRC,CL,CR,AlphaL,AlphaR,10.,10.,10.,10.,tf];x1oldL=Max[0.0,(x1L[t]/.solz/.ttf)[[1]]];x2oldL=Max[0.0,(x2L[t]/.solz/.ttf)[[1]]];x3oldL=Max[0.0,(x3L[t]/.solz/.ttf)[[1]]];x4oldL=Max[0.0,(x4L[t]/.solz/.ttf)[[1]]];x1oldR=Max[0.0,(x1R[t]/.solz/.ttf)[[1]]];x2oldR=Max[0.0,(x2R[t]/.solz/.ttf)[[1]]];x3oldR=Max[0.0,(x3R[t]/.solz/.ttf)[[1]]];x4oldR=Max[0.0,(x4R[t]/.solz/.ttf)[[1]]];vecs=Re[Eigenvectors[matoutsuper[Max[0.0,LRc],cL,cR,x1oldL,x2oldL,x3oldL,x4oldL,x1oldR,x2oldR,x3oldR,x4oldR]]];vals=Re[Eigenvalues[matoutsuper[Max[0.0,LRc],cL,cR,x1oldL,x2oldL,x3oldL,x4oldL,x1oldR,x2oldR,x3oldR,x4oldR]]];For[i=1;j=0,i<=8,i++,sum1=vecs[[i,1]]^2+vecs[[i,2]]^2+vecs[[i,3]]^2+vecs[[i,4]]^2;sum2=vecs[[i,5]]^2+vecs[[i,6]]^2;sum3=vecs[[i,7]]^2+vecs[[i,8]]^2;sum4=sum2+sum3;Which[(sum1>=0.8)&&(Sign[vecs[[i,1]]]==Sign[vecs[[i,2]]])&&(Sign[vecs[[i,1]]]==Sign[vecs[[i,3]]])&&(Sign[vecs[[i,1]]]==Sign[vecs[[i,4]]]),i1=i,(sum2>0.8)&&(Sign[vecs[[i,5]]]==Sign[vecs[[i,6]]]),i2=i,(sum3>0.8)&&(Sign[vecs[[i,7]]]==Sign[vecs[[i,8]]]),i3=i,sum4>0.8&&(Sign[vecs[[i,5]]]==Sign[vecs[[i,6]]])&&(Sign[vecs[[i,7]]]==Sign[vecs[[i,8]]]),j++;list[j]=i,True,i1=i2=i3=Ordering[vals,-1][[1]]]];If[j>=2&&Abs[vals[list[j]]]<0.01&&Abs[vals[list[j-1]]]<0.01,i2=i3=list[j]];(*incaseofdegeneracy*)Which[flag1,val=vals[[i1]],flag==2,val=vals[[i2]],flag==3,val=vals[[i3]];(*vecs=Re[Eigenvectors[matoutsuper[Max[0.0,cR[[3]]],Max[0.0,cR[[4]]],alpha11,alpha21,alpha31,alpha41,alpha12,alpha22,alpha32,alpha42,alpha13,alpha23,alpha33,alpha43,alpha14,alpha24,alpha34,alpha44,cR[[3]],cR[[4]],Alpha11,Alpha21,Alpha31,Alpha41,Alpha12,Alpha22,Alpha32,Alpha42,Alpha13,Alpha23,Alpha33,Alpha43,Alpha14,Alpha24,Alpha34,Alpha44,x1oldL,x2oldL,x3oldL,x4oldL,x1oldR,x2oldR,x3oldR,x4oldR]]]]*)];val);