Clear Filters
Clear Filters

How do you set up an ode solver with a more than one function in the ode?

2 views (last 30 days)
For example, if you have an ode such as a*df(t)/dt + b*dg(t)/dt + f(t) + g(t) + c = 0 that you'd like to solve and if you want to set one of the derivatives to be zero, how do you do that?
(By "ode solver" in the title I mean solve numerically.)
Sam Chak
Sam Chak on 9 Oct 2023
I would like to suggest providing a specific example of this generic "Exact differential equation" form:
This would allow us to test which ODE solver might be suitable for solving it.

Sign in to comment.

Answers (2)

Torsten on 7 Oct 2023
Moved: Torsten on 7 Oct 2023
If you have two unknown functions (f and g), you need two equations.

Walter Roberson
Walter Roberson on 8 Oct 2023
syms f(t) g(t) a b c
df = diff(f);
dg = diff(g);
eqn = a*df + b*dg + f(t) + g(t) + c == 0
eqn(t) = 
sol = dsolve(eqn, g)
sol = struct with fields:
g: 0 f: exp(-t/a)*(C1 - c*exp(t/a))
eqn2 = subs(diff(sol.f), t, 0) == 0
eqn2 = 
constant_of_integration = setdiff( symvar(eqn2), [a b c])
constant_of_integration = 
solution_for_constant = solve(eqn2, constant_of_integration)
solution_for_constant = 
subs(subs(eqn, sol), constant_of_integration, solution_for_constant)
ans(t) = 
subs(sol, constant_of_integration, solution_for_constant)
ans = struct with fields:
g: 0 f: -c





Community Treasure Hunt

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

Start Hunting!