matlab计算为什么会出错？

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

用matlab计算参数，但是其中有些参数明显不对

问题相关代码，请勿粘贴截图

运行结果及报错内容

这两步计算结果是错的，大的离谱
a1=1iZotan(a/2);%阻抗值
a2=-1iZo1/(sin(a));
a1=0.00000000000000 + 1.63312393531954e+16i
a2=0.00000000000000 - 8.16561967659769e+15i

我的解答思路和尝试过的方法

我想要达到的结果

请问各位为什么会这样呢？我注意了弧度和角度的区别了，程序中用的就是弧度角，不知道为什么会和正确结果相差那么多？

V=340;
f=20000;
T=V/f;
An=T/4;%电极宽度
Bn=An;%电极间距
Ln=An+Bn;
K2=0.15;%K2机电耦合系数
Wn=50*10^(-6);%孔径大小
Wmax=Wn;
Sxx=3400;%介电常数张量
Szz=3130;
Sxz=0;
Zo=1;
fn=V/(2*Ln);
a=2*pi*f*Ln/V;%越度角θ
a1=1i*Zo*tan(a/2);%阻抗值
a2=-1i*Zo*1/(sin(a));
qn0=sin((pi/2)*(An/Ln));%qn
qn1=(1-qn0^2)^(1/2);%qn‘
qn01=ellipke(qn0);%椭圆积分
qn11=ellipke(qn1);
q=ellipke(2^(-1/2));
Cn=(Wn/2)*((Sxx*Szz-(Sxz)^2)^(1/2))*qn01/qn11;
rn=((2*fn*K2*Cn*Zo)^(1/2))*((Wn/Wmax)^(1/2))*(q/qn01);
%仿真元件参数计算
L=a1/(2i*pi*f);
C=1/(1i*a2*2*pi*f);

你这里的运算符都没有了

在这里面随便选一个语言的就可以，然后再贴代码

a=2*pi*f*Ln/V/360.0;%越度角θ

你重新把代码按照代码格式发出来吧，你现在的代码缺少了好多符号，看不懂了

角度转弧度，你乘以2π之后，忘记除以360了