function result = Mie_tetado(x,fact) % Computation of Mie Power Scattering and diffraction functions % and gi coefficients of Legendre Polynomial decomposition % for complex refractive-index ratio m=m'+im", size parameters x=k0*a, % according to Bohren and Huffman (1983) BEWI:TDD122 % 1) polar diagram, linear or in dB scale with respect to minimum, with % SL in upper semicircle, SR in lower semicircle and 3 cartesian diagrams % 2) same for SL0 and SR0 without diffraction pattern, % 3) scattered intensity S (lin or log scale), and degree of polarisation % 4) scattered intensity without diffraction peak S0 (lin or log scale), % 5) beam efficiencies of S and S0, diffraction efficiency Qd % 6) gi-factors (coefficients of Legendre Polynomials of Phase Function). % C. Mätzler, April 2004. m1=[1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.1,2.2,2.3,2.4,2.5]; nj=length(m1); Q=[];gi=[];g0i=[]; m=m1+0.001i; for j=1:nj, y=mie_teta(m(j),x,fact); y.Q=[m1(j),y.Q]; Q=[Q;y.Q]; gi=[gi;y.gi]; g0i=[g0i;y.g0i]; end; result.Q=Q; result.gi=gi; result.g0i=g0i;