tspan=[0 20]; % time span
yo = [2; 0]; % initial condition
ode23(@vdpol,tspan,yo); % solve the ODE-IVP
%[t,y] = ode23(@vdpol,tspan,yo); % solve the ODE-IVP
%plot(t,y(:,1),t,y(:,2),'--')