[求助]求解一个化不开的函数的作图

[OKboy]

作者: watcher

文件名:f1.m
内容：
function F=f1(k,w)
F=tanh(pi/0.65*w*sqrt(k^2-1))-1/(17.36-0.715*i)*sqrt(k^2+17.36-0.715*i)/sqrt(k^2-1);
-----------------------------...

还有一个问题，如果我用双纵轴画图，怎么修改另外一条纵轴的scale和label呢？
clear all;
w=0:0.001:0.1;%这里设了步长0.001，可自行调节
for n=1:size(w,2)
k(n)=fsolve(@(k)f1(k,w(n)),6,optimset('Display','off','TolFun',1e-6));%用了fsolve，实际上是当作优化问题处理了
end
kr=real(k);
ki=imag(k);
plotyy(w,kr,w,ki);
axis([0,0.1,0,12,]);%另外一条纵轴ki我想设为scale为0：0.4，并加上标题