matlab解方程求分岔图的调用方程主程序

xiaoshihanlan

clear;
clc;
x0=[0.1;0;0.1;0];
global gamma
f=0.1:0.1:6;
for i=1:length(f)
disp(f(i));
gamma=f(i);
T=2*pi;
ts=0:T/100:200*T;
[t x]=ode45('507_2',ts,x0);
x0=x(end,:);
plot(f(i),x(400:100:end,1),'k','markersize',2);
xlabel('转速比(rad/s)');
ylabel('位移x');
hold on
end
请教前辈们，我这个程序有没有问题，调用的求解方程出来的图不对呢？急求啊