Wilson-θ法程序。威尔逊-θ法是在线性加速度法的一种拓展，其基本思路和实现方法是在时间段[t,t+θ△t]内线性变化，首先计算ti+θ△t时刻的运动，其中θ>1，然后通过内插得到时刻ti+△t的运动。可以证明当θ≥1.37时，威尔逊-θ法是无条件稳定的。