请教大家关于气动弹性的颤振问题

matlab007

各位大侠，小弟在用四阶龙格－库塔求解二元机翼的颤振速度时，结合机翼的结构方城，不会用时域方法模拟，关键是不知道怎样把这个方法运用到矩阵上面？用的是fortran语言，算出的结果，和我想要的误差比较大，不知是不是我得算法不对！！看哪位大侠懂，给小弟解释一下，好吗？谢谢了！
邮箱是：xiao_junzi@yahoo.com.cn
qq是：35361753