%无摄动力
clear;
clc;
%global omega
%f=0.2;
%for i=1:length(f)
%disp(f(i));
%omega=f(i);%变频率
%T=2*pi;
x0=[0;0];
tspan=[0,1000];
%[t,x]=ode45('duffing',tspan,x0);
%x0=x(end,:);%ba yi ge zhou qi de chu zhi geng xin
%tspan=0:0.001*2*pi:800;
options = odeset('RelTol', 1e-6,'AbsTol',1e-6);
[t,x]=ode45('duffing',tspan,x0);
plot(x(:,1),x(:,2));
hold on;
xlabel ('幅值');
ylabel ('位移x1');
%end