var tab=[];
var x0=xinit;
var v0=vinit;
tab.push([0,x0]);
nbPoints=1000;
for (var j=1;j<nbPoints;j++) {
x0+=v0*dt;
v0+=(-x0*omega*omega+epsilon*omega*(1-x0*x0)*v0-gamma*v0)*dt;
tab.push([j*dt,x0]);
};
tab