(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.0' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 146314, 3262]*) (*NotebookOutlinePosition[ 146959, 3284]*) (* CellTagsIndexPosition[ 146915, 3280]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{\( (*Riemann\ invarients\ include\ any\ function\ of\ w\ and\ [ u + v*the\ second\ entry\ of\ LEV1*stuff\ + \ the\ integral\ at\ the\ bottom]*) \), "\[IndentingNewLine]", RowBox[{\(Clear[Ahat]\), "\[IndentingNewLine]", \(Ahat = {{a[w], b[w], \(a'\)[w]*u + \(b'\)[w]*v}, {c[w], d[w], \(c'\)[w]*u + \(d'\)[w]*v}, {0, 0, \(f'\)[w]}}\), "\[IndentingNewLine]", \(MatrixForm[Ahat]\), "\[IndentingNewLine]", \({EV1, EV2, EV3} = Eigenvectors[Ahat];\), "\[IndentingNewLine]", \({LEV1, LEV2, LEV3} = Inverse[Transpose[{EV1, EV2, EV3}]];\), "\[IndentingNewLine]", \(MatrixForm[{LEV1, LEV2, LEV3}]\), "\[IndentingNewLine]", \(Dop[ f_] := {\[PartialD]\_u\ f, \[PartialD]\_v\ f, \[PartialD]\_w\ f}\), "\[IndentingNewLine]", \(LB[V1_, V2_] := Dop[V1] . V2 - Dop[V2] . V1\), "\[IndentingNewLine]", "\[IndentingNewLine]", "\"\\"", "\[IndentingNewLine]", \(MatrixForm[Simplify[LB[EV1, EV2]]]\), "\[IndentingNewLine]", \(MatrixForm[Simplify[LB[EV2, EV3]]]\), "\[IndentingNewLine]", \(MatrixForm[Simplify[LB[EV1, EV3]]]\), "\[IndentingNewLine]", "\[IndentingNewLine]", "\"\\"", "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{ RowBox[{\(Simplify[LEV1\ ]\), "/", RowBox[{"(", RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{\((a[w] - d[w] + \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ \ d[w] + d[w]\^2\))\), " ", RowBox[{"(", RowBox[{ RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}], \(c[ w]\)]}], "+", RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{\(-d[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "+", RowBox[{\(b[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}]}], ")"}]}], ")"}]}], "/", RowBox[{"(", RowBox[{ "2", " ", \(\@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] + d[w]\^2\)\), " ", RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}], ")"}]}], "]"}]}]}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{\(a[w]\), ",", \(b[w]\), ",", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{\(c[w]\), ",", \(d[w]\), ",", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "}"}]}], "}"}]], "Output"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {\(a[w]\), \(b[w]\), RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}]}, {\(c[w]\), \(d[w]\), RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}]}, {"0", "0", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {\(-\(c[ w]\/\@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] + \ d[w]\^2\)\)\), \(-\(\(c[ w]\ \((\(-\(a[w]\/\(2\ c[w]\)\)\) + d[w]\/\(2\ c[w]\) - \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ \ a[w]\ d[w] + d[w]\^2\)\/\(2\ c[w]\))\)\)\/\@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ \ a[w]\ d[w] + d[w]\^2\)\)\), RowBox[{"-", FractionBox[ RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{\((\(-a[w]\) + d[w] - \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ \ d[w] + d[w]\^2\))\), " ", RowBox[{"(", RowBox[{ RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}], RowBox[{"2", " ", \(c[w]\), " ", RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "-", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(d[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}]], "+", FractionBox[ RowBox[{ RowBox[{\(-u\), " ", \(d[w]\), " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{"v", " ", \(d[w]\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"u", " ", \(b[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", \(b[w]\), " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], RowBox[{\(b[w]\ c[w]\), "-", \(a[w]\ d[w]\), "+", RowBox[{\(a[w]\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(d[w]\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", SuperscriptBox[ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "2"]}]]}], ")"}]}], \(\@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] + d[w]\^2\)\)]}]}, {\(c[ w]\/\@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] + d[w]\^2\ \)\), \(-\(\(c[w]\ \((a[w]\/\(2\ c[w]\) - d[w]\/\(2\ c[w]\) - \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ \ a[w]\ d[w] + d[w]\^2\)\/\(2\ c[w]\))\)\)\/\@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ \ a[w]\ d[w] + d[w]\^2\)\)\), RowBox[{"-", FractionBox[ RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{\((\(-a[w]\) + d[w] + \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ \ a[w]\ d[w] + d[w]\^2\))\), " ", RowBox[{"(", RowBox[{ RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}], RowBox[{"2", " ", \(c[w]\), " ", RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "-", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(d[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}]]}], "-", FractionBox[ RowBox[{ RowBox[{\(-u\), " ", \(d[w]\), " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{"v", " ", \(d[w]\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"u", " ", \(b[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", \(b[w]\), " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], RowBox[{\(b[w]\ c[w]\), "-", \(a[w]\ d[w]\), "+", RowBox[{\(a[w]\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(d[w]\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", SuperscriptBox[ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "2"]}]]}], ")"}]}], \(\@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] + d[w]\^2\)\)]}]}, {"0", "0", "1"} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"], Cell[BoxData[ \("The Lie Brackets are"\)], "Output"], Cell[BoxData[ InterpretationBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0"}, {"0"}, { FractionBox[ RowBox[{ RowBox[{\(c[w]\), " ", \((\(-a[w]\) + d[w])\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(b[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{\((\(-a[w]\) + d[w])\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"2", " ", \(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "-", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}]}], \(c[w]\^2\ \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ \ a[w]\ d[w] + d[w]\^2\)\)]} }], "\[NoBreak]", ")"}], MatrixForm[ {0, 0, Times[ Power[ c[ w], -2], Power[ Plus[ Power[ a[ w], 2], Times[ 4, b[ w], c[ w]], Times[ -2, a[ w], d[ w]], Power[ d[ w], 2]], Rational[ -1, 2]], Plus[ Times[ c[ w], Plus[ Times[ -1, a[ w]], d[ w]], Derivative[ 1][ b][ w]], Times[ b[ w], Plus[ Times[ Plus[ Times[ -1, a[ w]], d[ w]], Derivative[ 1][ c][ w]], Times[ 2, c[ w], Plus[ Derivative[ 1][ a][ w], Times[ -1, Derivative[ 1][ d][ w]]]]]]]]}]]], "Output"], Cell[BoxData[ InterpretationBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { FractionBox[ RowBox[{ RowBox[{"2", " ", \(c[w]\^2\), " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"2", " ", \(c[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{"(", RowBox[{\(-a[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}], "-", RowBox[{\((a[w] - d[w] + \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] + \ d[w]\^2\))\), " ", RowBox[{"(", RowBox[{ RowBox[{\(d[w]\), " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{\(b[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{ RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}]}], RowBox[{"2", " ", \(c[w]\), " ", RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}]]}, { FractionBox[ RowBox[{ RowBox[{"2", " ", \(c[w]\^2\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"2", " ", \(c[w]\), " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{"(", RowBox[{\(-a[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}], "-", RowBox[{\((a[w] - d[w] + \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] + \ d[w]\^2\))\), " ", RowBox[{"(", RowBox[{ RowBox[{\(d[w]\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{\(b[w]\), " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{ RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}]}], RowBox[{"2", " ", \(c[w]\), " ", RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}]]}, { FractionBox[ RowBox[{ FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{\((a[w] - d[w] + \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ \ d[w] + d[w]\^2\))\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "+", FractionBox[ RowBox[{ RowBox[{\(-2\), " ", RowBox[{"(", RowBox[{ RowBox[{\(c[w]\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(b[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "+", RowBox[{\(d[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "-", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}], "+", RowBox[{\(a[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], \(\@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] + d[w]\^2\)\)]}], ")"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{\(-d[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "+", RowBox[{\(b[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}], \(c[w]\^2\)], "+", RowBox[{"2", " ", RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "-", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{\(c[w]\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(b[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{\(d[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "-", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "-", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}], "+", FractionBox[ RowBox[{\((\(-a[w]\) + d[w] - \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] \ + d[w]\^2\))\), " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{\(-d[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "+", RowBox[{\(b[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{\(c[w]\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(b[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{\(a[w]\), " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{ RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{ RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(a[w]\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{"2", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(d[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{\(-u\), " ", \(d[w]\), " ", RowBox[{ SuperscriptBox["a", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"u", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["a", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{"v", " ", \(d[w]\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["b", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"u", " ", \(b[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", \(b[w]\), " ", RowBox[{ SuperscriptBox["d", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}], "+", RowBox[{ RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}]}], ")"}]}]}], ")"}]}], \(c[w]\)]}], RowBox[{"2", " ", SuperscriptBox[ RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}], "2"]}]]} }], "\[NoBreak]", ")"}], MatrixForm[ { Times[ Rational[ 1, 2], Power[ c[ w], -1], Power[ Plus[ Times[ b[ w], c[ w]], Times[ Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ -1, d[ w]], Derivative[ 1][ f][ w]]]], -1], Plus[ Times[ 2, Power[ c[ w], 2], Derivative[ 1][ a][ w]], Times[ 2, c[ w], Derivative[ 1][ c][ w], Plus[ Times[ -1, a[ w]], Derivative[ 1][ f][ w]]], Times[ -1, Plus[ a[ w], Times[ -1, d[ w]], Power[ Plus[ Power[ a[ w], 2], Times[ 4, b[ w], c[ w]], Times[ -2, a[ w], d[ w]], Power[ d[ w], 2]], Rational[ 1, 2]]], Plus[ Times[ d[ w], Derivative[ 1][ a][ w]], Times[ -1, b[ w], Derivative[ 1][ c][ w]], Times[ -1, Derivative[ 1][ a][ w], Derivative[ 1][ f][ w]]]]]], Times[ Rational[ 1, 2], Power[ c[ w], -1], Power[ Plus[ Times[ b[ w], c[ w]], Times[ Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ -1, d[ w]], Derivative[ 1][ f][ w]]]], -1], Plus[ Times[ 2, Power[ c[ w], 2], Derivative[ 1][ b][ w]], Times[ 2, c[ w], Derivative[ 1][ d][ w], Plus[ Times[ -1, a[ w]], Derivative[ 1][ f][ w]]], Times[ -1, Plus[ a[ w], Times[ -1, d[ w]], Power[ Plus[ Power[ a[ w], 2], Times[ 4, b[ w], c[ w]], Times[ -2, a[ w], d[ w]], Power[ d[ w], 2]], Rational[ 1, 2]]], Plus[ Times[ d[ w], Derivative[ 1][ b][ w]], Times[ -1, b[ w], Derivative[ 1][ d][ w]], Times[ -1, Derivative[ 1][ b][ w], Derivative[ 1][ f][ w]]]]]], Times[ Rational[ 1, 2], Power[ Plus[ Times[ b[ w], c[ w]], Times[ Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ -1, d[ w]], Derivative[ 1][ f][ w]]]], -2], Plus[ Times[ Power[ c[ w], -2], Plus[ Times[ Plus[ a[ w], Times[ -1, d[ w]], Power[ Plus[ Power[ a[ w], 2], Times[ 4, b[ w], c[ w]], Times[ -2, a[ w], d[ w]], Power[ d[ w], 2]], Rational[ 1, 2]]], Derivative[ 1][ c][ w]], Times[ c[ w], Plus[ Times[ -1, Derivative[ 1][ a][ w]], Derivative[ 1][ d][ w], Times[ Power[ Plus[ Power[ a[ w], 2], Times[ 4, b[ w], c[ w]], Times[ -2, a[ w], d[ w]], Power[ d[ w], 2]], Rational[ -1, 2]], Plus[ Times[ -2, Plus[ Times[ c[ w], Derivative[ 1][ b][ w]], Times[ b[ w], Derivative[ 1][ c][ w]]]], Times[ d[ w], Plus[ Derivative[ 1][ a][ w], Times[ -1, Derivative[ 1][ d][ w]]]], Times[ a[ w], Plus[ Times[ -1, Derivative[ 1][ a][ w]], Derivative[ 1][ d][ w]]]]]]]], Plus[ Times[ -1, d[ w], Plus[ Times[ u, Derivative[ 1][ a][ w]], Times[ v, Derivative[ 1][ b][ w]]]], Times[ b[ w], Plus[ Times[ u, Derivative[ 1][ c][ w]], Times[ v, Derivative[ 1][ d][ w]]]], Times[ Plus[ Times[ u, Derivative[ 1][ a][ w]], Times[ v, Derivative[ 1][ b][ w]]], Derivative[ 1][ f][ w]]], Plus[ Times[ b[ w], c[ w]], Times[ Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ -1, d[ w]], Derivative[ 1][ f][ w]]]]], Times[ 2, Plus[ Times[ b[ w], c[ w]], Times[ Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ -1, d[ w]], Derivative[ 1][ f][ w]]]], Plus[ Times[ Plus[ Times[ u, Derivative[ 1][ a][ w]], Times[ v, Derivative[ 1][ b][ w]]], Derivative[ 1][ c][ w]], Times[ c[ w], Plus[ Times[ u, Derivative[ 2][ a][ w]], Times[ v, Derivative[ 2][ b][ w]]]], Times[ -1, Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ u, Derivative[ 2][ c][ w]], Times[ v, Derivative[ 2][ d][ w]]]], Times[ -1, Plus[ Times[ u, Derivative[ 1][ c][ w]], Times[ v, Derivative[ 1][ d][ w]]], Plus[ Derivative[ 1][ a][ w], Times[ -1, Derivative[ 2][ f][ w]]]]]], Times[ -2, Plus[ Times[ c[ w], Plus[ Times[ u, Derivative[ 1][ a][ w]], Times[ v, Derivative[ 1][ b][ w]]]], Times[ -1, Plus[ Times[ u, Derivative[ 1][ c][ w]], Times[ v, Derivative[ 1][ d][ w]]], Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]]]], Plus[ Times[ c[ w], Derivative[ 1][ b][ w]], Times[ b[ w], Derivative[ 1][ c][ w]], Times[ -1, Plus[ d[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Derivative[ 1][ a][ w], Times[ -1, Derivative[ 2][ f][ w]]]], Times[ -1, Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Derivative[ 1][ d][ w], Times[ -1, Derivative[ 2][ f][ w]]]]]], Times[ Power[ c[ w], -1], Plus[ Times[ -1, a[ w]], d[ w], Times[ -1, Power[ Plus[ Power[ a[ w], 2], Times[ 4, b[ w], c[ w]], Times[ -2, a[ w], d[ w]], Power[ d[ w], 2]], Rational[ 1, 2]]]], Plus[ Times[ Plus[ Times[ -1, d[ w], Plus[ Times[ u, Derivative[ 1][ a][ w]], Times[ v, Derivative[ 1][ b][ w]]]], Times[ b[ w], Plus[ Times[ u, Derivative[ 1][ c][ w]], Times[ v, Derivative[ 1][ d][ w]]]], Times[ Plus[ Times[ u, Derivative[ 1][ a][ w]], Times[ v, Derivative[ 1][ b][ w]]], Derivative[ 1][ f][ w]]], Plus[ Times[ c[ w], Derivative[ 1][ b][ w]], Times[ b[ w], Derivative[ 1][ c][ w]], Times[ -1, a[ w], Derivative[ 1][ d][ w]], Times[ Derivative[ 1][ a][ w], Derivative[ 1][ f][ w]], Times[ Derivative[ 1][ d][ w], Derivative[ 1][ f][ w]], Times[ a[ w], Derivative[ 2][ f][ w]], Times[ -2, Derivative[ 1][ f][ w], Derivative[ 2][ f][ w]], Times[ d[ w], Plus[ Times[ -1, Derivative[ 1][ a][ w]], Derivative[ 2][ f][ w]]]]], Times[ -1, Plus[ Times[ b[ w], c[ w]], Times[ Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ -1, d[ w]], Derivative[ 1][ f][ w]]]], Plus[ Times[ -1, u, d[ w], Derivative[ 2][ a][ w]], Times[ u, Derivative[ 1][ f][ w], Derivative[ 2][ a][ w]], Times[ -1, v, d[ w], Derivative[ 2][ b][ w]], Times[ v, Derivative[ 1][ f][ w], Derivative[ 2][ b][ w]], Times[ u, b[ w], Derivative[ 2][ c][ w]], Times[ v, b[ w], Derivative[ 2][ d][ w]], Times[ u, Derivative[ 1][ a][ w], Plus[ Times[ -1, Derivative[ 1][ d][ w]], Derivative[ 2][ f][ w]]], Times[ Derivative[ 1][ b][ w], Plus[ Times[ u, Derivative[ 1][ c][ w]], Times[ v, Derivative[ 2][ f][ w]]]]]]]]]]}]]], "Output"], Cell[BoxData[ InterpretationBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { FractionBox[ RowBox[{ RowBox[{"2", " ", \(c[w]\^2\), " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"2", " ", \(c[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{"(", RowBox[{\(-a[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}], "+", RowBox[{\((\(-a[w]\) + d[w] + \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] + \ d[w]\^2\))\), " ", RowBox[{"(", RowBox[{ RowBox[{\(d[w]\), " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{\(b[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{ RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}]}], RowBox[{"2", " ", \(c[w]\), " ", RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}]]}, { FractionBox[ RowBox[{ RowBox[{"2", " ", \(c[w]\^2\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"2", " ", \(c[w]\), " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{"(", RowBox[{\(-a[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}], "+", RowBox[{\((\(-a[w]\) + d[w] + \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] + \ d[w]\^2\))\), " ", RowBox[{"(", RowBox[{ RowBox[{\(d[w]\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{\(b[w]\), " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{ RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}]}], RowBox[{"2", " ", \(c[w]\), " ", RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}]]}, { FractionBox[ RowBox[{ RowBox[{"-", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{\((\(-a[w]\) + d[w] + \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ \ a[w]\ d[w] + d[w]\^2\))\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "+", FractionBox[ RowBox[{ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{\(c[w]\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(b[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "+", RowBox[{\(a[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "-", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}], "+", RowBox[{\(d[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], \(\@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] + d[w]\^2\)\)]}], ")"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{\(-d[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "+", RowBox[{\(b[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}], \(c[w]\^2\)]}], "+", RowBox[{"2", " ", RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "-", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{\(c[w]\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(b[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{\(d[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "-", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], "-", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}], "+", FractionBox[ RowBox[{\((\(-a[w]\) + d[w] + \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] \ + d[w]\^2\))\), " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{\(-d[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "+", RowBox[{\(b[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{\(c[w]\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(b[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{\(a[w]\), " ", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{ RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{ RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(a[w]\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{"2", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{\(d[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{\(-u\), " ", \(d[w]\), " ", RowBox[{ SuperscriptBox["a", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"u", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["a", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "-", RowBox[{"v", " ", \(d[w]\), " ", RowBox[{ SuperscriptBox["b", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{ SuperscriptBox["b", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"u", " ", \(b[w]\), " ", RowBox[{ SuperscriptBox["c", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", \(b[w]\), " ", RowBox[{ SuperscriptBox["d", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"u", " ", RowBox[{ SuperscriptBox["a", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox["d", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}], "+", RowBox[{ RowBox[{ SuperscriptBox["b", "\[Prime]", MultilineFunction->None], "[", "w", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], "+", RowBox[{"v", " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "w", "]"}]}]}], ")"}]}]}], ")"}]}]}], ")"}]}], \(c[w]\)]}], RowBox[{"2", " ", SuperscriptBox[ RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}], "2"]}]]} }], "\[NoBreak]", ")"}], MatrixForm[ { Times[ Rational[ 1, 2], Power[ c[ w], -1], Power[ Plus[ Times[ b[ w], c[ w]], Times[ Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ -1, d[ w]], Derivative[ 1][ f][ w]]]], -1], Plus[ Times[ 2, Power[ c[ w], 2], Derivative[ 1][ a][ w]], Times[ 2, c[ w], Derivative[ 1][ c][ w], Plus[ Times[ -1, a[ w]], Derivative[ 1][ f][ w]]], Times[ Plus[ Times[ -1, a[ w]], d[ w], Power[ Plus[ Power[ a[ w], 2], Times[ 4, b[ w], c[ w]], Times[ -2, a[ w], d[ w]], Power[ d[ w], 2]], Rational[ 1, 2]]], Plus[ Times[ d[ w], Derivative[ 1][ a][ w]], Times[ -1, b[ w], Derivative[ 1][ c][ w]], Times[ -1, Derivative[ 1][ a][ w], Derivative[ 1][ f][ w]]]]]], Times[ Rational[ 1, 2], Power[ c[ w], -1], Power[ Plus[ Times[ b[ w], c[ w]], Times[ Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ -1, d[ w]], Derivative[ 1][ f][ w]]]], -1], Plus[ Times[ 2, Power[ c[ w], 2], Derivative[ 1][ b][ w]], Times[ 2, c[ w], Derivative[ 1][ d][ w], Plus[ Times[ -1, a[ w]], Derivative[ 1][ f][ w]]], Times[ Plus[ Times[ -1, a[ w]], d[ w], Power[ Plus[ Power[ a[ w], 2], Times[ 4, b[ w], c[ w]], Times[ -2, a[ w], d[ w]], Power[ d[ w], 2]], Rational[ 1, 2]]], Plus[ Times[ d[ w], Derivative[ 1][ b][ w]], Times[ -1, b[ w], Derivative[ 1][ d][ w]], Times[ -1, Derivative[ 1][ b][ w], Derivative[ 1][ f][ w]]]]]], Times[ Rational[ 1, 2], Power[ Plus[ Times[ b[ w], c[ w]], Times[ Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ -1, d[ w]], Derivative[ 1][ f][ w]]]], -2], Plus[ Times[ -1, Power[ c[ w], -2], Plus[ Times[ Plus[ Times[ -1, a[ w]], d[ w], Power[ Plus[ Power[ a[ w], 2], Times[ 4, b[ w], c[ w]], Times[ -2, a[ w], d[ w]], Power[ d[ w], 2]], Rational[ 1, 2]]], Derivative[ 1][ c][ w]], Times[ -1, c[ w], Plus[ Times[ -1, Derivative[ 1][ a][ w]], Derivative[ 1][ d][ w], Times[ Power[ Plus[ Power[ a[ w], 2], Times[ 4, b[ w], c[ w]], Times[ -2, a[ w], d[ w]], Power[ d[ w], 2]], Rational[ -1, 2]], Plus[ Times[ 2, Plus[ Times[ c[ w], Derivative[ 1][ b][ w]], Times[ b[ w], Derivative[ 1][ c][ w]]]], Times[ a[ w], Plus[ Derivative[ 1][ a][ w], Times[ -1, Derivative[ 1][ d][ w]]]], Times[ d[ w], Plus[ Times[ -1, Derivative[ 1][ a][ w]], Derivative[ 1][ d][ w]]]]]]]], Plus[ Times[ -1, d[ w], Plus[ Times[ u, Derivative[ 1][ a][ w]], Times[ v, Derivative[ 1][ b][ w]]]], Times[ b[ w], Plus[ Times[ u, Derivative[ 1][ c][ w]], Times[ v, Derivative[ 1][ d][ w]]]], Times[ Plus[ Times[ u, Derivative[ 1][ a][ w]], Times[ v, Derivative[ 1][ b][ w]]], Derivative[ 1][ f][ w]]], Plus[ Times[ b[ w], c[ w]], Times[ Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ -1, d[ w]], Derivative[ 1][ f][ w]]]]], Times[ 2, Plus[ Times[ b[ w], c[ w]], Times[ Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ -1, d[ w]], Derivative[ 1][ f][ w]]]], Plus[ Times[ Plus[ Times[ u, Derivative[ 1][ a][ w]], Times[ v, Derivative[ 1][ b][ w]]], Derivative[ 1][ c][ w]], Times[ c[ w], Plus[ Times[ u, Derivative[ 2][ a][ w]], Times[ v, Derivative[ 2][ b][ w]]]], Times[ -1, Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ u, Derivative[ 2][ c][ w]], Times[ v, Derivative[ 2][ d][ w]]]], Times[ -1, Plus[ Times[ u, Derivative[ 1][ c][ w]], Times[ v, Derivative[ 1][ d][ w]]], Plus[ Derivative[ 1][ a][ w], Times[ -1, Derivative[ 2][ f][ w]]]]]], Times[ -2, Plus[ Times[ c[ w], Plus[ Times[ u, Derivative[ 1][ a][ w]], Times[ v, Derivative[ 1][ b][ w]]]], Times[ -1, Plus[ Times[ u, Derivative[ 1][ c][ w]], Times[ v, Derivative[ 1][ d][ w]]], Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]]]], Plus[ Times[ c[ w], Derivative[ 1][ b][ w]], Times[ b[ w], Derivative[ 1][ c][ w]], Times[ -1, Plus[ d[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Derivative[ 1][ a][ w], Times[ -1, Derivative[ 2][ f][ w]]]], Times[ -1, Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Derivative[ 1][ d][ w], Times[ -1, Derivative[ 2][ f][ w]]]]]], Times[ Power[ c[ w], -1], Plus[ Times[ -1, a[ w]], d[ w], Power[ Plus[ Power[ a[ w], 2], Times[ 4, b[ w], c[ w]], Times[ -2, a[ w], d[ w]], Power[ d[ w], 2]], Rational[ 1, 2]]], Plus[ Times[ Plus[ Times[ -1, d[ w], Plus[ Times[ u, Derivative[ 1][ a][ w]], Times[ v, Derivative[ 1][ b][ w]]]], Times[ b[ w], Plus[ Times[ u, Derivative[ 1][ c][ w]], Times[ v, Derivative[ 1][ d][ w]]]], Times[ Plus[ Times[ u, Derivative[ 1][ a][ w]], Times[ v, Derivative[ 1][ b][ w]]], Derivative[ 1][ f][ w]]], Plus[ Times[ c[ w], Derivative[ 1][ b][ w]], Times[ b[ w], Derivative[ 1][ c][ w]], Times[ -1, a[ w], Derivative[ 1][ d][ w]], Times[ Derivative[ 1][ a][ w], Derivative[ 1][ f][ w]], Times[ Derivative[ 1][ d][ w], Derivative[ 1][ f][ w]], Times[ a[ w], Derivative[ 2][ f][ w]], Times[ -2, Derivative[ 1][ f][ w], Derivative[ 2][ f][ w]], Times[ d[ w], Plus[ Times[ -1, Derivative[ 1][ a][ w]], Derivative[ 2][ f][ w]]]]], Times[ -1, Plus[ Times[ b[ w], c[ w]], Times[ Plus[ a[ w], Times[ -1, Derivative[ 1][ f][ w]]], Plus[ Times[ -1, d[ w]], Derivative[ 1][ f][ w]]]], Plus[ Times[ -1, u, d[ w], Derivative[ 2][ a][ w]], Times[ u, Derivative[ 1][ f][ w], Derivative[ 2][ a][ w]], Times[ -1, v, d[ w], Derivative[ 2][ b][ w]], Times[ v, Derivative[ 1][ f][ w], Derivative[ 2][ b][ w]], Times[ u, b[ w], Derivative[ 2][ c][ w]], Times[ v, b[ w], Derivative[ 2][ d][ w]], Times[ u, Derivative[ 1][ a][ w], Plus[ Times[ -1, Derivative[ 1][ d][ w]], Derivative[ 2][ f][ w]]], Times[ Derivative[ 1][ b][ w], Plus[ Times[ u, Derivative[ 1][ c][ w]], Times[ v, Derivative[ 2][ f][ w]]]]]]]]]]}]]], "Output"], Cell[BoxData[ \("Here is a way to find the riemann invarients of LEV1,2"\)], "Output"], Cell[BoxData[ InterpretationBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"-", FractionBox["1", RowBox[{ "2", " ", \((a[w]\^2 + 4\ b[w]\ c[w] - 2\ a[w]\ d[w] + d[w]\^2)\), " ", RowBox[{"(", RowBox[{\(b[w]\ c[w]\), "+", RowBox[{ RowBox[{"(", RowBox[{\(a[w]\), "-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{\(-d[w]\), "+", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "w", "]"}]}], ")"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{\((a[w] - d[w] + \@\(a[w]\^2 + 4\ b[w]\ c[w] - 2\ \ a[w]\ d[w] + d[w]\^2\))\), " ", RowBox[{"(", RowBox[{ RowBox[{\(c[w]\), " ", RowBox[{"(", RowBox[{ RowBox[{"u", " ", RowBox[{