Re: DAE daspk solver problem

[Sebastian Schöps]

Lazy_Tom wrote
> the output was :
> 
> warning: Option "RelTol" not set, new value 1.0e-006 is used
> warning: called from
>     radau5_test_1_V04 at line 19 column 6
> warning: Option "AbsTol" not set, new value 1.0e-006 is used
> warning: Option "NormControl" will be ignored by this solver
> warning: Option "InitialStep" not set, new value 1.0e-006 is used
> warning: Option "MaxStep" not set, new value 1.4e+000 is used
> warning: Option "Mass" only supports constant mass matrices M() and not
> M(t,y)

Did you see the last line? The solver does not like function handles for the
mass matrix. However, you can just use the matrix itself:

vmass =  diag([1, 1, 1, 1, 0]);
vopt = odeset ('Mass', vmass, 'NormControl', 'on');

I will try to look to your example in more detail next week.

Sebastian



--
View this message in context: 
http://octave.1599824.n4.nabble.com/DAE-daspk-solver-problem-tp4680377p4680419.html
Sent from the Octave - General mailing list archive at Nabble.com.