function wzp

for k3=-20:20
    y0=ones(3,1);
    Y(:,k3+21)=fsolve(@(y)fun(y,k3),y0);
end
Y1=Y(1,:)
plot(-20:20,Y1)
end

function F=fun(y,k3)
k1=1;
k2=3;
r=1;
C=4;
Q=pi/3;
F=[(k2-k3)*i*y(1)-i*y(2)*y(3)-i*C*y(3)*exp(-i*Q);
    (-r-i*k2)*y(2)-i*C-i*y(3)*y(1);
    (-r-i*k1)*y(3)-i*8-i*y(2)*y(1)*exp(i*Q);];
end

如果我想k3以步长为0.1来画图，怎么控制？
请大虾指点！

k3=-20:0.1:20

没有你想的那样简单了
谢谢回复！