求助高手看下程序怎么修改(对隐函数画图)

cypro

a0=2,deltae=0.01,ro0/a0=0.8;
c0=0.01,epson0=5.8e7,epsone=10*epson0;
mu0=pi*4e7,mue=mu0,mus=mu0;
deltas=1e5,eta0=120*pi;
for f=0:1e7:1e9                omiga=2*pi*f;
    gama=alfa+j*beta;
    gamae=j*mue*omiga*(deltae+j*epsone*omiga);
    gama0=sqrt(-omiga^2*mu0*epson0);
    zs=(sqrt(mus*omiga/(2*deltas)))*(1+j)/(2*pi*c0);
    v=sqrt(gama0^2-gama^2);
    u=sqrt(gamae^2-gama^2);
    eps=0;
    for m=-20:20            %这个小的for循环是对eps进行累加
        zm=-j*mue*omiga*diff('besselk(m,u*a0)',u)/(u*besselk(m,u*a0));
        ym=j*gamae^2*diff('besselk(m,u*a0)',u)/(u*besselk(m,u*a0));
        delta0=(j*m*gama/a0)^2*(1/v^2-1/u^2)^2/(((gama0/v)*diff('besseli(m,v*a0)',u)/besseli(m,v*a0))+(zm/eta0));
        rm=(((gama0/v)*diff('besselk(m,u*a0)',u)/besselk(m,u*a0))+ym+delta0)/(((gama0/v)*diff('besseli(m,v*a0)',u)/besseli(m,v*a0))+ym+delta0);
        eps=eps+rm*besselk(m,v*a0)/besseli(m,v*a0)*besseli(m,v*ro0)*besseli(m,v*(c0+ro0));
    end
    -j*mu0*omiga*v^2/(2*pi*gama0)*(besselk(0,v*c0)-eps)=zs;
    plot(f,alfa);
end
如上，倒数第三行是一个关于gama和omiga的方程，其他参数都由前述表示了。中间的小循环是对eps的求和。gama和omiga分别是f和alfa的函数，其他的函数关系如上书。现在要画f和alfa的函数关系图，感觉这样编程可能有点问题，本人才疏学浅，有同学能帮帮忙么。谢谢了。。。