请教糖葫芦低通滤波器的设计

大侠们，
我准备设计糖葫芦低通滤波器，但是由于阶跃电容的计算公式总是算的不正确，怀疑公式可能有问题，请你们帮我看看我的公式正确吗?
% spec cd=0.113e-12=1.13e-13
pi=3.14;
a=0.073/2;
b=0.897/2;
c=0.389/2;
rf=(b-c)/(b-a)
gm=b/a;
e0=8.854186e-12;
er=1;
e=e0*er;
cd1=[(e/(100*pi))*(((rf+1)/rf)*log((1+rf)/(1-rf))-2*log(4*rf/(1-rf)))+1.11e-15*(1-rf)*(gm-1)]*pi*b
cd2=(e/(100*pi))*[((rf*rf+1)/rf)*log((1+rf)/(1-rf))-2*log(4*rf/(1-rf*rf))]+1.11e-15*(1-rf)*(gm-1);
两个公式是两篇论文上的，但没一个准确计算的

你怎么知道不准呢?当然，公式不可能很准，实际做的时候要慢慢调

问题已经解决，公式有一点错误