接上段程序
elseif j>3*n-1 && j<=4*n-1 %12-15列
for i=1:4*(2*n-1)
if i<=n-1
h(i,j)=-(m-1)*besselk(m-1,v*r1)*sin(o2*(m-1))/r1;
elseif i>n-1 && i<=2*n-2
h(i,j)=w1^2*(m-1)*besselk(m-1,v*r2)*sin((m-1)*o2)/(v^2*r2);
elseif i>2*n-2 && i<=3*n-2
h(i,j)=-k^2*n2^2*v*((m-1)*besselk(m-1,v*r1)-besselk(m,v*r1))*...
cos((m-1)*o1)/(w*u0*beta);
elseif i>3*n-2 && i<=4*n-2
h(i,j)=w1^2*k^2*n2^2*v*((m-1)*besselk(m-1,v*r2)/(v*r2)-besselk(m,v*r2))*...
cos(o1*(m-1))/(u^2*w*u0*beta);
elseif i>4*n-2 && i<=5*n-2
h(i,j)=-besselk(m-1,v*r1)*cos((m-1)*o1);
elseif i>5*n-2 && i<=6*n-2
h(i,j)=besselk(m-1,v*r2)*cos(o1*(m-1));
elseif i>6*n-2 && i<=7*n-3
h(i,j)=0;