[snipercl]

以下是我的程序

b=5;
r=1000;
k=10;
T=0;
for n=0:1:80
v=n+0.5;
T=T+(-j)^(n+1)*(n+0.5)*besselj(v,k*b)*besselh(v,k*r)*lerangde(n,cos(theta))/besselh(v,k*b);
end


勒让德函数是我自己写的

function pn=lerangde(l,x)
pn=0;
for i=0:l/2
pn=pn+(-1)^i*jiecheng(2*l-2*i)*x.^(l-2*i)/(2^l*jiecheng(i)*jiecheng(l-i)*jiecheng(l-2*i));
end

阶乘函数
function m=jiecheng(n);
if (n==0)
m=1;
else
m=1;
for i=1:n;
m=m*i;
end
end