Acoustic Scattering from a Sphere - Steve Turley

% It always uses a minimum of 5 terms.x=-10:.1:10;y=x;k=2*pi; % wave number with wavelength of 1nmax=2*k*a+5;nary=0:nmax;nterm=2*nary+1;jn=sphj(nary,k*a);hn=sphh(nary,k*a);% legendre calculates the associated legendre polynomials.% I'm only interested in the m=0 onefor ix=1:size(x,2)for iy=1:size(y,2)r=sqrt(x(ix)^2+y(iy)^2);if(r