Setting of mass matrix

Hi, I have tried to set mass matrix

options = odeset('RelTol',1e-9,'AbsTol',1e-9*ones(608,1),'Mass', @(t)diag([ones(1,604) sin(t) 1 1 0]));

ode15s(@moven2,[0 L],vector,options,m,H0,x3,Ib,number,number1,number2,number3,mu);

but have obtained "Error using @(t)diag([ones(1,604),sin(t),1,1,0])

Too many input arguments"

Where is mistake?

Answer 1

[Torsten]

options = odeset('RelTol',1e-9,'AbsTol',1e-9*ones(608,1),'Mass', @(t,y)diag([ones(1,604) sin(t) 1 1 0]));

or

options = odeset('RelTol',1e-9,'AbsTol',1e-9*ones(608,1),'MStateDependence','none','Mass', @(t)diag([ones(1,604) sin(t) 1 1 0]));

[Oleksandr Maksimenko]

Thank you for the help. Unfortunately, nothing has changed

Error using @(t)diag([ones(1,604),sin(t),1,1,0])

Too many input arguments.

Error in odemass (line 59)

massM = feval(massFcn,t0,massArgs{:});

Error in ode15s (line 188)

[Mtype, Mt, Mfun, Margs, dMoptions] = odemass(FcnHandlesUsed,odeFcn,t0,y0,...

[Oleksandr Maksimenko]

I understood where I had a mistake. I rewrote the function of the right parts and this error disappeared.

options = odeset('RelTol',1e-9,'AbsTol',1e-9*ones(608,1),'Mass', @(t)diag([ones(1,604) sin(t) 1 1 0]));

ode15s(@(X,Y)moven2(X,Y,m,H0,x3,Ib,number,number1,number2,number3,mu),[0 L],vector,options);

However, a new error has occurred and I do not understand what it means and how to solve it.

"Error using daeic12 (line 77)

This DAE appears to be of index greater than 1."

[Torsten]

I'd check equation 605. It switches from algebraic at time t=0 to differential at t>0.