%Duffing間歇性混沌產生
%調用duffingpinlv函數
clear all
tspan=[0:0.004:600];%龍格-庫塔法求解方程，積分步長選0.04
x0=[1,1];%初始值
[t,z]=ode45('duffingpinlv',tspan,x0);
plot(t,z(:,1))%時序圖
%plot(z(:,1),z(:,2))%相圖