matlab解含空气阻力的动力学微分方程组时的奇怪现象

由于我想对下面的微分方程组进行数值求解：
u'＝-au²+bv
v'＝-av²-bv-c
（为上旋球的动力学方程，其中a为约化的空气阻力因子，u,v为x,y方向速度）
故现学matlab，试图给出xy平面内的轨迹图像。但是作出的u,v关于t的函数图像非常奇怪，基本上是一开始为零（不论初始条件如何），然后在t超过某一阈值后开始极速下降为负值。当我把空气阻力项删除后发现图像很正常；当我只保留空气阻力项时，也呈现优美的指数衰减图样。但是一旦在空气阻力项之外加上其他项就会发生上述奇怪现象。请问是什么原因？

刚性，可以用ode23s函数试一试