求高手帮忙改一下圆口径泰勒的程序

新手刚编程，希望大侠们指点一番
function I=Ip(rho,a)
%Ip为圆口径泰勒分布电流幅度
%a为圆口径半径
R0dB=30;
a=0.25;
R0=10^(R0dB/20);
A=log(R0+sqrt(R0^2-1))/pi;
n1=6;
r1=[1.2197 ;2.2331; 3.2383; 4.2411; 5.2428 ;6.2439;];%一阶bessel的前6个根/pi
phi=6.244/sqrt(A^2+(n1-0.5)^2);%波瓣展宽因子
for n=1:n1
    un=phi*sqrt(A^2+(n-0.5)^2);
end
rho=1:0.01:a;
g=1;
if m==0
    Bm=1;
else m=1:n1-1
       R11=1;
       R22=1;
    for n=1:n1-1
        R11=1-(r1(m,:)/un)^2;
        if n=m
           R22=1-(r1(m,:)/r1(1,n))^2;
        end
           R2=R22*R;
    end
        R1=R11*R;
       Bm(m)=-R1/R2/besselj(0,r1(1,m)*pi);
       g=g+Bm(m)*besselj(0,r1(1,m)*rho*pi/a);
end

已经解决了

能发那程序给我吗