function result = Mie_xscan(m, nsteps, dx)
nx=(1:nsteps)';
x=(nx-1)*dx;
for j = 1:nsteps
a(j,:)=Mie(m,x(j));
end
output_parameters='Qext, Qsca, Qabs, Qb, <costeta>, Qb/Qsca';
% plotting the results
m1=real(m);m2=imag(m);
plot(x,a(:,1:6))
legend('Qext','Qsca','Qabs','Qb','<costeta>','Qb/Qsca')
title(sprintf('Mie Efficiencies, m=%g+%gi',m1,m2))
xlabel('x')
result=a;

轉載請註明出處,本文鏈接:https://www.uj5u.com/qita/29960.html
標籤:其他開發語言
