How can I define and solve an ordinary differential equation with a random input (random inhomogeneous part)?

d2X(t)+c*dX(t)+k*X(t)=F(t)

where f(t) is a random function with normal distribution.

model1.mw

When I run the program it shows an error about one of the dependant variables.

In the attached file, first, I use dsolve to analytically solve a system of equations. Maple solves it easily.

But, when I change the frequency (in the Sine function) of the second equation, Maple cannot solve it. Why?

question.mw