如果我想以步长为0.1画图,怎么控制?
function wzpfor 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来画图,怎么控制?
请大虾指点!