Correction
加速度 a=dv/dt =-bv/m;
dv/v=-b/m dt;
log(v)=c-bt/m;
v=exp(c-bt/m);
当t=0, v=sqrt(2gh/m); 所以exp(c)=sqrt(2gh/m)=v0;
v=v0*exp(-bt/m)==dy/dt;
dy=v0*exp(-bt/m)dt;
y=-mv0/b *exp(-bt/m)+y0;
当t=0, y=0, 所以y0=mv0/b;
y=mv0/b(1-exp(-bt/m));
y=mv0/b-mv/b=m(v0-v)/b;
v=v0+yb/m