I've been trying for weeks to solve a system writen in its matrix form like:
M*y(t)'' + C*y(t)' + K*y(t) = F(t)
My system can have like 300 dof. If I generate all de equations and I use the dsolve/numeric command, it works but not all the time. I think that there is a problem with the Heaviside function that I have in some coefficients of my matrix.
I've been reading in forums that it's better to have a first order system to solve complex problems. For that we can create a procedure as explained in the dsolve/numeric/IVP help.
I would like to work with my matrix and vectors without generate all the equations. I think that that is possible but I'm not able to do it, even with a case of 2 dog with initial conditions set to 0.
I will be very tankfull if someone can help me with this problem.