用matlab画图Cp与f皆为变量,请附上程序指令代码

该给个解的区间的,害的我还推了半天.
clear;
clc;
C1=6330;
Cs=2720;
x=@(y)fsolve(@(x)4*Cs^3*sqrt((y^2-C1^2)*(y^2-Cs^2))*tan(pi*x*sqrt(y^2-Cs^2)/y/Cs)+C1*(2*Cs^2-y^2)^2*tan(pi*x*sqrt(y^2-C1^2)/y/C1),2000,optimset('display','off'));
yy=6350:0.5:8000;
for i=1:length(yy)
xx(i)=x(yy(i));
end
plot(xx,yy);
%区间取得是0.5,0.005的话电脑扛不住,感觉也不需要.Cp区间是6350到8000,你可以改,但必须Kp>6330,你自己看看吧