Plot 3rd differential equation
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Hello everyone, I could use a little help with plotting a 3rd differential equation usind "eqn" and "dsolve" syntax.
Here is my equation:
y''' - 3y' + 2y = x^2*e^x
0 commentaires
Réponses (2)
Star Strider
le 16 Jan 2022
(It will be necessary to determine if the left-hand-side are in terms of or .)
2 commentaires
Star Strider
le 16 Jan 2022
Thank you. That’s what I thought, since the Newton ‘prime’ notation implies derivatives of time.
It would also be nice to know if ‘x’ is a function of time (or anything else) or is simply a value. If it is a funciton of time, integrating the differential equation symbolically may not be possible (at least in a closed form, since the RHS would have to be integrated separately, then substituted back into the integrated LHS). That would slightly complicate things.
With what I currently know about this, it can just be plugged into the demonstration code I cited to in my original Answer to produce —
syms x y(t) t y0 D1y0 D2y0
sympref('AbbreviateOutput',false);
D1y = diff(y);
D2y = diff(D1y);
D3y = diff(D2y);
Eqn = D3y - 3*D1y + 2*y == x^2*exp(x)
ics = [y(0) == y0, D1y(0) == D1y0, D2y(0) == D2y0]; % Initial Conditions
Sy = dsolve(Eqn, ics);
Sy = simplify(Sy, 500)
I leave the rest to you.
.
Friedel Hartmann
le 16 Jan 2022
syms x y(x)
D1y = diff(y,1)
D3y = diff(y,3)
dsolve(D3y - 3*D1y + 2*y == x^2 * exp(x),x)
0 commentaires
Voir également
Catégories
En savoir plus sur Symbolic Math Toolbox dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!