ab=[];
P=[0,100];
tab.push(P);
nbPoints=4000;
for (var j=1;j<nbPoints;j++) {
Q=[j*dt,P[1]+0.04*(20-P[1])*dt];
tab.push(Q);
P=Q;
}
tab