l1=3.5*(1/10)^6                                      
l2=1.5*(1/10)^6            %电感

f=9                        %频率

u=4*pi*(1/10)^7            %真空磁导率
p=1.75*(1/10)^8            %铜的电导率
e=1/(36*pi)*((1/10)^(9))
b=(u*e)^(0.5)*f.*2*pi*(10^6)       %波印廷矢量
m=0.57*(1/10)^6                    %互感

r11=1280*(pi^6)*(0.6^2)*((f/300).^4)
r12=(pi*(f*(10^6))*u*p).^(0.5)*0.6*125
r21=1280*(pi^6)*(0.3^2)*((f/300).^4)
r22=(pi*(f*(10^6))*u*p).^(0.5)*0.3*125
r1=r11+r12
r2=r21+r22                        %电阻

u1=b*pi*0.6*0.6  
u2=b*pi*0.3*0.3                   %感应电压

w=f.*2*pi*10^6                   %角频率

z=90:1:110

c=1./(z.*w*i)                   %电容
a=(i*w*m-r2-i*w*l2)./(r1+i*w*l1+i*w*m) % i1=a*i2+b
b=(u1-u2)./(r1+i*w*l1+i*w*m)
g=(i*w*l1+r1+z*i).*b-u1         %i2表达式中的分子
d=(i*w*m-z*i)-a.*(r1+j*w*l1+z*i)%i2表达式中的分母

i2=g./d
i1=a.*i2+b

uo=abs((i1+i2).*z*i)

plot(z,uo)

上面是我的程序，感觉不对，求指点

Z是电容