[t,y] = ode45(@vdp1,[0 20],[2;0]);
plot(t,y(:,1), '-',t,y(:,2),':');
title('Solution of van der Pol Equation, \mu = 1')
xlabel('time t');
ylabel('solution y');
legend('y_1','y_2');
