v_y(t)= \sqrt{\cfrac{g}{\alpha}} \tanh\left(-t\sqrt -\mbox{arctanh}\left(v_0\sqrt{\cfrac{\alpha}{g}}\right) \right) & \epsilon \le 0