%% Evaluate van der Pol Equation
% This example solves the system y' = vdp1(t,y) using ode45, then plots the first component of the solution.
% Solve the system using ode45.

sol = ode45(@vdp1, [0 20], [2 0]);

%
% Evaluate the first component of the solution at 100 points in the interval [0 20].

x = linspace(0,20,100);
y = deval(sol,x,1);

%
% Plot the solution.

plot(x,y)

% Copyright 2012 The MathWorks, Inc.