Plot 3rd differential equation
2 次查看(过去 30 天)
显示 更早的评论
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 个评论
回答(2 个)
Star Strider
2022-1-16
(It will be necessary to determine if the left-hand-side are in terms of or .)
2 个评论
Star Strider
2022-1-16
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
2022-1-16
syms x y(x)
D1y = diff(y,1)
D3y = diff(y,3)
dsolve(D3y - 3*D1y + 2*y == x^2 * exp(x),x)
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Symbolic Math Toolbox 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!