解四阶杆振动方程的辛算法

本文用中心差商代替高阶偏导数，将四阶杆振动方程转化成正则方程组，并利用辛欧拉中点格式数值求解。数值结果与理论分析相符。