function M=MieMueller(x,Qscat,TetaMie,S1,S2,Teta)
% normalized (to 4pi) Scattering Matrix for Mie scattering
% based on size parameter x, Qescat, S1,S2 Mie scatering ampliude and
% scattering angle Teta
S1new=interp1(TetaMie,S1,Teta);
S2new=interp1(TetaMie,S2,Teta);
S3=0;
S4=0;
M=4*pi/(pi*Qscat*x^2)*mueller(S1new,S2new,S3,S4);
