Edit: The answer below addressed the question of how to get results from nested functions, but doesn't adequately take into account the need to provide the correct function handle argument to ode15s. See the comments. Anyone more expert in using the ode solvers: please contribute!
One method is to return t and x from odestribeck in the normal way:
function Parent
...
[t,x,dx] = odestribeck(...);
plot(t,x);
...
function [t, x, dx] = odestribeck (t,x)
...
alg;
function alg
[t,x]=ode15s(ode,[0,5],[0,5]);
end
end
end
Note that the call to plot can go after the nested function definition if you like, but it's probably clearer to leave it with the other code in Parent, before the nested definitions - it doesn't affect the order of execution.
Alternatively, you can make t and x global to both the nested functions, by initialising them in Parent and not passing them as arguments:
function Parent
...
t = ...;
x = ...;
dx = odestribeck;
plot(t,x);
...
function dx = odestribeck
...
alg;
function alg
[t,x]=ode15s(ode,[0,5],[0,5]);
end
end
end
Again, you can put code after the nested function definitions if you wish.
