Wrong answer for ODE

2 vues (au cours des 30 derniers jours)
Umer
Umer le 29 Avr 2023
Commenté : Umer le 29 Avr 2023
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)
sol = 
MATLAB says:
sol =
(C4*exp(-(p*t)/m) - g*m)/p
What am I doing wrong?
  3 commentaires
Image Analyst
Image Analyst le 29 Avr 2023
Modifié(e) : Image Analyst le 29 Avr 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.
Umer
Umer le 29 Avr 2023
Oh yes!!! p IS a constant, that makes sense. sorry my bad. just learning this stuff for the first time. Thank youu so muchh.

Connectez-vous pour commenter.

Réponse acceptée

John D'Errico
John D'Errico le 29 Avr 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
ODE(t) = 
vsol = dsolve(ODE)
vsol = 
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?
  1 commentaire
Umer
Umer le 29 Avr 2023
Thank you so much!!

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Programming dans Help Center et File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by