Solving constrained differential equations

I would like to solve coupled differential equations with some constraints. The detailed question appear below

[Torsten]

The people in this forum are not LATEX interpreters.

[Salil Kulkarni]

I am sorry. I could not figure out how to type matrices in the Mathworks forum. I have posted my question again

Answer 1

[Torsten]

Solving the above system is only possible in the form

M = [1 0 0 0 0; ...
     0 m1 0 0 0; ...
     0 0 1 0 0; ...
     0 0 0 m2 0; ...
     0 m2*tan(theta) 0 -m2*tan(theta) 0];
f = [y2 ; -sin(theta)*N ; y4 ; sin(theta)*N ; -m2*g+cos(theta)*N];
M*[y1dot ; y2dot ; y3dot ; y4dot ; Ndot] = f(y1,y2,y3,y4,N)

But you might get the message from the integrator that the index of the system is >1.

Test it.

[Salil Kulkarni]

Thank you very much. It works.