def equations(Y, t, m1, m2, c1, c2, c3, k1, k2, k3, F1, omega):
x1, v1, x2, v2 = Y
dx1dt = v1
dx1dt = v1
dv1dt = (-k1*x1 + k2*(x2 - x1) - c1*v1 - c2*(v1 - v2) + F1*np.sin(omega*t)) / m1
dx2dt = v2
dv2dt = (-k3*x2 - k2*(x2 - x1) - c3*v2 - c2*(v2 - v1)) / m2
return...