ASKSAGE: Sage Q&A Forum - RSS feedhttps://ask.sagemath.org/questions/Q&A Forum for SageenCopyright Sage, 2010. Some rights reserved under creative commons license.Wed, 21 Jun 2023 02:07:10 +0200Euler Lagrange Equations on Sagehttps://ask.sagemath.org/question/69312/euler-lagrange-equations-on-sage/I was wondering if it's possible to graph Euler Lagrange Equations on Sage. Especially if one of the equation has 3 parameters? There would be 2 equations. The 3 parameters below would be x, $\theta_1, \theta_2$.
I know how it's done in Mathematica, but wondering how to do it in Sage. Thanks! RedQueenWed, 21 Jun 2023 02:07:10 +0200https://ask.sagemath.org/question/69312/lagranian mechanicshttps://ask.sagemath.org/question/7856/lagranian-mechanics/I'm working on using sage to help with the Euler-Lagrange equation in my mechanics class. I have this worked up so far for a simple pendulum.
var('m,l,g,th,thdot,thdotdot,t')
th = function('th',t)
thdot = th.diff(t)
thdotdot = thdot.diff(t)
L = 1/2*m*l^2*thdot^2 - m*g*l*(1-cos(th))
diff(L,thdot)
but that gives off the error:
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "_sage_input_23.py", line 10, in <module>
exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("dmFyKCdtLGwsZyx0aCx0aGRvdCx0aGRvdGRvdCx0JykKdGggPSBmdW5jdGlvbigndGgnLHQpCnRoZG90ID0gdGguZGlmZih0KQp0aGRvdGRvdCA9IHRoZG90LmRpZmYodCkKTCA9IDEvMiptKmxeMip0aGRvdF4yIC0gbSpnKmwqKDEtY29zKHRoKSkKZGlmZihMLHRoZG90KQ=="),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))
File "", line 1, in <module>
File "/tmp/tmpTPoOKC/___code___.py", line 8, in <module>
exec compile(u'diff(L,thdot)
File "", line 1, in <module>
File "/home/eric/sage/local/lib/python2.6/site-packages/sage/calculus/functional.py", line 130, in derivative
return f.derivative(*args, **kwds)
File "expression.pyx", line 2502, in sage.symbolic.expression.Expression.derivative (sage/symbolic/expression.cpp:11917)
File "derivative.pyx", line 216, in sage.misc.derivative.multi_derivative (sage/misc/derivative.c:2191)
File "expression.pyx", line 2570, in sage.symbolic.expression.Expression._derivative (sage/symbolic/expression.cpp:12263)
TypeError: argument symb must be a symbol
I don't have any idea what all this means other than I'm guessing it doesn't like that I am trying to take the derivative with respect to a function?
Any help would be appreciated.ehremingtonSat, 08 Jan 2011 17:02:19 +0100https://ask.sagemath.org/question/7856/