Maple code for a rectangular pulse
m: mass, r: damping constant, k: spring constant a: constant in the impulse function
eps: width of a rectangular pulse , n: number of frames