|
fricas (1) -> jb := IJB('x,'u,'p,4,3);
Type: Type
fricas jbe := JBE jb;
Type: Type
fricas de := JDE(jb,jbe);
Type: Type
fricas ck := CKP(jb,jbe);
Type: Type
fricas eq1:jbe := P(1,[4]) + U(1)*P(1,[1]) + U(2)*P(1,[2]) + U(3)*P(1,[3])
Type: JetBundleExpression ?(IndexedJetBundle ?(x, u, p, 4, 3))
fricas eq2:jbe := P(2,[4]) + U(1)*P(2,[1]) + U(2)*P(2,[2]) + U(3)*P(2,[3])
Type: JetBundleExpression ?(IndexedJetBundle ?(x, u, p, 4, 3))
fricas eq3:jbe := P(3,[4]) + U(1)*P(3,[1]) + U(2)*P(3,[2]) + U(3)*P(3,[3])
Type: JetBundleExpression ?(IndexedJetBundle ?(x, u, p, 4, 3))
fricas eq4:jbe := P(1,[1]) + P(2,[2]) + P(3,[3])
Type: JetBundleExpression ?(IndexedJetBundle ?(x, u, p, 4, 3))
fricas euler:de := makeSystem([eq1,eq2,eq3,eq4])
Type: JetDifferentialEquation ?(IndexedJetBundle ?(x, u, p, 4, 3), JetBundleExpression ?(IndexedJetBundle ?(x, u, p, 4, 3)))
fricas setOutMode(14)$ck
Type: NonNegativeInteger ?
fricas setSimpMode(1)$ck
Type: NonNegativeInteger ?
Fix::
fricas complete(euler)$ck
Type: Void
|
![\label{eq5}\begin{array}{c}
{\ }
\
{{{p_{4}^{1}}+{{u^{3}}\ {p_{3}^{1}}}+{{u^{2}}\ {p_{2}^{1}}}+{{u^{1}}\ {p_{1}^{1}}}}= 0}
\
{\ }
\
{{{p_{4}^{2}}+{{u^{3}}\ {p_{3}^{2}}}+{{u^{2}}\ {p_{2}^{2}}}+{{u^{1}}\ {p_{1}^{2}}}}= 0}
\
{\ }
\
{{{p_{4}^{3}}+{{u^{3}}\ {p_{3}^{3}}}+{{u^{2}}\ {p_{2}^{3}}}+{{u^{1}}\ {p_{1}^{3}}}}= 0}
\
{\ }
\
{{{p_{3}^{3}}+{p_{2}^{2}}+{p_{1}^{1}}}= 0}
\](images/6647003692909907635-16.0px.png "\label{eq5}\begin{array}{c}
{\ }
\
{{{p_{4}^{1}}+{{u^{3}}\ {p_{3}^{1}}}+{{u^{2}}\ {p_{2}^{1}}}+{{u^{1}}\ {p_{1}^{1}}}}= 0}
\
{\ }
\
{{{p_{4}^{2}}+{{u^{3}}\ {p_{3}^{2}}}+{{u^{2}}\ {p_{2}^{2}}}+{{u^{1}}\ {p_{1}^{2}}}}= 0}
\
{\ }
\
{{{p_{4}^{3}}+{{u^{3}}\ {p_{3}^{3}}}+{{u^{2}}\ {p_{2}^{3}}}+{{u^{1}}\ {p_{1}^{3}}}}= 0}
\
{\ }
\
{{{p_{3}^{3}}+{p_{2}^{2}}+{p_{1}^{1}}}= 0}
\")
![\label{eq8}\](images/563200947207232556-16.0px.png "\label{eq8}\")
![\label{eq10}\](images/7480555134030225570-16.0px.png "\label{eq10}\")
![\label{eq12}\](images/4218519603404177416-16.0px.png "\label{eq12}\")
![\label{eq15}\begin{array}{c}
{\ }
\
{{{2 \ {p_{2}^{3}}\ {p_{3}^{2}}}+{2 \ {p_{1}^{3}}\ {p_{3}^{1}}}+{2 \ {{p_{2}^{2}}^{2}}}+{2 \ {p_{1}^{1}}\ {p_{2}^{2}}}+{2 \ {p_{1}^{2}}\ {p_{2}^{1}}}+{2 \ {{p_{1}^{1}}^{2}}}}= 0}
\](images/7926749830839756851-16.0px.png "\label{eq15}\begin{array}{c}
{\ }
\
{{{2 \ {p_{2}^{3}}\ {p_{3}^{2}}}+{2 \ {p_{1}^{3}}\ {p_{3}^{1}}}+{2 \ {{p_{2}^{2}}^{2}}}+{2 \ {p_{1}^{1}}\ {p_{2}^{2}}}+{2 \ {p_{1}^{2}}\ {p_{2}^{1}}}+{2 \ {{p_{1}^{1}}^{2}}}}= 0}
\")
![\label{eq16}\begin{array}{@{}l}
\displaystyle
= = = = = = = = = = = = = = = = = = =
\
\
\displaystyle
= = = = = = = = = =](images/3620217381178742251-16.0px.png "\label{eq16}\begin{array}{@{}l}
\displaystyle
= = = = = = = = = = = = = = = = = = =
\
\
\displaystyle
= = = = = = = = = =")
![\label{eq17}\](images/4534520014727666991-16.0px.png "\label{eq17}\")
![\label{eq19}\](images/7816524684722880981-16.0px.png "\label{eq19}\")
![\label{eq21}\](images/7607606247266326842-16.0px.png "\label{eq21}\")
![\label{eq24}\begin{array}{c}
{\ }
\
{{{p_{4}^{3}}+{{u^{3}}\ {p_{3}^{3}}}+{{u^{2}}\ {p_{2}^{3}}}+{{u^{1}}\ {p_{1}^{3}}}}= 0}
\
{\ }
\
{{{p_{4}^{2}}+{{u^{3}}\ {p_{3}^{2}}}+{{u^{2}}\ {p_{2}^{2}}}+{{u^{1}}\ {p_{1}^{2}}}}= 0}
\
{\ }
\
{{{p_{4}^{1}}+{{u^{3}}\ {p_{3}^{1}}}+{{u^{2}}\ {p_{2}^{1}}}+{{u^{1}}\ {p_{1}^{1}}}}= 0}
\
{\ }
\
{{{p_{3}^{3}}+{p_{2}^{2}}+{p_{1}^{1}}}= 0}
\
{\ }
\
{{{2 \ {p_{2}^{3}}\ {p_{3}^{2}}}+{2 \ {p_{1}^{3}}\ {p_{3}^{1}}}+{2 \ {{p_{2}^{2}}^{2}}}+{2 \ {p_{1}^{1}}\ {p_{2}^{2}}}+{2 \ {p_{1}^{2}}\ {p_{2}^{1}}}+{2 \ {{p_{1}^{1}}^{2}}}}= 0}
\](images/9094346296002367213-16.0px.png "\label{eq24}\begin{array}{c}
{\ }
\
{{{p_{4}^{3}}+{{u^{3}}\ {p_{3}^{3}}}+{{u^{2}}\ {p_{2}^{3}}}+{{u^{1}}\ {p_{1}^{3}}}}= 0}
\
{\ }
\
{{{p_{4}^{2}}+{{u^{3}}\ {p_{3}^{2}}}+{{u^{2}}\ {p_{2}^{2}}}+{{u^{1}}\ {p_{1}^{2}}}}= 0}
\
{\ }
\
{{{p_{4}^{1}}+{{u^{3}}\ {p_{3}^{1}}}+{{u^{2}}\ {p_{2}^{1}}}+{{u^{1}}\ {p_{1}^{1}}}}= 0}
\
{\ }
\
{{{p_{3}^{3}}+{p_{2}^{2}}+{p_{1}^{1}}}= 0}
\
{\ }
\
{{{2 \ {p_{2}^{3}}\ {p_{3}^{2}}}+{2 \ {p_{1}^{3}}\ {p_{3}^{1}}}+{2 \ {{p_{2}^{2}}^{2}}}+{2 \ {p_{1}^{1}}\ {p_{2}^{2}}}+{2 \ {p_{1}^{2}}\ {p_{2}^{1}}}+{2 \ {{p_{1}^{1}}^{2}}}}= 0}
\")
![\label{eq25}\](images/2614916660719994254-16.0px.png "\label{eq25}\")
