介质谐振器天线谐振频率计算

本人对程序基本算是不懂，那位做过DRA，请教下：
上面的公式如何用matlab编程实现，求出f0，主要是想求出kx = m*pi/a, ky = n*pi/b;中m，n取不同值的时候，多个模式的谐振频率，谢谢了！