DSolve[{Derivative[2][x][t] == c1 + a1*x[t] + b1*y[t], Derivative[2][y][t] == c2 + a2*x[t] + b2*y[t]},{x[t], y[t]},t]