timestep=[0.00,10.00];
qinit=[2.00;0.00];
[time,q]=ode45(@van_der_Pol,timestep,qinit);

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