% ファンデルポール方程式を数値的に解く（固定ステップ幅）
interval = 0.00:0.01:10.00;
qinit = [2.00;0.00];
[time,q] = ode45(@van_der_Pol,interval,qinit);

% line style   solid -   broken -.   chain --   dotted :
plot(time,q(:,1),'-', time,q(:,2),'-.');
