Time Integration scheme for non constant M, C and K matrices?

Can anyone suggest a time integration scheme for non constant mass (M), stiffness(K) and damping (C) matrices? I am trying to solve a dynamic system (Ma+Cv+Ku=R) where the matrices M,C and K are time dependent.

Any thoughts/ideas will be highly aprreciated. Thank you.