to=0;
tf=1000;
x0=[0,32000];
[t,x]=ode23('vdpol', to, tf, x0, 0.00001,0);
