% calculate w everywhere
% based on T1200470

LIM4=0.413;

q=-0.2325+13.3832i; % at MC3 HR
MC3IM2=[1,1.8445;0,1];
IM2=[1,0;-2/12.8,1];
IM2IM3=[1,1.1703;0,1];
IM3=[1,0;-2/(-6.24),1];
IM3IM4=[1,1.5889-LIM4;0,1];
IM4=[1,0;0,1];
IM4PRMAR=[1,LIM4;0,1];

%IM4PRMAR*IM4*IM3IM4;
M=IM3IM4*IM3*IM2IM3*IM2*MC3IM2;

v=[1;1/q];

%printMode(v);
printMode(M*v);
%a=printMode(M*v);
