请教有关ode15s或ode23t解带质量矩阵的微分方程的问题！

sky8754180

如题，其中质量矩阵M(t,y)不是常数矩阵，而是时间t和状态变量y的函数且为奇异矩阵，形式大概为：
[m11(t,y),0,0,0;0,m21(t,y),0,0;0,0,m33(t,y),0;0,0,0,0]
有做过类似题目的大侠请指点一下，谢谢！