主要想观察a与s之间的关系，并画出二者之间的函数关系。
麻烦大家帮我看看，这个代码要怎么样编写呢？谢谢各位！

我的程序如下：但结果是错误的。
x1即为b，x2为c，k就是s，
function F=myfun(x,k)   
F=[x(1)-(((0.6/x(2))+x(1)*(k/(1-k)))/((1/x(2))+0.6*x(1)*(k/(1-k))))^0.5;x(2)-((0.6*k/((1-k)*x(1))+1)/(k/((1-k)*x(1))+0.6))];

x0=[1,1];%给x赋予初始值
options = optimset('Display','off');
s=0.1:0.05:0.95;
for i=1:1:length(s);
    kk=s(i);
    x=fsolve(@(x) myfun(x,kk), x0,options);
    x1(i)=x(1)
    x2(i)=x(2)
end
[m,n]=size(s)
s1=(s./(ones(m,n)-s)).^(((-0.28)/1.7)*ones(m,n))
a=x1.^(0.6*ones(m,n)).*s1.*x2.^(-0.6*ones(m,n))
plot(s,a)