I just checked with 3 other places and their answers match with mine (by hand solution) but matlab can't be wrong, can it? i just wanna know if there's any coding error i made.
Wrong answer for ODE
2 views (last 30 days)
Show older comments
MATLAB is giving the wrong answer to my ordinary differential equation.
m*diff(v) + p*v(t) == -g*m
The solution I am getting by solving it by hand is
C*exp(-(p*t)/m) - (g*m)/p
This is part of my code:
syms v(t) p m g
sol = dsolve (m*diff(v) + p*v(t) == -g*m)
MATLAB says:
sol =
(C4*exp(-(p*t)/m) - g*m)/p
What am I doing wrong?
3 Comments
Image Analyst
on 29 Apr 2023
Edited: Image Analyst
on 29 Apr 2023
How did you get the one you got by hand? Why is your exp() not divided by p? Is p a constant? If so then it could be bundled into the constant.
Accepted Answer
John D'Errico
on 29 Apr 2023
As is so often the case, a DIFFERENT answer is not always a WRONG answer. It may just look different. But different is not always a bad thing. (Why do I think these are words we all should take more seriously? Skipping the editorial comments...)
We can look at the ODE. You give no initial condition, but that should not be relevant. As well, you do not say which release of MATLAB you did that with.
syms m p g v(t)
ODE = m*diff(v) + p*v(t) == -g*m
vsol = dsolve(ODE)
And that is the same as what you got from dsolve. But is it different? Remember that C1 in this solution is a COMPLETELY unknown parameter. As well, p is a CONSTANT, but it could be any constant.
So the only difference is that IF we decided to write
C = C1/p
then do you agree that the two solutions would now look the same?
More Answers (0)
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!