Nx=(1: nsteps) ';
X=(nx - 1) * dx;
For j=1: nsteps
A (j)=Mie (m, x (j));
End
Output_parameters='Qext, Qsca Qabs, Qb,
% plotting the results
M1=real (m); M2=imag (m);
Plot (x, a (:, 1:6))
Legend (' Qext ', 'Qsca', 'Qabs',' Qb ', '
The title (sprintf (' Mie Efficiencies, m=% g + % gi ', m1, m2))
Xlabel (' x ')
Result=a;