looking for completion of pendulum (constrained motion) equations using multiindex notation (1 independent variable [t] and 3 dependent variables [x,y,F]) fricas (1) -> jb:=IJB('t, Type: Type
fricas jbe:=JBE jb; Type: Type
fricas de := JDE(jb, Type: Type
fricas ck := CKP(jb, Type: Type
fricas eq1:jbe := m*P(4,
fricas eq2:jbe := m*P(5,
fricas eq3:jbe := U(1)^2 + U(2)^2 - L^2
fricas eq4:jbe := P(1,
fricas eq5:jbe := P(2,
fricas printSys([eq1,
Type: OutputForm?
fricas pendulum:de := makeSystem [eq1,
Type: JetDifferentialEquation?(IndexedJetBundle?(t,
fricas setOutMode(14)$ck
Type: NonNegativeInteger?
fricas setRedMode(1)$ck
Type: NonNegativeInteger?
fricas complete(pendulum)$ck
Type: Void
|
|
|
|
last edited 10 years ago by test1 |
![\label{eq6}\begin{array}{c}
{\ }
\
{{\frac{{L \ m \ {p_{1}^{4}}}+{{u^{1}}\ {u^{3}}}}{L}}= 0}
\
{\ }
\
{{\frac{{L \ m \ {p_{1}^{5}}}+{{u^{2}}\ {u^{3}}}+{L \ g \ m}}{L}}= 0}
\
{\ }
\
{{{{u^{2}}^{2}}+{{u^{1}}^{2}}-{{L}^{2}}}= 0}
\
{\ }
\
{{-{u^{4}}+{p_{1}^{1}}}= 0}
\
{\ }
\
{{-{u^{5}}+{p_{1}^{2}}}= 0}
\](images/8675463992067101343-16.0px.png "\label{eq6}\begin{array}{c}
{\ }
\
{{\frac{{L \ m \ {p_{1}^{4}}}+{{u^{1}}\ {u^{3}}}}{L}}= 0}
\
{\ }
\
{{\frac{{L \ m \ {p_{1}^{5}}}+{{u^{2}}\ {u^{3}}}+{L \ g \ m}}{L}}= 0}
\
{\ }
\
{{{{u^{2}}^{2}}+{{u^{1}}^{2}}-{{L}^{2}}}= 0}
\
{\ }
\
{{-{u^{4}}+{p_{1}^{1}}}= 0}
\
{\ }
\
{{-{u^{5}}+{p_{1}^{2}}}= 0}
\")
![\label{eq7}\begin{array}{c}
{\ }
\
{{{L \ m \ {p_{1}^{4}}}+{{u^{1}}\ {u^{3}}}}= 0}
\
{\ }
\
{{{L \ m \ {p_{1}^{5}}}+{{u^{2}}\ {u^{3}}}+{L \ g \ m}}= 0}
\
{\ }
\
{{-{u^{4}}+{p_{1}^{1}}}= 0}
\
{\ }
\
{{-{u^{5}}+{p_{1}^{2}}}= 0}
\
{\ }
\
{{{{u^{2}}^{2}}+{{u^{1}}^{2}}-{{L}^{2}}}= 0}
\](images/6844709053551123256-16.0px.png "\label{eq7}\begin{array}{c}
{\ }
\
{{{L \ m \ {p_{1}^{4}}}+{{u^{1}}\ {u^{3}}}}= 0}
\
{\ }
\
{{{L \ m \ {p_{1}^{5}}}+{{u^{2}}\ {u^{3}}}+{L \ g \ m}}= 0}
\
{\ }
\
{{-{u^{4}}+{p_{1}^{1}}}= 0}
\
{\ }
\
{{-{u^{5}}+{p_{1}^{2}}}= 0}
\
{\ }
\
{{{{u^{2}}^{2}}+{{u^{1}}^{2}}-{{L}^{2}}}= 0}
\")
![\label{eq10}\](images/7480555134030225570-16.0px.png "\label{eq10}\")
![\label{eq12}\](images/4218519603404177416-16.0px.png "\label{eq12}\")
![\label{eq15}\begin{array}{c}
{\ }
\
{{{L \ m \ {p_{1}^{5}}}+{{u^{2}}\ {u^{3}}}+{L \ g \ m}}= 0}
\
{\ }
\
{{{L \ m \ {p_{1}^{4}}}+{{u^{1}}\ {u^{3}}}}= 0}
\
{\ }
\
{{-{u^{5}}+{p_{1}^{2}}}= 0}
\
{\ }
\
{{-{u^{4}}+{p_{1}^{1}}}= 0}
\
{\ }
\
{{{{u^{2}}^{2}}+{{u^{1}}^{2}}-{{L}^{2}}}= 0}
\](images/4679836929869205428-16.0px.png "\label{eq15}\begin{array}{c}
{\ }
\
{{{L \ m \ {p_{1}^{5}}}+{{u^{2}}\ {u^{3}}}+{L \ g \ m}}= 0}
\
{\ }
\
{{{L \ m \ {p_{1}^{4}}}+{{u^{1}}\ {u^{3}}}}= 0}
\
{\ }
\
{{-{u^{5}}+{p_{1}^{2}}}= 0}
\
{\ }
\
{{-{u^{4}}+{p_{1}^{1}}}= 0}
\
{\ }
\
{{{{u^{2}}^{2}}+{{u^{1}}^{2}}-{{L}^{2}}}= 0}
\")
![\label{eq16}\](images/8817300687797269412-16.0px.png "\label{eq16}\")
