%解范德波爾(Van der Pol)微分方程
tspan=[0,30];  %起始值0和終止值30
y0=[1;0];  %初始值
[t,y]=ode45(@vdpol,tspan,y0);  %解微分方程
y1=y(:,1);
y2=y(:,2);
figure(1)
plot(t,y1,':b',t,y2,'-r')  %畫微分方程解
figure(2)
plot(y1,y2)  %畫相平面圖