Re: matlab compatibility of od45 with adaptive stepping

[Carlo De Falco]

On 11 Oct 2015, at 09:36, Sebastian Schöps <address@hidden> wrote:

>> I don't have access to matlab right now, so could
>> someone please try the following and report the results?
>> 
>> fun = @(t, x) cos (t);
>> o = odeset ('AbsTol', 1e-11, 'RelTol', 1e-11);
>> [t, x] = ode45 (fun, [0, 2*pi], 0, o);
>> size (t), size (x)
>> [t, x] = ode45 (fun, [0, pi, 2*pi], 0, o);
>> size (t), size (x)
> 
> ans =
>   457     1
> ans =
>   457     1
> ans =
>     3     1
> ans =
>     3     1
> 
>> My interpratation of the documentation is that size (t) should be [3, 1] in 
>> the second case ...
> 
> It seems you are right :)

Unfortunately, this also means that integrate_adaptive needs an almost complete 
rewrite :(
c.