how to change time integration in pdepe

7 views (last 30 days)
ode15s is used in pdepe. How to change it to say ode23 etc when in need?
Torsten on 8 Feb 2024
I don't know how "pdepe" is arranged. If you can find the call to "ode15s" somewhere while editing the code, you could try to change it to a call to "ode23".
feynman feynman
feynman feynman on 8 Feb 2024
not really, because in pdepe there are functions like pdentrp that we can't use out of pdepe.

Sign in to comment.

Accepted Answer

Sanju on 28 Feb 2024
I understand that you want to change ODE solver whenever required,
‘pdepe’ is a function specifically designed for solving systems of partial differential equations in MATLAB. ‘pdepe’ internally selects an appropriate ODE solver based on the characteristics of the problem, and it does not allow you to specify ‘ode23’ or any other specific ODE solver.
If you need more control over the ODE solver used within ‘pdepe’, you can try to adjust the options provided to ‘pdepe’. For instance, you can specify tolerances, time steps, and other parameters that might affect the behaviour of the internal ODE solver.
Also, there’s an option for directly selecting the desired ODE solver for which you can directly use ‘ode23’ function call to solve the ODE,
Here’s an example code you can refer to,
% Define the ODE function
odefun = @(t, y) -y;
% Define the initial condition
y0 = 1;
% Define the time span
tspan = [0 10];
% Solve the ODE using ode23
[t, y] = ode23(odefun, tspan, y0);
Similarly, you can use other ODE solvers available in MATLAB such as ode45, ode23s, etc., by replacing ode23 with the desired solver function handle.
You can also refer to the following documentation links if required,
Hope this helps!
Torsten on 3 Mar 2024
Edited: Torsten on 3 Mar 2024
Because the other options are not related to the structure of the DAE problem, but to the performance of its solution.
If the solver cannot cope with the structure of the problem, this is a Ko criterion right at the beginning.

Sign in to comment.

More Answers (0)


Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!