% STANDARD APPROACH 

clear all;
p=vteamparam(1);
[tplot,wplot]=ode45(@vteamRHS,[0 0.03],p.winit,odeset('MaxStep',1e-4),p);
clf
plot(tplot,1e9*wplot(:),'-b')

clear all;
p=vteamparam(2);
[tplot,wplot]=ode45(@vteamRHS,[0 0.03],p.winit,odeset('MaxStep',1e-4),p);
clf
plot(tplot,1e9*wplot(:),'-b')

% POINCARE MAP APPROACH
clear all;
[tplot,wplot]=vteamIntegrate(vteamparam(1),odeset('MaxStep',1e-4),0.03);
clf
plot(tplot,1e9*wplot(:),'-b')
