请各路大神帮忙：求不同m1、m2、f值时，四重积分求阻抗Zm的值,最后画出不同m1-m2下的Zm-f图

m1=[1,3,10];m2= [1,3,10]; lx=4;ly=3;      %m1、m2取同一值。

f=0:2000;  %频率f间隔可适当改变

rho0=1000;

k=2*pi*f/1500; %波数

omega=2*pi*f;

%积分范围：0<x<lx;0<y<ly;-20<kx<20;-20<ky<20;