> Then yes, you're absolutely right. I've never investigated complex problems
> before, so I expect there must be more changes that need to be done in the
> m-files to handle these kind of problems.
>> I think I will have some time to further investigate the codes (at least,
>> the m-files).
> Great, thanks, then I'll wait with that single modification to the m-file
> solvers until I either get a patch or modified m-files or whatever from you
> for all the other modifications.
So, after the change ' -> .' everywhere, I modified this block
--> vdelta = abs(y5 - y4);
--> vtau = max (vodeoptions.RelTol * abs(vu.'), vodeoptions.AbsTol);
vdelta = norm (y5 - y4, Inf);
vtau = max (vodeoptions.RelTol * max (norm (vu.', Inf), 1.0), ...
with the absolute values where pointed. I have to say that I don't
max (norm (vu.', Inf), 1.0)
Why the maximum between the norm and 1.0?
I modified in a similar way ode23, ode54 and
ode78 and tested the equation y'(t)=sqrt(i)*y(t) and a more complicated
Schroedinger equation and they work without any problem.
The change in vdelta makes the codes to behave like with "NormControl" =
"on", since the infinity-norm is used. I think the 2-norm should be used in