% vdpolRrun.m, to run vdpolReqns.m by ralph abraham 14 oct 2003
% this creates two arrays by ode45 integration
% see help ode45 for the args (time span, inits, options, params)
options = odeset('OutputFcn',@odephas2)
[t,y]=ode45(@vdpolReqns,[0 20],[2 0],options,1.0);
% 3rd arg is OPTIONS placeholder
% and this plots the first column of y vs t
%%% plot(t,y(:,1));
% end