% Mueller matrix manipulations
syms m11 m12 m13 m14 
syms m21 m22 m23 m24
syms m31 m32 m33 m34
syms m41 m42 m43 m44

syms s1 s2 s3 s4

M = [m11 m12 m13 m14;...
     m21 m22 m23 m24;...
     m31 m32 m33 m34;...
     m41 m42 m43 m44];...
     
S = [s1; s2; s3; s4];

Full = M*S;

s(1).s = [1;  1;  0;  0];   % LPH
s(2).s = [1; -1;  0;  0];   % LPV
s(3).s = [1;  0;  1;  0];   % LP+45
s(4).s = [1;  0; -1;  0];   % LP-45
s(5).s = [1;  0;  0;  1];   % RHC
s(6).s = [1;  0;  0; -1];   % LHC
 
for ns = 1:6
    so(ns).so = M*s(ns).s;
end
