非线性用fslove函数解，matlab求解非线性用的是迭代逼近的方法
function x = f1( k1,k2,k3,a,b,l)
f=@(x)([k1*x(1)-x(3)*(1-x(2))*b+l;k2*x(2)-x(3)*(1-x(1))*b;k3*x(3)-a*(1-x(1)*(1-x(2)))]) ;

x=fsolve(f,[1;1;1]);
end