这个求功率耦合系数的，该怎么表达

问题遇到的现象和发生背景

img

img

img

问题相关代码，请勿粘贴截图
clear 
close all
V=5.0;
U=3.1527;
W=3.8808;
Npoint=501;
w0=17e-6;
lamble=1.55e-6;
a=25e-6;
k=2*pi/lamble;
x=linspace(-5,5,Npoint);
y=linspace(-5,5,Npoint);
deta=pi/4;
X=meshgrid(x,y);
Y=meshgrid(y,x);
r=sqrt(x.^2+y.^2);
EI=1/w0*exp(-r.^2/w0^2).*exp(-i*k*deta.*r);

if  (r<a)
  Em=besselj(1,U.*r/a)./besselj(1,U);
else
  Em=besselk(1,W.*r/a)./besselk(1,W);
end

up1=dblquad(@(x,y) EI.*Em,-5,5,-5,5);
up=abs(up1).^2;
down11=abs(EI).^2
down22=abs(Em).^2;
down1=dblquad(@(x,y) down11,-5,5,-5,5);
down2=dblquad(@(x,y) down12,-5,5,-5,5);
down=down1.*down2;
eta=real(up./down);
运行结果及报错内容
我的解答思路和尝试过的方法
我想要达到的结果