function c1ex4
[t,x]=ode45('myvdpeq',[0,10],[-1;1]);  % 直接求微分方程數值解
% 下面的函數描述 Van de Pol 方程本身
function y=myvdpeq(t,x)
y=[x(2); -(x(1)^2-1)*x(2)-x(1)];

%延遲微分方程可以用 dde23() 函數求解，也可以用 Simulink 求解，后者更直觀
% 下面繪制出 Simulink 模型，選擇 Simulation/Start 菜單可以啟動求解程序
c1ex4mod