function [T,X] = vdp(X0,Tf,mu) ;
TSPAN = [0 Tf] ;
options = odeset('RelTol',1e-4,'AbsTol',1e-6) ;
[T,X] = ode45(@vdpODE,TSPAN,X0,options,mu) ;