请求帮助！（知道lyapunov exponent的请点击）

cailiang

各位好：
       我是一垃圾学校研究生，明年毕业，现在准备毕业论文。导师带四个学生，有三个女生。确定论文方向时，导师把他的方向给她们三个了，叫我自己回去想想。我的妈呀，我要是想好了，我就成他导师了呢？时赶汶川大地震，我问他能否做点与地震有关的，我是学理论物理的。随即在电脑上给我打了一份文献，叫我回去看，然后写个调研报告。好的，上学期，我查阅相当多的文献，最终写了个报告。昨晚报告之后，他说，你看看你能做点什么吗？我的妈，这句话又来了。我说就做地震中的非线性、混沌之类的东西吧，他说那就试着做做吧。就这样，我开始在这地震模型的基础上，研究起了混混啊。于是把混沌的书看了些，于是就学者用fortran语言编程画相图，分岔图，截面图，微扰的轨道发散图什么的，这些都还差不多，可是当计算李雅普诺夫指数的时候，弄了我大半年的时间这么就是做不出来，去找导师，他说这东西我没算过，你找找文献，看看书。我但是无语。我的系统也不难啊，就是两个二阶的微分方程，在加上一个非线性的摩擦力但是系统是非自治的，两方程右端都含有时间t，眼瞅着就要毕业，心里急得慌啊。来到这里，我知道有很多会计算指数的高人，也在这里找到了几种计算的方法，我用的是格莱姆-斯密特正交化矩阵的方法来求指数的。如果有人能为我帮忙，在下感激不尽啊。真诚的希望大家能够帮助我！谢谢！谢谢！
      我的QQ:254875089(加我时注明指数);EMAIL:cailiang1984811@163.com
     附：运动方程

默默打酱油

我也最近在计算李雅普诺夫指数