Mathematica NDSolve different solutions same eq's

  • Context: Mathematica 
  • Thread starter Thread starter Vrbic
  • Start date Start date
  • Tags Tags
    Mathematica
Click For Summary
SUMMARY

The discussion centers on the discrepancies encountered when using Mathematica's NDSolve function to solve second-order differential equations derived from Lagrangian mechanics. The user reports that solving the equations directly yields different results compared to transforming the equations by introducing the constant of angular momentum, L. The key issue identified is the inconsistency in initial conditions, which leads to different solutions despite the equations appearing equivalent. This highlights the importance of ensuring that initial conditions are correctly set when using computational tools for solving differential equations.

PREREQUISITES
  • Understanding of Lagrangian mechanics and cyclic coordinates
  • Familiarity with Mathematica 12.3 and its NDSolve function
  • Knowledge of second-order differential equations
  • Concept of angular momentum in mechanical systems
NEXT STEPS
  • Explore the use of Mathematica's NDSolve for solving differential equations
  • Study the implications of cyclic coordinates in Lagrangian mechanics
  • Investigate how initial conditions affect solutions in differential equations
  • Learn about the transformation of equations in classical mechanics
USEFUL FOR

This discussion is beneficial for physicists, mathematicians, and engineers working with differential equations in Lagrangian mechanics, particularly those using Mathematica for computational analysis.

Vrbic
Messages
400
Reaction score
18
Hello, I have peculiar trouble with mathematica. I'm playing with lagrangians of some systems and I hit on "central problem". If I use my procedure in mathematica it solve 2 eq. of second order. But if you solve it by hand, you know, ##\phi## coord is "cyclic". And second eq. is reduced to ## \delta_t (mr^2 \phi)=0 ## and we introduce constant of angular momentum ##L##.
And the problem: If I let solve mathematica general equation without introduction ##L## I get another solution than if I transform second equation by introduction ##L##
(different eq.'s: ##2mr\phi'+mr^2\phi''=0## and ##mr^2\phi'=L##)
I used same initial condition.
Could you somebody advise where is the problem or how to resolve it?
 
Physics news on Phys.org
My bad, I found out that, the initial conditions aren't same. But it is for me still mystery. Why for this two same eq.'s in different forms, the initial conditions setted in same way aren't same?
 
As far as I can tell they are not the same equation.
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 5 ·
Replies
5
Views
4K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K