x0star160.nb

(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     74452,       2192]
NotebookOptionsPosition[     70488,       2058]
NotebookOutlinePosition[     70831,       2073]
CellTagsIndexPosition[     70788,       2070]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f1", "=", 
    RowBox[{"q", "-", 
     RowBox[{"2", "*", 
      RowBox[{"q", "^", "3"}]}], "-", 
     RowBox[{"q", "^", "5"}], "+", 
     RowBox[{"2", "*", 
      RowBox[{"q", "^", "7"}]}], "+", 
     RowBox[{"q", "^", "9"}], "+", 
     RowBox[{"2", "*", 
      RowBox[{"q", "^", "13"}]}], "+", 
     RowBox[{"2", "*", 
      RowBox[{"q", "^", "15"}]}], "-", 
     RowBox[{"6", "*", 
      RowBox[{"q", "^", "17"}]}], "-", 
     RowBox[{"4", "*", 
      RowBox[{"q", "^", "19"}]}], "-", 
     RowBox[{"4", "*", 
      RowBox[{"q", "^", "21"}]}], "+", 
     RowBox[{"6", "*", 
      RowBox[{"q", "^", "23"}]}], "+", 
     RowBox[{"q", "^", "25"}], "+", 
     RowBox[{"4", "*", 
      RowBox[{"q", "^", "27"}]}], "+", 
     RowBox[{"6", "*", 
      RowBox[{"q", "^", "29"}]}]}]}], ";"}], " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"f3", "=", 
   RowBox[{"q", "+", 
    RowBox[{"2", "*", 
     RowBox[{"q", "^", "3"}]}], "-", 
    RowBox[{"q", "^", "5"}], "-", 
    RowBox[{"2", "*", 
     RowBox[{"q", "^", "7"}]}], "+", 
    RowBox[{"q", "^", "9"}], "+", 
    RowBox[{"2", "*", 
     RowBox[{"q", "^", "13"}]}], "-", 
    RowBox[{"2", "*", 
     RowBox[{"q", "^", "15"}]}], "-", 
    RowBox[{"6", "*", 
     RowBox[{"q", "^", "17"}]}], "+", 
    RowBox[{"4", "*", 
     RowBox[{"q", "^", "19"}]}], "-", 
    RowBox[{"4", "*", 
     RowBox[{"q", "^", "21"}]}], "-", 
    RowBox[{"6", "*", 
     RowBox[{"q", "^", "23"}]}], "+", 
    RowBox[{"q", "^", "25"}], "-", 
    RowBox[{"4", "*", 
     RowBox[{"q", "^", "27"}]}], "+", 
    RowBox[{"6", "*", 
     RowBox[{"q", "^", "29"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"f4", "=", 
   RowBox[{"q", "-", 
    RowBox[{"2", "*", 
     RowBox[{"q", "^", "3"}]}], "-", 
    RowBox[{"q", "^", "5"}], "-", 
    RowBox[{"2", "*", 
     RowBox[{"q", "^", "7"}]}], "+", 
    RowBox[{"q", "^", "9"}], "-", 
    RowBox[{"4", "*", 
     RowBox[{"q", "^", "11"}]}], "-", 
    RowBox[{"6", "*", 
     RowBox[{"q", "^", "13"}]}], "+", 
    RowBox[{"2", "*", 
     RowBox[{"q", "^", "15"}]}], "+", 
    RowBox[{"2", "*", 
     RowBox[{"q", "^", "17"}]}], "+", 
    RowBox[{"8", "*", 
     RowBox[{"q", "^", "19"}]}], "+", 
    RowBox[{"4", "*", 
     RowBox[{"q", "^", "21"}]}], "-", 
    RowBox[{"6", "*", 
     RowBox[{"q", "^", "23"}]}], "+", 
    RowBox[{"q", "^", "25"}], "+", 
    RowBox[{"4", "*", 
     RowBox[{"q", "^", "27"}]}], "-", 
    RowBox[{"2", "*", 
     RowBox[{"q", "^", "29"}]}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.907595485755*^9, 3.9075956206140003`*^9}, {
  3.9075959643970003`*^9, 3.907595966207*^9}, {3.907600174572*^9, 
  3.907600183672*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.9093369950081997`*^9, 3.9093371116582003`*^9}, {
   3.9093371566982*^9, 3.9093372123382*^9}, {3.9093372465882*^9, 
   3.9093372571282*^9}, 3.9093670614456*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"g1", "=", 
   RowBox[{
    RowBox[{"f1", "-", 
     RowBox[{"2", 
      RowBox[{"(", 
       RowBox[{"f1", "/.", 
        RowBox[{"q", "\[Rule]", 
         RowBox[{"q", "^", "2"}]}]}], ")"}]}]}], "//", "Expand"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.907496426975*^9, 3.907496441227*^9}, {3.907496485135*^9,
   3.907496508999*^9}, {3.907595664492*^9, 3.907595668887*^9}, {
  3.907595994172*^9, 3.9075960087320004`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"g11", "=", 
   RowBox[{
    RowBox[{"g1", "+", 
     RowBox[{"2", 
      RowBox[{"(", 
       RowBox[{"g1", "/.", 
        RowBox[{"q", "\[Rule]", 
         RowBox[{"q", "^", "2"}]}]}], ")"}]}]}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"g111", "=", 
   RowBox[{
    RowBox[{"g11", "+", 
     RowBox[{"2", 
      RowBox[{"(", 
       RowBox[{"g11", "/.", 
        RowBox[{"q", "\[Rule]", 
         RowBox[{"q", "^", "2"}]}]}], ")"}]}]}], "//", "Expand"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.907496426975*^9, 3.907496441227*^9}, {3.907496485135*^9,
   3.907496508999*^9}, {3.907595678142*^9, 3.907595722707*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"g12", "=", 
   RowBox[{
    RowBox[{"g1", "-", 
     RowBox[{"2", 
      RowBox[{"(", 
       RowBox[{"g1", "/.", 
        RowBox[{"q", "\[Rule]", 
         RowBox[{"q", "^", "2"}]}]}], ")"}]}]}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"g121", "=", 
   RowBox[{
    RowBox[{"g12", "-", 
     RowBox[{"2", 
      RowBox[{"(", 
       RowBox[{"g12", "/.", 
        RowBox[{"q", "\[Rule]", 
         RowBox[{"q", "^", "2"}]}]}], ")"}]}]}], "//", "Expand"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.907595738802*^9, 3.907595752632*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"g3", "=", 
   RowBox[{
    RowBox[{"f3", "-", 
     RowBox[{"2", 
      RowBox[{"(", 
       RowBox[{"f3", "/.", 
        RowBox[{"q", "\[Rule]", 
         RowBox[{"q", "^", "2"}]}]}], ")"}]}]}], "//", "Expand"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.907595819977*^9, 3.907595851822*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"m", "=", "29"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h1", "=", " ", 
   RowBox[{"Series", "[", 
    RowBox[{"g111", ",", 
     RowBox[{"{", 
      RowBox[{"q", ",", "0", ",", "m"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"h2", "=", " ", 
    RowBox[{"Series", "[", 
     RowBox[{"g121", ",", 
      RowBox[{"{", 
       RowBox[{"q", ",", "0", ",", "m"}], "}"}]}], "]"}]}], ";"}], 
  " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h3", "=", 
   RowBox[{"Series", "[", 
    RowBox[{"g3", ",", 
     RowBox[{"{", 
      RowBox[{"q", ",", "0", ",", "m"}], "}"}]}], "]"}]}], ";", 
  RowBox[{"h4", "=", 
   RowBox[{"Series", "[", 
    RowBox[{"f4", ",", 
     RowBox[{"{", 
      RowBox[{"q", ",", "0", ",", "m"}], "}"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.907496653381*^9, 3.907496685008*^9}, {3.907595778612*^9,
   3.907595794427*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"P2", "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "a1", ",", "a2", ",", "a3", ",", "a4", ",", "a5", ",", "a6", ",", "a7", 
     ",", "a8", ",", "a9", ",", "a10"}], "}"}], ".", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"x", "^", "2"}], ",", 
     RowBox[{"y", "^", "2"}], ",", 
     RowBox[{"z", "^", "2"}], ",", 
     RowBox[{"t", "^", "2"}], ",", " ", 
     RowBox[{"x", " ", "y"}], ",", " ", 
     RowBox[{"x", " ", "z"}], ",", " ", 
     RowBox[{"x", " ", "t"}], ",", " ", 
     RowBox[{"y", " ", "z"}], ",", " ", 
     RowBox[{"y", " ", "t"}], ",", 
     RowBox[{"z", " ", "t"}]}], "}"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"a4", " ", 
   SuperscriptBox["t", "2"]}], "+", 
  RowBox[{"a7", " ", "t", " ", "x"}], "+", 
  RowBox[{"a1", " ", 
   SuperscriptBox["x", "2"]}], "+", 
  RowBox[{"a9", " ", "t", " ", "y"}], "+", 
  RowBox[{"a5", " ", "x", " ", "y"}], "+", 
  RowBox[{"a2", " ", 
   SuperscriptBox["y", "2"]}], "+", 
  RowBox[{"a10", " ", "t", " ", "z"}], "+", 
  RowBox[{"a6", " ", "x", " ", "z"}], "+", 
  RowBox[{"a8", " ", "y", " ", "z"}], "+", 
  RowBox[{"a3", " ", 
   SuperscriptBox["z", "2"]}]}]], "Output",
 CellChangeTimes->{
  3.907595924962*^9, 3.907600194317*^9, {3.9076002356070004`*^9, 
   3.907600261657*^9}, 3.907600593467*^9, 3.9076066980422*^9, {
   3.9076067822472*^9, 3.9076067916222*^9}, {3.9076069043822002`*^9, 
   3.9076069237072*^9}, {3.9093370114982*^9, 3.9093371174082003`*^9}, {
   3.9093371763282003`*^9, 3.9093372170481997`*^9}, {3.9093372512782*^9, 
   3.9093372639282*^9}, 3.909338695502*^9, 3.909338874277*^9, 
   3.909339063737*^9, 3.909339247071*^9, 3.909339413851*^9, 
   3.9093394702460003`*^9, 3.909339519046*^9, 3.909339569231*^9, {
   3.909339610401*^9, 3.909339624256*^9}, 3.9093670681296*^9, 
   3.909367199869*^9, 3.909367310849*^9, 3.909411060802*^9, 3.909411351484*^9,
    3.909411408325*^9, {3.909411441257*^9, 3.909411474918*^9}, 
   3.909412595526*^9, 3.90941264079*^9, {3.9094127845290003`*^9, 
   3.909412807179*^9}, {3.909412865039*^9, 3.909412892289*^9}, 
   3.9094129420290003`*^9, 3.90947935176*^9, 3.916843590072*^9, 
   3.91684404323625*^9, 3.91684408652625*^9, 3.91684427419025*^9, 
   3.91684444854025*^9, 3.91684460434125*^9, 3.9168448537922497`*^9, 
   3.9168455774352503`*^9, 3.9168456597252502`*^9, 3.91684679954525*^9, 
   3.9168473527372503`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Q", "=", 
   RowBox[{"P2", "/.", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x", "\[Rule]", "h1"}], ",", 
      RowBox[{"y", "\[Rule]", "h2"}], ",", 
      RowBox[{"z", "\[Rule]", "h3"}], ",", 
      RowBox[{"t", "\[Rule]", "h4"}]}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"l", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Coefficient", "[", 
       RowBox[{"Q", ",", "q", ",", "i"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "2", ",", "14"}], "}"}]}], "]"}], "//", "Factor"}]}],
   ";", 
  RowBox[{"T", "=", 
   RowBox[{
    RowBox[{"Solve", "[", " ", 
     RowBox[{
      RowBox[{"l", "\[Equal]", 
       RowBox[{"{", 
        RowBox[{
        "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
         "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]}], 
      ",", 
      RowBox[{"{", 
       RowBox[{
       "a1", ",", "a2", ",", "a3", ",", "a4", ",", "a5", ",", "a6", ",", "a7",
         ",", "a8", ",", "a9", ",", "a10"}], "}"}]}], "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]}], "Input",
 CellChangeTimes->{{3.907600224292*^9, 3.907600257242*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"a2", "\[Rule]", 
    RowBox[{"-", 
     FractionBox["a1", "3"]}]}], ",", 
   RowBox[{"a3", "\[Rule]", 
    RowBox[{"-", 
     FractionBox[
      RowBox[{"8", " ", "a1"}], "3"]}]}], ",", 
   RowBox[{"a4", "\[Rule]", 
    RowBox[{
     RowBox[{"-", "16"}], " ", "a1"}]}], ",", 
   RowBox[{"a5", "\[Rule]", 
    RowBox[{"2", " ", "a1"}]}], ",", 
   RowBox[{"a6", "\[Rule]", 
    RowBox[{"12", " ", "a1"}]}], ",", 
   RowBox[{"a7", "\[Rule]", 
    FractionBox[
     RowBox[{"8", " ", "a1"}], "3"]}], ",", 
   RowBox[{"a8", "\[Rule]", 
    RowBox[{"4", " ", "a1"}]}], ",", 
   RowBox[{"a9", "\[Rule]", 
    RowBox[{"-", 
     FractionBox[
      RowBox[{"8", " ", "a1"}], "3"]}]}], ",", 
   RowBox[{"a10", "\[Rule]", "0"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.907595925762*^9, 3.907600194487*^9, {3.907600235757*^9, 
   3.907600261817*^9}, 3.9076005936070004`*^9, 3.9076066982922*^9, {
   3.9076067824072*^9, 3.9076067918022003`*^9}, {3.9076069045422*^9, 
   3.9076069238772*^9}, {3.9093370116582003`*^9, 3.9093371175782003`*^9}, {
   3.9093371764882*^9, 3.9093372171982*^9}, {3.9093372514282*^9, 
   3.9093372640782003`*^9}, 3.909338695797*^9, 3.909338874457*^9, 
   3.909339063957*^9, 3.9093392472609997`*^9, 3.9093394140360003`*^9, 
   3.909339470451*^9, 3.909339519211*^9, 3.9093395694110003`*^9, {
   3.909339610576*^9, 3.909339624481*^9}, 3.9093670683496*^9, 
   3.9093672000290003`*^9, 3.909367310939*^9, 3.909411060922*^9, 
   3.909411351652*^9, 3.909411408465*^9, {3.9094114414370003`*^9, 
   3.9094114750880003`*^9}, 3.9094125957060003`*^9, 3.9094126409300003`*^9, {
   3.909412784699*^9, 3.9094128073190002`*^9}, {3.909412865219*^9, 
   3.909412892449*^9}, 3.909412942189*^9, 3.90947935201*^9, 3.916843590251*^9,
    3.91684404345625*^9, 3.91684408667625*^9, 3.91684427436025*^9, 
   3.91684444882025*^9, 3.91684460454125*^9, 3.91684485403225*^9, 
   3.91684557761525*^9, 3.91684565980525*^9, 3.91684679966525*^9, 
   3.91684735286725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Q2", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{"P2", "/.", "T"}], ")"}], "/.", 
     RowBox[{"a1", "\[Rule]", "1"}]}], "//", "Factor"}], "//", 
   "Numerator"}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "48"}], " ", 
   SuperscriptBox["t", "2"]}], "+", 
  RowBox[{"8", " ", "t", " ", "x"}], "+", 
  RowBox[{"3", " ", 
   SuperscriptBox["x", "2"]}], "-", 
  RowBox[{"8", " ", "t", " ", "y"}], "+", 
  RowBox[{"6", " ", "x", " ", "y"}], "-", 
  SuperscriptBox["y", "2"], "+", 
  RowBox[{"36", " ", "x", " ", "z"}], "+", 
  RowBox[{"12", " ", "y", " ", "z"}], "-", 
  RowBox[{"8", " ", 
   SuperscriptBox["z", "2"]}]}]], "Output",
 CellChangeTimes->{
  3.9075959259519997`*^9, 3.907600194517*^9, {3.907600235797*^9, 
   3.9076002618570004`*^9}, 3.907600593637*^9, 3.9076066983622*^9, {
   3.9076067825372*^9, 3.9076067918322*^9}, {3.9076069045822*^9, 
   3.9076069239072*^9}, {3.9093370117282*^9, 3.9093371176082*^9}, {
   3.9093371765182*^9, 3.9093372172382*^9}, {3.9093372514681997`*^9, 
   3.9093372641182003`*^9}, 3.909338695867*^9, 3.909338874487*^9, 
   3.9093390639969997`*^9, 3.9093392472860003`*^9, 3.909339414066*^9, 
   3.909339470491*^9, 3.909339519241*^9, 3.909339569441*^9, {
   3.909339610611*^9, 3.9093396246210003`*^9}, 3.9093670683796*^9, 
   3.909367200134*^9, 3.909367310984*^9, 3.9094110609519997`*^9, 
   3.909411351692*^9, 3.909411408495*^9, {3.909411441477*^9, 
   3.9094114751280003`*^9}, 3.9094125957460003`*^9, 3.90941264096*^9, {
   3.9094127847390003`*^9, 3.909412807349*^9}, {3.909412865249*^9, 
   3.909412892479*^9}, 3.909412942229*^9, 3.90947935204*^9, 3.916843590299*^9,
    3.91684404349625*^9, 3.91684408670625*^9, 3.91684427439025*^9, 
   3.9168444488502502`*^9, 3.91684460458125*^9, 3.91684485407225*^9, 
   3.91684557765525*^9, 3.91684565983525*^9, 3.91684679969525*^9, 
   3.91684735290725*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"P1", "=", "Q2"}], ";"}]], "Input",
 CellChangeTimes->{{3.907600582297*^9, 3.9076005849820004`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Q2", "/.", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"x", "\[Rule]", "h1"}], ",", 
    RowBox[{"y", "\[Rule]", "h2"}], ",", 
    RowBox[{"z", "\[Rule]", "h3"}], ",", 
    RowBox[{"t", "\[Rule]", "h4"}]}], "}"}]}]], "Input"],

Cell[BoxData[
 InterpretationBox[
  SuperscriptBox[
   RowBox[{"O", "[", "q", "]"}], "31"],
  SeriesData[$CellContext`q, 0, {}, 31, 31, 1],
  Editable->False]], "Output",
 CellChangeTimes->{
  3.907595926052*^9, 3.9076001946070004`*^9, {3.907600235847*^9, 
   3.9076002618970003`*^9}, 3.9076005936870003`*^9, 3.9076066984622*^9, {
   3.9076067825972*^9, 3.9076067918922*^9}, {3.9076069046422*^9, 
   3.9076069239672003`*^9}, {3.9093370118282003`*^9, 3.9093371176782*^9}, {
   3.9093371765782003`*^9, 3.9093372172981997`*^9}, {3.9093372515282*^9, 
   3.9093372641782*^9}, 3.909338695967*^9, 3.909338874542*^9, 
   3.909339064047*^9, 3.909339247336*^9, 3.909339414116*^9, 3.909339470546*^9,
    3.909339519291*^9, 3.909339569491*^9, {3.9093396106610003`*^9, 
   3.909339624671*^9}, 3.9093670685095997`*^9, 3.909367200189*^9, 
   3.909367311064*^9, 3.909411060992*^9, 3.909411351762*^9, 
   3.9094114085550003`*^9, {3.9094114415369997`*^9, 3.909411475198*^9}, 
   3.9094125958059998`*^9, 3.90941264102*^9, {3.909412784799*^9, 
   3.909412807399*^9}, {3.9094128653190002`*^9, 3.9094128925290003`*^9}, 
   3.909412942289*^9, 3.9094793521400003`*^9, 3.9168435904049997`*^9, 
   3.9168440435862503`*^9, 3.91684408675625*^9, 3.91684427444025*^9, 
   3.91684444890025*^9, 3.9168446046412497`*^9, 3.91684485414225*^9, 
   3.91684557777525*^9, 3.91684565988525*^9, 3.91684679973525*^9, 
   3.91684735295725*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.907600574462*^9, 3.907600575902*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{"P1", "-", 
     RowBox[{"(", 
      RowBox[{"P1", "/.", 
       RowBox[{"x", "->", 
        RowBox[{"-", "x"}]}]}], ")"}]}], ",", 
    RowBox[{"P1", "-", 
     RowBox[{"(", 
      RowBox[{"P1", "/.", 
       RowBox[{"y", "->", 
        RowBox[{"-", "y"}]}]}], ")"}]}], ",", 
    RowBox[{"P1", "-", 
     RowBox[{"(", 
      RowBox[{"P1", "/.", 
       RowBox[{"z", "->", 
        RowBox[{"-", "z"}]}]}], ")"}]}], ",", 
    RowBox[{"P1", "-", 
     RowBox[{"(", 
      RowBox[{"P1", "/.", 
       RowBox[{"t", "->", 
        RowBox[{"-", "t"}]}]}], ")"}]}]}], "}"}], "//", "Factor"}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"4", " ", "x", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", "t"}], "+", 
      RowBox[{"3", " ", "y"}], "+", 
      RowBox[{"18", " ", "z"}]}], ")"}]}], ",", 
   RowBox[{
    RowBox[{"-", "4"}], " ", "y", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", "t"}], "-", 
      RowBox[{"3", " ", "x"}], "-", 
      RowBox[{"6", " ", "z"}]}], ")"}]}], ",", 
   RowBox[{"24", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"3", " ", "x"}], "+", "y"}], ")"}], " ", "z"}], ",", 
   RowBox[{"16", " ", "t", " ", 
    RowBox[{"(", 
     RowBox[{"x", "-", "y"}], ")"}]}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.9076005937469997`*^9, 3.9076066985222*^9, {3.9076067826572*^9, 
   3.9076067919622*^9}, {3.9076069047022*^9, 3.9076069240372*^9}, {
   3.9093370118882*^9, 3.9093371177382*^9}, {3.9093371766482*^9, 
   3.9093372173582*^9}, {3.9093372515882*^9, 3.9093372642382*^9}, 
   3.9093386960220003`*^9, 3.909338874597*^9, 3.909339064097*^9, 
   3.909339247391*^9, 3.909339414171*^9, 3.909339470606*^9, 3.909339519351*^9,
    3.909339569541*^9, {3.909339610716*^9, 3.909339624726*^9}, 
   3.9093670685596*^9, 3.9093672002390003`*^9, 3.909367311119*^9, 
   3.909411061032*^9, 3.909411351822*^9, 3.909411408615*^9, {
   3.9094114416070004`*^9, 3.909411475258*^9}, 3.909412595876*^9, 
   3.9094126410699997`*^9, {3.909412784869*^9, 3.909412807459*^9}, {
   3.909412865399*^9, 3.909412892579*^9}, 3.909412942359*^9, 3.90947935219*^9,
    3.9168435904630003`*^9, 3.91684404373625*^9, 3.91684408681625*^9, 
   3.91684427450025*^9, 3.9168444489502497`*^9, 3.91684460470125*^9, 
   3.91684485420225*^9, 3.91684557783525*^9, 3.9168456599352503`*^9, 
   3.9168467997852497`*^9, 3.91684735301725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{"P1", "-", 
     RowBox[{"(", 
      RowBox[{"P1", "/.", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"x", "->", 
          RowBox[{"-", "x"}]}], ",", " ", 
         RowBox[{"y", "->", 
          RowBox[{"-", "y"}]}]}], "}"}]}], ")"}]}], ",", 
    RowBox[{"P1", "-", 
     RowBox[{"(", 
      RowBox[{"P1", "/.", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"z", "->", 
          RowBox[{"-", "z"}]}], ",", " ", 
         RowBox[{"t", "->", 
          RowBox[{"-", "t"}]}]}], "}"}]}], ")"}]}]}], "}"}], "//", 
  "Factor"}]], "Input",
 CellChangeTimes->{{3.9076045158599997`*^9, 3.907604547815*^9}, {
  3.90760458121*^9, 3.907604650225*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"8", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"2", " ", "t", " ", "x"}], "-", 
      RowBox[{"2", " ", "t", " ", "y"}], "+", 
      RowBox[{"9", " ", "x", " ", "z"}], "+", 
      RowBox[{"3", " ", "y", " ", "z"}]}], ")"}]}], ",", 
   RowBox[{"8", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"2", " ", "t", " ", "x"}], "-", 
      RowBox[{"2", " ", "t", " ", "y"}], "+", 
      RowBox[{"9", " ", "x", " ", "z"}], "+", 
      RowBox[{"3", " ", "y", " ", "z"}]}], ")"}]}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.90760465417*^9, 3.9076066985522003`*^9, {3.9076067826872*^9, 
   3.9076067919922*^9}, {3.9076069047422*^9, 3.9076069240772*^9}, {
   3.9093370119282*^9, 3.9093371177782*^9}, {3.9093371766782*^9, 
   3.9093372173882*^9}, {3.9093372516182003`*^9, 3.9093372642782*^9}, 
   3.909338696057*^9, 3.909338874627*^9, 3.909339064127*^9, 3.909339247421*^9,
    3.909339414201*^9, 3.909339470646*^9, 3.9093395194560003`*^9, 
   3.909339569576*^9, {3.9093396107609997`*^9, 3.9093396247609997`*^9}, 
   3.9093670685996*^9, 3.909367200274*^9, 3.909367311149*^9, 
   3.9094110610620003`*^9, 3.909411351862*^9, 3.909411408635*^9, {
   3.9094114416470003`*^9, 3.9094114752980003`*^9}, 3.909412595926*^9, 
   3.9094126411*^9, {3.909412784909*^9, 3.909412807499*^9}, {
   3.909412865439*^9, 3.909412892619*^9}, 3.909412942399*^9, 
   3.9094793522200003`*^9, 3.916843590488*^9, 3.9168440437962503`*^9, 
   3.91684408684625*^9, 3.91684427453025*^9, 3.91684444899025*^9, 
   3.91684460473125*^9, 3.91684485424225*^9, 3.91684557788525*^9, 
   3.91684565996525*^9, 3.91684679981525*^9, 3.91684735304725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{"P1", "-", 
     RowBox[{"(", 
      RowBox[{"P1", "/.", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"x", "->", 
          RowBox[{"-", "x"}]}], ",", " ", 
         RowBox[{"t", "->", 
          RowBox[{"-", "t"}]}]}], "}"}]}], ")"}]}], ",", 
    RowBox[{"P1", "-", 
     RowBox[{"(", 
      RowBox[{"P1", "/.", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"y", "->", 
          RowBox[{"-", "y"}]}], ",", " ", 
         RowBox[{"t", "->", 
          RowBox[{"-", "t"}]}]}], "}"}]}], ")"}]}]}], "}"}], "//", 
  "Factor"}]], "Input",
 CellChangeTimes->{{3.90760475462*^9, 3.9076047647650003`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "4"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", "t", " ", "y"}], "-", 
      RowBox[{"3", " ", "x", " ", "y"}], "-", 
      RowBox[{"18", " ", "x", " ", "z"}]}], ")"}]}], ",", 
   RowBox[{"4", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", "t", " ", "x"}], "+", 
      RowBox[{"3", " ", "x", " ", "y"}], "+", 
      RowBox[{"6", " ", "y", " ", "z"}]}], ")"}]}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.90760476937*^9, 3.9076066985922003`*^9, {3.9076067827272*^9, 
   3.9076067920322*^9}, {3.9076069047822*^9, 3.9076069241071997`*^9}, {
   3.9093370119681997`*^9, 3.9093371178082*^9}, {3.9093371767181997`*^9, 
   3.9093372174282*^9}, {3.9093372516582003`*^9, 3.9093372643082*^9}, 
   3.9093386961070004`*^9, 3.9093388746619997`*^9, 3.909339064157*^9, 
   3.9093392474560003`*^9, 3.909339414236*^9, 3.9093394706809998`*^9, 
   3.909339519616*^9, 3.909339569606*^9, {3.9093396108710003`*^9, 
   3.909339624796*^9}, 3.9093670686296*^9, 3.909367200304*^9, 
   3.909367311184*^9, 3.909411061092*^9, 3.909411351902*^9, 
   3.909411408675*^9, {3.9094114416870003`*^9, 3.909411475348*^9}, 
   3.909412595966*^9, 3.9094126411400003`*^9, {3.909412784949*^9, 
   3.909412807539*^9}, {3.909412865479*^9, 3.909412892649*^9}, 
   3.909412942439*^9, 3.90947935225*^9, 3.91684359052*^9, 3.91684404384625*^9,
    3.9168440868762503`*^9, 3.9168442745602503`*^9, 3.91684444903025*^9, 
   3.91684460477125*^9, 3.91684485428225*^9, 3.91684557791525*^9, 
   3.91684565998525*^9, 3.9168467999352503`*^9, 3.91684735307725*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.9076064977022*^9, 3.9076065120872*^9}, 
   3.9076066050422*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PP1", "=", 
  RowBox[{
   RowBox[{"P1", "/.", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x", "\[Rule]", 
       RowBox[{
        RowBox[{"aa1", " ", "x"}], "+", " ", 
        RowBox[{"aa2", " ", "y"}]}]}], ",", 
      RowBox[{"y", "\[Rule]", " ", 
       RowBox[{
        RowBox[{"bb1", " ", "x"}], "+", " ", 
        RowBox[{"bb2", " ", "y"}]}]}], ",", 
      RowBox[{"t", "\[Rule]", 
       RowBox[{"-", "t"}]}]}], "}"}]}], "//", "Expand"}]}]], "Input",
 CellChangeTimes->{{3.9076067277672*^9, 3.9076067317072*^9}, {
   3.9076067722372*^9, 3.9076067756822*^9}, {3.9076068970772*^9, 
   3.9076069186422*^9}, 3.909367189199*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "48"}], " ", 
   SuperscriptBox["t", "2"]}], "-", 
  RowBox[{"8", " ", "aa1", " ", "t", " ", "x"}], "+", 
  RowBox[{"8", " ", "bb1", " ", "t", " ", "x"}], "+", 
  RowBox[{"3", " ", 
   SuperscriptBox["aa1", "2"], " ", 
   SuperscriptBox["x", "2"]}], "+", 
  RowBox[{"6", " ", "aa1", " ", "bb1", " ", 
   SuperscriptBox["x", "2"]}], "-", 
  RowBox[{
   SuperscriptBox["bb1", "2"], " ", 
   SuperscriptBox["x", "2"]}], "-", 
  RowBox[{"8", " ", "aa2", " ", "t", " ", "y"}], "+", 
  RowBox[{"8", " ", "bb2", " ", "t", " ", "y"}], "+", 
  RowBox[{"6", " ", "aa1", " ", "aa2", " ", "x", " ", "y"}], "+", 
  RowBox[{"6", " ", "aa2", " ", "bb1", " ", "x", " ", "y"}], "+", 
  RowBox[{"6", " ", "aa1", " ", "bb2", " ", "x", " ", "y"}], "-", 
  RowBox[{"2", " ", "bb1", " ", "bb2", " ", "x", " ", "y"}], "+", 
  RowBox[{"3", " ", 
   SuperscriptBox["aa2", "2"], " ", 
   SuperscriptBox["y", "2"]}], "+", 
  RowBox[{"6", " ", "aa2", " ", "bb2", " ", 
   SuperscriptBox["y", "2"]}], "-", 
  RowBox[{
   SuperscriptBox["bb2", "2"], " ", 
   SuperscriptBox["y", "2"]}], "+", 
  RowBox[{"36", " ", "aa1", " ", "x", " ", "z"}], "+", 
  RowBox[{"12", " ", "bb1", " ", "x", " ", "z"}], "+", 
  RowBox[{"36", " ", "aa2", " ", "y", " ", "z"}], "+", 
  RowBox[{"12", " ", "bb2", " ", "y", " ", "z"}], "-", 
  RowBox[{"8", " ", 
   SuperscriptBox["z", "2"]}]}]], "Output",
 CellChangeTimes->{
  3.9076066986622*^9, {3.9076067827872*^9, 3.9076067921022*^9}, {
   3.9076069048422003`*^9, 3.9076069241772003`*^9}, {3.9093370120382*^9, 
   3.9093371179582*^9}, {3.9093371767882*^9, 3.9093372174982*^9}, {
   3.9093372517181997`*^9, 3.9093372643782*^9}, 3.909338696382*^9, 
   3.909338874722*^9, 3.909339064217*^9, 3.9093392475109997`*^9, 
   3.909339414291*^9, 3.909339470736*^9, 3.909339519686*^9, 
   3.909339569756*^9, {3.909339610936*^9, 3.9093396249709997`*^9}, 
   3.9093670686896*^9, 3.909367200359*^9, 3.9093673112390003`*^9, 
   3.909411061132*^9, 3.909411351972*^9, 3.909411408725*^9, {
   3.909411441757*^9, 3.909411475418*^9}, 3.909412596043*^9, 
   3.9094126412*^9, {3.909412785019*^9, 3.909412807599*^9}, {
   3.909412865549*^9, 3.909412892699*^9}, 3.909412942509*^9, 3.90947935231*^9,
    3.916843590574*^9, 3.9168440439362497`*^9, 3.9168440869362497`*^9, 
   3.9168442746202497`*^9, 3.91684444909025*^9, 3.91684460484125*^9, 
   3.91684485435225*^9, 3.91684557796525*^9, 3.9168456600352497`*^9, 
   3.91684679998525*^9, 3.91684735313725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"PP1", "-", 
      RowBox[{"(", 
       RowBox[{"PP1", "/.", 
        RowBox[{"x", "\[Rule]", 
         RowBox[{"-", "x"}]}]}], ")"}]}], "//", "Expand"}], ")"}], "/", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "2"}], " ", "x"}], ")"}]}], "//", "Factor"}]], "Input"],

Cell[BoxData[
 RowBox[{"2", " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"4", " ", "aa1", " ", "t"}], "-", 
    RowBox[{"4", " ", "bb1", " ", "t"}], "-", 
    RowBox[{"3", " ", "aa1", " ", "aa2", " ", "y"}], "-", 
    RowBox[{"3", " ", "aa2", " ", "bb1", " ", "y"}], "-", 
    RowBox[{"3", " ", "aa1", " ", "bb2", " ", "y"}], "+", 
    RowBox[{"bb1", " ", "bb2", " ", "y"}], "-", 
    RowBox[{"18", " ", "aa1", " ", "z"}], "-", 
    RowBox[{"6", " ", "bb1", " ", "z"}]}], ")"}]}]], "Output",
 CellChangeTimes->{
  3.9076066986921997`*^9, {3.9076067828272*^9, 3.9076067921322002`*^9}, {
   3.9076069048822002`*^9, 3.9076069242072*^9}, {3.9093370120782003`*^9, 
   3.9093371179982*^9}, {3.9093371768182*^9, 3.9093372175282*^9}, {
   3.9093372517581997`*^9, 3.9093372644182*^9}, 3.909338696422*^9, 
   3.909338874757*^9, 3.9093390642469997`*^9, 3.909339247546*^9, 
   3.909339414326*^9, 3.909339470771*^9, 3.9093395197209997`*^9, 
   3.9093395697860003`*^9, {3.9093396109709997`*^9, 3.909339625006*^9}, 
   3.9093670687296*^9, 3.909367200394*^9, 3.909367311269*^9, 
   3.9094110611619997`*^9, 3.909411352012*^9, 3.9094114087650003`*^9, {
   3.909411441797*^9, 3.909411475458*^9}, 3.909412596083*^9, 
   3.90941264123*^9, {3.909412785059*^9, 3.909412807639*^9}, {
   3.9094128655889997`*^9, 3.909412892729*^9}, 3.909412942549*^9, 
   3.90947935234*^9, 3.916843590601*^9, 3.91684404398625*^9, 
   3.91684408696625*^9, 3.91684427465025*^9, 3.9168444491202497`*^9, 
   3.91684460488125*^9, 3.91684485439225*^9, 3.9168455779952497`*^9, 
   3.91684566006525*^9, 3.91684680010525*^9, 3.91684735317725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ll", "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Coefficient", "[", 
      RowBox[{"%", ",", "y", ",", "1"}], "]"}], ",", " ", 
     RowBox[{"Coefficient", "[", 
      RowBox[{"%", ",", "z", ",", "1"}], "]"}], ",", " ", 
     RowBox[{"Coefficient", "[", 
      RowBox[{"%", ",", "t", ",", "1"}], "]"}]}], "}"}], "//", 
   "Factor"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "2"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"3", " ", "aa1", " ", "aa2"}], "+", 
      RowBox[{"3", " ", "aa2", " ", "bb1"}], "+", 
      RowBox[{"3", " ", "aa1", " ", "bb2"}], "-", 
      RowBox[{"bb1", " ", "bb2"}]}], ")"}]}], ",", 
   RowBox[{
    RowBox[{"-", "12"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"3", " ", "aa1"}], "+", "bb1"}], ")"}]}], ",", 
   RowBox[{"8", " ", 
    RowBox[{"(", 
     RowBox[{"aa1", "-", "bb1"}], ")"}]}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.9076066987321997`*^9, {3.9076067828672*^9, 3.9076067921722*^9}, {
   3.9076069049222*^9, 3.9076069242472*^9}, {3.9093370122382*^9, 
   3.9093371180382*^9}, {3.9093371768582*^9, 3.9093372175682*^9}, {
   3.9093372517981997`*^9, 3.9093372644582*^9}, 3.909338696462*^9, 
   3.909338874792*^9, 3.909339064282*^9, 3.909339247576*^9, 3.909339414391*^9,
    3.9093394708059998`*^9, 3.9093395197609997`*^9, 3.909339569821*^9, {
   3.9093396110109997`*^9, 3.909339625046*^9}, 3.9093670689196*^9, 
   3.909367200424*^9, 3.909367311304*^9, 3.909411061192*^9, 3.909411352052*^9,
    3.909411408795*^9, {3.909411441837*^9, 3.909411475498*^9}, 
   3.909412596123*^9, 3.90941264126*^9, {3.909412785099*^9, 
   3.909412807679*^9}, {3.9094128657390003`*^9, 3.909412892769*^9}, 
   3.9094129425889997`*^9, 3.90947935237*^9, 3.9168435906289997`*^9, 
   3.91684404402625*^9, 3.91684408700625*^9, 3.91684427469025*^9, 
   3.91684444915025*^9, 3.91684460491125*^9, 3.9168448544622498`*^9, 
   3.9168455780352497`*^9, 3.91684566009525*^9, 3.91684680013525*^9, 
   3.91684735320725*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"%", "Idem", " ", "t"}], "\[Rule]", 
   RowBox[{"-", "t"}]}], ",", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"or", " ", "z"}], "\[Rule]", 
    RowBox[{
     RowBox[{"-", "z"}], " ", "need", " ", "impose", " ", "aa1"}]}], "=", 
   RowBox[{"bb1", "=", "0"}]}], ",", " ", 
  RowBox[{"thus", " ", "no", " ", "invertible", " ", "matrix"}]}]], "Input",
 CellChangeTimes->{{3.9076069323772*^9, 3.9076069685122004`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PP1", "=", 
  RowBox[{
   RowBox[{"P1", "/.", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x", "\[Rule]", 
       RowBox[{
        RowBox[{"aa1", " ", "x"}], "+", " ", 
        RowBox[{"aa2", " ", "y"}]}]}], ",", 
      RowBox[{"y", "\[Rule]", " ", 
       RowBox[{
        RowBox[{"bb1", " ", "x"}], "+", " ", 
        RowBox[{"bb2", " ", "y"}]}]}], ",", 
      RowBox[{"z", "\[Rule]", 
       RowBox[{"-", "z"}]}]}], "}"}]}], "//", "Expand"}]}]], "Input",
 CellChangeTimes->{3.909338599145*^9, {3.909339050192*^9, 3.909339052547*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "48"}], " ", 
   SuperscriptBox["t", "2"]}], "+", 
  RowBox[{"8", " ", "aa1", " ", "t", " ", "x"}], "-", 
  RowBox[{"8", " ", "bb1", " ", "t", " ", "x"}], "+", 
  RowBox[{"3", " ", 
   SuperscriptBox["aa1", "2"], " ", 
   SuperscriptBox["x", "2"]}], "+", 
  RowBox[{"6", " ", "aa1", " ", "bb1", " ", 
   SuperscriptBox["x", "2"]}], "-", 
  RowBox[{
   SuperscriptBox["bb1", "2"], " ", 
   SuperscriptBox["x", "2"]}], "+", 
  RowBox[{"8", " ", "aa2", " ", "t", " ", "y"}], "-", 
  RowBox[{"8", " ", "bb2", " ", "t", " ", "y"}], "+", 
  RowBox[{"6", " ", "aa1", " ", "aa2", " ", "x", " ", "y"}], "+", 
  RowBox[{"6", " ", "aa2", " ", "bb1", " ", "x", " ", "y"}], "+", 
  RowBox[{"6", " ", "aa1", " ", "bb2", " ", "x", " ", "y"}], "-", 
  RowBox[{"2", " ", "bb1", " ", "bb2", " ", "x", " ", "y"}], "+", 
  RowBox[{"3", " ", 
   SuperscriptBox["aa2", "2"], " ", 
   SuperscriptBox["y", "2"]}], "+", 
  RowBox[{"6", " ", "aa2", " ", "bb2", " ", 
   SuperscriptBox["y", "2"]}], "-", 
  RowBox[{
   SuperscriptBox["bb2", "2"], " ", 
   SuperscriptBox["y", "2"]}], "-", 
  RowBox[{"36", " ", "aa1", " ", "x", " ", "z"}], "-", 
  RowBox[{"12", " ", "bb1", " ", "x", " ", "z"}], "-", 
  RowBox[{"36", " ", "aa2", " ", "y", " ", "z"}], "-", 
  RowBox[{"12", " ", "bb2", " ", "y", " ", "z"}], "-", 
  RowBox[{"8", " ", 
   SuperscriptBox["z", "2"]}]}]], "Output",
 CellChangeTimes->{
  3.909338696597*^9, 3.9093388748570004`*^9, 3.909339064337*^9, 
   3.909339247791*^9, 3.909339414521*^9, 3.909339470866*^9, 3.909339519821*^9,
    3.909339569891*^9, {3.909339611071*^9, 3.909339625106*^9}, 
   3.9093670690696*^9, 3.909367200484*^9, 3.909367311359*^9, 
   3.909411061252*^9, 3.9094113521219997`*^9, 3.909411408855*^9, {
   3.909411442007*^9, 3.909411475678*^9}, 3.909412596193*^9, 
   3.9094126413199997`*^9, {3.9094127851689997`*^9, 3.909412807749*^9}, {
   3.909412865809*^9, 3.9094128928190002`*^9}, 3.9094129426689997`*^9, 
   3.90947935248*^9, 3.916843590748*^9, 3.91684404409625*^9, 
   3.91684408705625*^9, 3.91684427474025*^9, 3.91684444921025*^9, 
   3.91684460498125*^9, 3.9168448545422497`*^9, 3.91684557823525*^9, 
   3.91684566013525*^9, 3.9168468001852503`*^9, 3.91684735335725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"PP1", "-", 
      RowBox[{"(", 
       RowBox[{"PP1", "/.", 
        RowBox[{"x", "\[Rule]", 
         RowBox[{"-", "x"}]}]}], ")"}]}], "//", "Expand"}], ")"}], "/", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "2"}], " ", "x"}], ")"}]}], "//", "Factor"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "2"}], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"4", " ", "aa1", " ", "t"}], "-", 
    RowBox[{"4", " ", "bb1", " ", "t"}], "+", 
    RowBox[{"3", " ", "aa1", " ", "aa2", " ", "y"}], "+", 
    RowBox[{"3", " ", "aa2", " ", "bb1", " ", "y"}], "+", 
    RowBox[{"3", " ", "aa1", " ", "bb2", " ", "y"}], "-", 
    RowBox[{"bb1", " ", "bb2", " ", "y"}], "-", 
    RowBox[{"18", " ", "aa1", " ", "z"}], "-", 
    RowBox[{"6", " ", "bb1", " ", "z"}]}], ")"}]}]], "Output",
 CellChangeTimes->{
  3.909367311389*^9, 3.909411061282*^9, 3.909411352152*^9, 
   3.909411408885*^9, {3.909411442047*^9, 3.909411475708*^9}, 
   3.909412596363*^9, 3.9094126413599997`*^9, {3.909412785209*^9, 
   3.909412807789*^9}, {3.909412865849*^9, 3.909412892859*^9}, 
   3.909412942709*^9, 3.90947935269*^9, 3.916843590888*^9, 
   3.91684404413625*^9, 3.91684408709625*^9, 3.91684427478025*^9, 
   3.91684444925025*^9, 3.91684460502125*^9, 3.91684485458225*^9, 
   3.91684557826525*^9, 3.91684566016525*^9, 3.91684680021525*^9, 
   3.91684735339725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ll", "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Coefficient", "[", 
      RowBox[{"%", ",", "y", ",", "1"}], "]"}], ",", " ", 
     RowBox[{"Coefficient", "[", 
      RowBox[{"%", ",", "z", ",", "1"}], "]"}], ",", " ", 
     RowBox[{"Coefficient", "[", 
      RowBox[{"%", ",", "t", ",", "1"}], "]"}]}], "}"}], "//", 
   "Factor"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "2"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"3", " ", "aa1", " ", "aa2"}], "+", 
      RowBox[{"3", " ", "aa2", " ", "bb1"}], "+", 
      RowBox[{"3", " ", "aa1", " ", "bb2"}], "-", 
      RowBox[{"bb1", " ", "bb2"}]}], ")"}]}], ",", 
   RowBox[{"12", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"3", " ", "aa1"}], "+", "bb1"}], ")"}]}], ",", 
   RowBox[{
    RowBox[{"-", "8"}], " ", 
    RowBox[{"(", 
     RowBox[{"aa1", "-", "bb1"}], ")"}]}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.909338838292*^9, 3.909338874892*^9, 3.9093390643719997`*^9, 
   3.909339247826*^9, 3.909339414566*^9, 3.909339470906*^9, 3.909339519966*^9,
    3.909339569926*^9, {3.909339611111*^9, 3.909339625146*^9}, 
   3.9093670691096*^9, 3.909367200554*^9, 3.909367311479*^9, 
   3.9094110613120003`*^9, 3.909411352292*^9, 3.909411408915*^9, {
   3.909411442087*^9, 3.909411475748*^9}, 3.909412596403*^9, 
   3.9094126413900003`*^9, {3.9094127853789997`*^9, 3.909412807829*^9}, {
   3.909412865889*^9, 3.909412892889*^9}, 3.909412942749*^9, 3.90947935291*^9,
    3.916843590918*^9, 3.91684404417625*^9, 3.91684408713625*^9, 
   3.9168442748102503`*^9, 3.91684444930025*^9, 3.91684460508125*^9, 
   3.9168448546122503`*^9, 3.91684557830525*^9, 3.91684566027525*^9, 
   3.9168468002452497`*^9, 3.91684735342725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"PP1", "-", 
      RowBox[{"(", 
       RowBox[{"PP1", "/.", 
        RowBox[{"y", "\[Rule]", 
         RowBox[{"-", "y"}]}]}], ")"}]}], "//", "Expand"}], ")"}], "/", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "2"}], " ", "y"}], ")"}]}], "//", "Factor"}]], "Input",
 CellChangeTimes->{{3.909339233009*^9, 3.909339241619*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "2"}], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"4", " ", "aa2", " ", "t"}], "-", 
    RowBox[{"4", " ", "bb2", " ", "t"}], "+", 
    RowBox[{"3", " ", "aa1", " ", "aa2", " ", "x"}], "+", 
    RowBox[{"3", " ", "aa2", " ", "bb1", " ", "x"}], "+", 
    RowBox[{"3", " ", "aa1", " ", "bb2", " ", "x"}], "-", 
    RowBox[{"bb1", " ", "bb2", " ", "x"}], "-", 
    RowBox[{"18", " ", "aa2", " ", "z"}], "-", 
    RowBox[{"6", " ", "bb2", " ", "z"}]}], ")"}]}]], "Output",
 CellChangeTimes->{
  3.909338696632*^9, 3.9093388749370003`*^9, 3.909339064407*^9, 
   3.909339247861*^9, 3.909339414606*^9, 3.909339470946*^9, 3.909339520001*^9,
    3.909339570086*^9, {3.909339611266*^9, 3.9093396251809998`*^9}, 
   3.9093670691496*^9, 3.9093672005889997`*^9, 3.909367311514*^9, 
   3.909411061332*^9, 3.909411352332*^9, 3.9094114089449997`*^9, {
   3.909411442127*^9, 3.909411475788*^9}, 3.909412596443*^9, 
   3.90941264142*^9, {3.9094127854189997`*^9, 3.909412807869*^9}, {
   3.9094128659189997`*^9, 3.9094128929189997`*^9}, 3.909412942889*^9, 
   3.90947935295*^9, 3.916843590955*^9, 3.91684404421625*^9, 
   3.91684408716625*^9, 3.9168442748502502`*^9, 3.91684444933025*^9, 
   3.91684460512125*^9, 3.9168448546522503`*^9, 3.91684557833525*^9, 
   3.91684566029525*^9, 3.91684680026525*^9, 3.91684735346725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{"P1", "-", 
     RowBox[{"(", 
      RowBox[{"P1", "/.", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"x", "->", 
          RowBox[{
           RowBox[{"-", "x"}], "-", 
           RowBox[{"2", "*", "y"}]}]}], ",", " ", 
         RowBox[{"z", "->", 
          RowBox[{"-", "z"}]}]}], "}"}]}], ")"}]}], ",", 
    RowBox[{"P1", "-", 
     RowBox[{"(", 
      RowBox[{"P1", "/.", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"x", "->", 
          RowBox[{
           RowBox[{"-", "x"}], "+", 
           RowBox[{"b", "*", "y"}]}]}], ",", " ", 
         RowBox[{"t", "->", 
          RowBox[{"-", "t"}]}]}], "}"}]}], ")"}]}]}], "}"}], "//", 
  "Factor"}]], "Input",
 CellChangeTimes->{{3.9093386540150003`*^9, 3.909338665115*^9}, {
  3.909339396651*^9, 3.909339407311*^9}, {3.909339449691*^9, 
  3.9093394643059998`*^9}, {3.909339509896*^9, 3.909339512376*^9}, {
  3.909339563066*^9, 3.909339563606*^9}, {3.909339602071*^9, 
  3.909339618901*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"16", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"t", " ", "x"}], "+", 
      RowBox[{"t", " ", "y"}], "-", 
      RowBox[{"3", " ", "y", " ", "z"}]}], ")"}]}], ",", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "16"}], " ", "t", " ", "y"}], "+", 
    RowBox[{"8", " ", "b", " ", "t", " ", "y"}], "+", 
    RowBox[{"12", " ", "x", " ", "y"}], "+", 
    RowBox[{"6", " ", "b", " ", "x", " ", "y"}], "-", 
    RowBox[{"6", " ", "b", " ", 
     SuperscriptBox["y", "2"]}], "-", 
    RowBox[{"3", " ", 
     SuperscriptBox["b", "2"], " ", 
     SuperscriptBox["y", "2"]}], "+", 
    RowBox[{"72", " ", "x", " ", "z"}], "-", 
    RowBox[{"36", " ", "b", " ", "y", " ", "z"}]}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.909338696667*^9, 3.909338875047*^9, 3.9093390644370003`*^9, 
   3.909339247891*^9, 3.9093394147060003`*^9, 3.909339470976*^9, 
   3.909339520026*^9, 3.909339570116*^9, {3.909339611296*^9, 
   3.909339625211*^9}, 3.9093670691896*^9, 3.909367200619*^9, 
   3.909367311549*^9, 3.909411061362*^9, 3.909411352362*^9, 
   3.909411408975*^9, {3.909411442167*^9, 3.909411475828*^9}, 
   3.909412596493*^9, 3.90941264146*^9, {3.909412785459*^9, 
   3.909412807909*^9}, {3.909412865969*^9, 3.909412892959*^9}, 
   3.909412942929*^9, 3.90947935298*^9, 3.916843590991*^9, 
   3.91684404425625*^9, 3.91684408719625*^9, 3.91684427488025*^9, 
   3.9168444493702497`*^9, 3.91684460516125*^9, 3.9168448546922503`*^9, 
   3.91684557836525*^9, 3.9168456603252497`*^9, 3.91684680029525*^9, 
   3.91684735349725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PP1", "=", 
  RowBox[{
   RowBox[{"P1", "/.", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x", "\[Rule]", 
       RowBox[{
        RowBox[{"aa1", " ", "x"}], "+", " ", 
        RowBox[{"aa2", " ", "y"}]}]}], ",", 
      RowBox[{"y", "\[Rule]", " ", 
       RowBox[{
        RowBox[{"bb1", " ", "x"}], "+", " ", 
        RowBox[{"bb2", " ", "y"}]}]}], ",", 
      RowBox[{"z", "\[Rule]", "z"}], ",", 
      RowBox[{"t", "\[Rule]", 
       RowBox[{"-", "t"}]}]}], "}"}]}], "//", "Expand"}]}]], "Input",
 CellChangeTimes->{{3.909410977612*^9, 3.909410986602*^9}, 3.909411402771*^9, 
   3.909412583841*^9, {3.909412858209*^9, 3.9094128611289997`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "48"}], " ", 
   SuperscriptBox["t", "2"]}], "-", 
  RowBox[{"8", " ", "aa1", " ", "t", " ", "x"}], "+", 
  RowBox[{"8", " ", "bb1", " ", "t", " ", "x"}], "+", 
  RowBox[{"3", " ", 
   SuperscriptBox["aa1", "2"], " ", 
   SuperscriptBox["x", "2"]}], "+", 
  RowBox[{"6", " ", "aa1", " ", "bb1", " ", 
   SuperscriptBox["x", "2"]}], "-", 
  RowBox[{
   SuperscriptBox["bb1", "2"], " ", 
   SuperscriptBox["x", "2"]}], "-", 
  RowBox[{"8", " ", "aa2", " ", "t", " ", "y"}], "+", 
  RowBox[{"8", " ", "bb2", " ", "t", " ", "y"}], "+", 
  RowBox[{"6", " ", "aa1", " ", "aa2", " ", "x", " ", "y"}], "+", 
  RowBox[{"6", " ", "aa2", " ", "bb1", " ", "x", " ", "y"}], "+", 
  RowBox[{"6", " ", "aa1", " ", "bb2", " ", "x", " ", "y"}], "-", 
  RowBox[{"2", " ", "bb1", " ", "bb2", " ", "x", " ", "y"}], "+", 
  RowBox[{"3", " ", 
   SuperscriptBox["aa2", "2"], " ", 
   SuperscriptBox["y", "2"]}], "+", 
  RowBox[{"6", " ", "aa2", " ", "bb2", " ", 
   SuperscriptBox["y", "2"]}], "-", 
  RowBox[{
   SuperscriptBox["bb2", "2"], " ", 
   SuperscriptBox["y", "2"]}], "+", 
  RowBox[{"36", " ", "aa1", " ", "x", " ", "z"}], "+", 
  RowBox[{"12", " ", "bb1", " ", "x", " ", "z"}], "+", 
  RowBox[{"36", " ", "aa2", " ", "y", " ", "z"}], "+", 
  RowBox[{"12", " ", "bb2", " ", "y", " ", "z"}], "-", 
  RowBox[{"8", " ", 
   SuperscriptBox["z", "2"]}]}]], "Output",
 CellChangeTimes->{
  3.909411061392*^9, 3.909411352402*^9, 3.909411409005*^9, {3.909411442207*^9,
    3.909411475868*^9}, 3.909412596533*^9, 3.9094126415*^9, {
   3.909412785499*^9, 3.909412808039*^9}, {3.909412866009*^9, 
   3.9094128929890003`*^9}, 3.909412942969*^9, 3.90947935301*^9, 
   3.916843591027*^9, 3.9168440442962503`*^9, 3.91684408723625*^9, 
   3.9168442749102497`*^9, 3.91684444940025*^9, 3.91684460523125*^9, 
   3.91684485474225*^9, 3.91684557840525*^9, 3.91684566035525*^9, 
   3.9168468004352503`*^9, 3.9168473536672497`*^9}]
}, Open  ]],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{{3.909410994042*^9, 3.9094110145150003`*^9}, {
  3.91684727439725*^9, 3.91684727993725*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.909411046785*^9, 3.909411050645*^9}, 3.909411296524*^9}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.909411284238*^9, 3.9094112876610003`*^9}, {
   3.90941133929*^9, 3.9094113772320004`*^9}, {3.909411425815*^9, 
   3.909411470738*^9}, 3.90941263392*^9, {3.909412772959*^9, 
   3.909412802539*^9}, {3.909412881559*^9, 3.909412908479*^9}, 
   3.9168472852772503`*^9}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.909412930849*^9, 3.909412936269*^9}, 
   3.91684730594725*^9}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.90947921525*^9, 3.90947925605*^9}, {3.90947933899*^9, 
   3.9094793455*^9}, 3.91684730921725*^9}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{3.91684732367725*^9}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{3.90947932198*^9, 3.91684732715725*^9}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{3.91684734747725*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"PP2", "=", 
   RowBox[{
    RowBox[{"P1", "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "\[Rule]", 
        RowBox[{
         RowBox[{"aa1", " ", "x"}], "+", " ", 
         RowBox[{"aa2", " ", "y"}]}]}], ",", 
       RowBox[{"y", "\[Rule]", " ", 
        RowBox[{
         RowBox[{"bb1", " ", "x"}], "+", " ", 
         RowBox[{"bb2", " ", "y"}]}]}], ",", 
       RowBox[{"z", "\[Rule]", 
        RowBox[{"-", "z"}]}], ",", 
       RowBox[{"t", "\[Rule]", 
        RowBox[{"-", "t"}]}]}], "}"}]}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"RP1", "=", "PP2"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"l", "=", " ", 
   RowBox[{
    RowBox[{"Table", "[", " ", 
     RowBox[{
      RowBox[{"Coefficient", "[", 
       RowBox[{
        RowBox[{"RP1", "/.", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"x", "\[Rule]", "h1"}], ",", 
           RowBox[{"y", "\[Rule]", "h2"}], ",", 
           RowBox[{"z", "\[Rule]", "h3"}], ",", 
           RowBox[{"t", "\[Rule]", "h4"}]}], "}"}]}], ",", "q", ",", "i"}], 
       "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "2", ",", "5"}], "}"}]}], "]"}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{"l", "\[Equal]", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"aa1", ",", "aa2", ",", "bb1", ",", "bb2"}], "}"}]}], "]"}], "[", 
  RowBox[{"[", "1", "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.91684398052625*^9, 3.9168440252662497`*^9}, {
  3.9168440793362503`*^9, 3.9168440798362503`*^9}, {3.91684411938625*^9, 
  3.91684412875625*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"aa1", "\[Rule]", 
    RowBox[{"-", 
     FractionBox["3", "2"]}]}], ",", 
   RowBox[{"aa2", "\[Rule]", 
    RowBox[{"-", 
     FractionBox["5", "6"]}]}], ",", 
   RowBox[{"bb1", "\[Rule]", 
    FractionBox["3", "2"]}], ",", 
   RowBox[{"bb2", "\[Rule]", 
    FractionBox["3", "2"]}]}], "}"}]], "Output",
 CellChangeTimes->{3.91684404638625*^9, 3.91684408949625*^9, 
  3.91684413890625*^9, 3.91684427688025*^9, 3.91684445199025*^9, 
  3.91684460794125*^9, 3.91684485706225*^9, 3.91684558048525*^9, 
  3.9168456619752502`*^9, 3.91684680217525*^9, 3.91684735401725*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"PP3", "=", 
   RowBox[{
    RowBox[{"P1", "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "\[Rule]", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "3"}], "/", "2"}], ")"}], " ", "x"}], "+", " ", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "5"}], "/", "6"}], ")"}], " ", "y"}]}]}], ",", 
       RowBox[{"y", "\[Rule]", " ", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"3", "/", "2"}], ")"}], " ", "x"}], "+", " ", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"3", "/", "2"}], ")"}], " ", "y"}]}]}], ",", 
       RowBox[{"z", "\[Rule]", 
        RowBox[{"-", "z"}]}], ",", 
       RowBox[{"t", "\[Rule]", 
        RowBox[{"-", "t"}]}]}], "}"}]}], "//", "Expand"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.91684419315725*^9, 3.91684424465025*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PP3", "/.", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"x", "\[Rule]", "h1"}], ",", 
    RowBox[{"y", "\[Rule]", "h2"}], ",", 
    RowBox[{"z", "\[Rule]", "h3"}], ",", 
    RowBox[{"t", "\[Rule]", "h4"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.9168442623702497`*^9, 3.91684426690025*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   RowBox[{"-", 
    FractionBox[
     RowBox[{"512", " ", 
      SuperscriptBox["q", "8"]}], "3"]}], "+", 
   FractionBox[
    RowBox[{"2560", " ", 
     SuperscriptBox["q", "9"]}], "3"], "-", 
   FractionBox[
    RowBox[{"2048", " ", 
     SuperscriptBox["q", "10"]}], "3"], "-", 
   FractionBox[
    RowBox[{"5120", " ", 
     SuperscriptBox["q", "11"]}], "3"], "+", 
   RowBox[{"2048", " ", 
    SuperscriptBox["q", "12"]}], "-", 
   FractionBox[
    RowBox[{"2560", " ", 
     SuperscriptBox["q", "13"]}], "3"], "+", 
   FractionBox[
    RowBox[{"5120", " ", 
     SuperscriptBox["q", "14"]}], "3"], "+", 
   RowBox[{"2048", " ", 
    SuperscriptBox["q", "15"]}], "-", 
   RowBox[{"2048", " ", 
    SuperscriptBox["q", "16"]}], "-", 
   FractionBox[
    RowBox[{"9728", " ", 
     SuperscriptBox["q", "17"]}], "3"], "+", 
   FractionBox[
    RowBox[{"4096", " ", 
     SuperscriptBox["q", "18"]}], "3"], "+", 
   FractionBox[
    RowBox[{"1024", " ", 
     SuperscriptBox["q", "19"]}], "3"], "-", 
   FractionBox[
    RowBox[{"16384", " ", 
     SuperscriptBox["q", "20"]}], "3"], "+", 
   FractionBox[
    RowBox[{"13312", " ", 
     SuperscriptBox["q", "21"]}], "3"], "-", 
   FractionBox[
    RowBox[{"11264", " ", 
     SuperscriptBox["q", "22"]}], "3"], "+", 
   FractionBox[
    RowBox[{"17408", " ", 
     SuperscriptBox["q", "23"]}], "3"], "-", 
   FractionBox[
    RowBox[{"1024", " ", 
     SuperscriptBox["q", "24"]}], "3"], "+", 
   FractionBox[
    RowBox[{"20480", " ", 
     SuperscriptBox["q", "25"]}], "3"], "+", 
   FractionBox[
    RowBox[{"8192", " ", 
     SuperscriptBox["q", "26"]}], "3"], "-", 
   RowBox[{"1024", " ", 
    SuperscriptBox["q", "27"]}], "-", 
   FractionBox[
    RowBox[{"10240", " ", 
     SuperscriptBox["q", "28"]}], "3"], "-", 
   FractionBox[
    RowBox[{"25600", " ", 
     SuperscriptBox["q", "29"]}], "3"], "-", 
   FractionBox[
    RowBox[{"8192", " ", 
     SuperscriptBox["q", "30"]}], "3"], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", "q", "]"}], "31"],
    SeriesData[$CellContext`q, 0, {}, 8, 31, 1],
    Editable->False]}],
  SeriesData[$CellContext`q, 0, {
    Rational[-512, 3], 
    Rational[2560, 3], 
    Rational[-2048, 3], 
    Rational[-5120, 3], 2048, 
    Rational[-2560, 3], 
    Rational[5120, 3], 2048, -2048, 
    Rational[-9728, 3], 
    Rational[4096, 3], 
    Rational[1024, 3], 
    Rational[-16384, 3], 
    Rational[13312, 3], 
    Rational[-11264, 3], 
    Rational[17408, 3], 
    Rational[-1024, 3], 
    Rational[20480, 3], 
    Rational[8192, 3], -1024, 
    Rational[-10240, 3], 
    Rational[-25600, 3], 
    Rational[-8192, 3]}, 8, 31, 1],
  Editable->False]], "Output",
 CellChangeTimes->{3.91684427692025*^9, 3.9168444520602503`*^9, 
  3.9168446080012503`*^9, 3.9168448571122503`*^9, 3.9168455805352497`*^9, 
  3.91684566205525*^9, 3.91684680226525*^9, 3.91684735407725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"PP3", "=", 
   RowBox[{
    RowBox[{"P1", "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "\[Rule]", 
        RowBox[{
         RowBox[{"aa1", " ", "x"}], "+", " ", 
         RowBox[{"aa2", " ", "y"}]}]}], ",", 
       RowBox[{"y", "\[Rule]", " ", 
        RowBox[{
         RowBox[{"bb1", " ", "x"}], "+", " ", 
         RowBox[{"bb2", " ", "y"}]}]}], ",", 
       RowBox[{"z", "\[Rule]", "z"}], ",", 
       RowBox[{"t", "\[Rule]", 
        RowBox[{"-", "t"}]}]}], "}"}]}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"RP1", "=", "PP3"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"l", "=", " ", 
   RowBox[{
    RowBox[{"Table", "[", " ", 
     RowBox[{
      RowBox[{"Coefficient", "[", 
       RowBox[{
        RowBox[{"RP1", "/.", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"x", "\[Rule]", "h1"}], ",", 
           RowBox[{"y", "\[Rule]", "h2"}], ",", 
           RowBox[{"z", "\[Rule]", "h3"}], ",", 
           RowBox[{"t", "\[Rule]", "h4"}]}], "}"}]}], ",", "q", ",", "i"}], 
       "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "2", ",", "5"}], "}"}]}], "]"}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{"l", "\[Equal]", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"aa1", ",", "aa2", ",", "bb1", ",", "bb2"}], "}"}]}], "]"}], "[", 
  RowBox[{"[", "1", "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.91684441609025*^9, 3.91684443118025*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"aa1", "\[Rule]", "0"}], ",", 
   RowBox[{"aa2", "\[Rule]", 
    RowBox[{"-", 
     FractionBox["1", "3"]}]}], ",", 
   RowBox[{"bb1", "\[Rule]", "3"}], ",", 
   RowBox[{"bb2", "\[Rule]", "2"}]}], "}"}]], "Output",
 CellChangeTimes->{3.91684445224025*^9, 3.91684460815125*^9, 
  3.91684485722225*^9, 3.91684558063525*^9, 3.91684566213525*^9, 
  3.91684680242525*^9, 3.9168473542372503`*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"PP4", "=", 
   RowBox[{
    RowBox[{"P1", "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "\[Rule]", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "/", "3"}], ")"}], " ", "y"}]}], ",", 
       RowBox[{"y", "\[Rule]", " ", 
        RowBox[{
         RowBox[{"3", " ", "x"}], "+", " ", 
         RowBox[{"2", " ", "y"}]}]}], ",", 
       RowBox[{"z", "\[Rule]", "z"}], ",", 
       RowBox[{"t", "\[Rule]", 
        RowBox[{"-", "t"}]}]}], "}"}]}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"PP4", "/.", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x", "\[Rule]", "h1"}], ",", 
      RowBox[{"y", "\[Rule]", "h2"}], ",", 
      RowBox[{"z", "\[Rule]", "h3"}], ",", 
      RowBox[{"t", "\[Rule]", "h4"}]}], "}"}]}], "//", "Expand"}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.91684451628025*^9, 3.91684457727125*^9}, {
  3.91684463934125*^9, 3.91684463991125*^9}, {3.91684476134225*^9, 
  3.91684476363225*^9}, {3.91684556123525*^9, 3.91684556642525*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "48"}], " ", 
   SuperscriptBox["t", "2"]}], "+", 
  RowBox[{"24", " ", "t", " ", "x"}], "-", 
  RowBox[{"9", " ", 
   SuperscriptBox["x", "2"]}], "+", 
  FractionBox[
   RowBox[{"56", " ", "t", " ", "y"}], "3"], "-", 
  RowBox[{"18", " ", "x", " ", "y"}], "-", 
  FractionBox[
   RowBox[{"23", " ", 
    SuperscriptBox["y", "2"]}], "3"], "+", 
  RowBox[{"36", " ", "x", " ", "z"}], "+", 
  RowBox[{"12", " ", "y", " ", "z"}], "-", 
  RowBox[{"8", " ", 
   SuperscriptBox["z", "2"]}]}]], "Input",
 NumberMarks->False],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "48"}], " ", 
   SuperscriptBox["t", "2"]}], "+", 
  RowBox[{"24", " ", "t", " ", "x"}], "-", 
  RowBox[{"9", " ", 
   SuperscriptBox["x", "2"]}], "+", 
  FractionBox[
   RowBox[{"56", " ", "t", " ", "y"}], "3"], "-", 
  RowBox[{"18", " ", "x", " ", "y"}], "-", 
  FractionBox[
   RowBox[{"23", " ", 
    SuperscriptBox["y", "2"]}], "3"], "+", 
  RowBox[{"36", " ", "x", " ", "z"}], "+", 
  RowBox[{"12", " ", "y", " ", "z"}], "-", 
  RowBox[{"8", " ", 
   SuperscriptBox["z", "2"]}]}]], "Output",
 CellChangeTimes->{3.9168447960422497`*^9, 3.91684485726225*^9, 
  3.91684558067525*^9, 3.91684566217525*^9, 3.91684680259525*^9, 
  3.91684735428725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"PP5", "=", 
   RowBox[{
    RowBox[{"P1", "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "\[Rule]", 
        RowBox[{
         RowBox[{"aa1", " ", "x"}], "+", " ", 
         RowBox[{"aa2", " ", "y"}]}]}], ",", 
       RowBox[{"y", "\[Rule]", " ", 
        RowBox[{
         RowBox[{"bb1", " ", "x"}], "+", " ", 
         RowBox[{"bb2", " ", "y"}]}]}], ",", 
       RowBox[{"z", "\[Rule]", 
        RowBox[{"-", "z"}]}], ",", 
       RowBox[{"t", "\[Rule]", "t"}]}], "}"}]}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"RP1", "=", "PP5"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"l", "=", " ", 
   RowBox[{
    RowBox[{"Table", "[", " ", 
     RowBox[{
      RowBox[{"Coefficient", "[", 
       RowBox[{
        RowBox[{"RP1", "/.", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"x", "\[Rule]", "h1"}], ",", 
           RowBox[{"y", "\[Rule]", "h2"}], ",", 
           RowBox[{"z", "\[Rule]", "h3"}], ",", 
           RowBox[{"t", "\[Rule]", "h4"}]}], "}"}]}], ",", "q", ",", "i"}], 
       "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "2", ",", "5"}], "}"}]}], "]"}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{"l", "\[Equal]", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"aa1", ",", "aa2", ",", "bb1", ",", "bb2"}], "}"}]}], "]"}], "[", 
  RowBox[{"[", "1", "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.9168456464752502`*^9, 3.9168456907852497`*^9}, 
   3.91684574701525*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"aa1", "\[Rule]", 
    RowBox[{"-", 
     FractionBox["1", "2"]}]}], ",", 
   RowBox[{"aa2", "\[Rule]", 
    RowBox[{"-", 
     FractionBox["1", "2"]}]}], ",", 
   RowBox[{"bb1", "\[Rule]", 
    RowBox[{"-", 
     FractionBox["3", "2"]}]}], ",", 
   RowBox[{"bb2", "\[Rule]", 
    FractionBox["1", "2"]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.9168456622452497`*^9, 3.91684569575525*^9}, 
   3.9168457519752502`*^9, 3.91684680273525*^9, 3.91684735442725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"PP6", "=", 
   RowBox[{
    RowBox[{"P1", "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "\[Rule]", 
        RowBox[{
         RowBox[{"aa1", " ", "x"}], "+", " ", 
         RowBox[{"aa2", " ", "y"}]}]}], ",", 
       RowBox[{"y", "\[Rule]", " ", 
        RowBox[{
         RowBox[{"bb1", " ", "x"}], "+", " ", 
         RowBox[{"bb2", " ", "y"}]}]}], ",", 
       RowBox[{"z", "\[Rule]", 
        RowBox[{"-", "z"}]}], ",", 
       RowBox[{"t", "\[Rule]", 
        RowBox[{"-", "t"}]}]}], "}"}]}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"RP1", "=", "PP6"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"l", "=", " ", 
   RowBox[{
    RowBox[{"Table", "[", " ", 
     RowBox[{
      RowBox[{"Coefficient", "[", 
       RowBox[{
        RowBox[{"RP1", "/.", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"x", "\[Rule]", "h1"}], ",", 
           RowBox[{"y", "\[Rule]", "h2"}], ",", 
           RowBox[{"z", "\[Rule]", "h3"}], ",", 
           RowBox[{"t", "\[Rule]", "h4"}]}], "}"}]}], ",", "q", ",", "i"}], 
       "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "2", ",", "5"}], "}"}]}], "]"}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{"l", "\[Equal]", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"aa1", ",", "aa2", ",", "bb1", ",", "bb2"}], "}"}]}], "]"}], "[", 
  RowBox[{"[", "1", "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.91684581598525*^9, 3.91684582834525*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"aa1", "\[Rule]", 
    RowBox[{"-", 
     FractionBox["3", "2"]}]}], ",", 
   RowBox[{"aa2", "\[Rule]", 
    RowBox[{"-", 
     FractionBox["5", "6"]}]}], ",", 
   RowBox[{"bb1", "\[Rule]", 
    FractionBox["3", "2"]}], ",", 
   RowBox[{"bb2", "\[Rule]", 
    FractionBox["3", "2"]}]}], "}"}]], "Output",
 CellChangeTimes->{3.9168458357852497`*^9, 3.91684680287525*^9, 
  3.91684735455725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"PP8", "=", 
   RowBox[{
    RowBox[{"P1", "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "\[Rule]", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "3"}], "/", "2"}], ")"}], " ", "x"}], "-", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"5", "/", "6"}], ")"}], " ", "y"}]}]}], ",", 
       RowBox[{"y", "\[Rule]", " ", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"3", "/", "2"}], ")"}], " ", "x"}], "+", " ", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"3", "/", "2"}], ")"}], " ", "y"}]}]}], ",", 
       RowBox[{"z", "\[Rule]", 
        RowBox[{"-", "z"}]}], ",", 
       RowBox[{"t", "\[Rule]", 
        RowBox[{"-", "t"}]}]}], "}"}]}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"PP8", "/.", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"x", "\[Rule]", "h1"}], ",", 
    RowBox[{"y", "\[Rule]", "h2"}], ",", 
    RowBox[{"z", "\[Rule]", "h3"}], ",", 
    RowBox[{"t", "\[Rule]", "h4"}]}], "}"}]}]}], "Input",
 CellChangeTimes->{{3.9168461119752502`*^9, 3.9168461814352503`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   RowBox[{"-", 
    FractionBox[
     RowBox[{"512", " ", 
      SuperscriptBox["q", "8"]}], "3"]}], "+", 
   FractionBox[
    RowBox[{"2560", " ", 
     SuperscriptBox["q", "9"]}], "3"], "-", 
   FractionBox[
    RowBox[{"2048", " ", 
     SuperscriptBox["q", "10"]}], "3"], "-", 
   FractionBox[
    RowBox[{"5120", " ", 
     SuperscriptBox["q", "11"]}], "3"], "+", 
   RowBox[{"2048", " ", 
    SuperscriptBox["q", "12"]}], "-", 
   FractionBox[
    RowBox[{"2560", " ", 
     SuperscriptBox["q", "13"]}], "3"], "+", 
   FractionBox[
    RowBox[{"5120", " ", 
     SuperscriptBox["q", "14"]}], "3"], "+", 
   RowBox[{"2048", " ", 
    SuperscriptBox["q", "15"]}], "-", 
   RowBox[{"2048", " ", 
    SuperscriptBox["q", "16"]}], "-", 
   FractionBox[
    RowBox[{"9728", " ", 
     SuperscriptBox["q", "17"]}], "3"], "+", 
   FractionBox[
    RowBox[{"4096", " ", 
     SuperscriptBox["q", "18"]}], "3"], "+", 
   FractionBox[
    RowBox[{"1024", " ", 
     SuperscriptBox["q", "19"]}], "3"], "-", 
   FractionBox[
    RowBox[{"16384", " ", 
     SuperscriptBox["q", "20"]}], "3"], "+", 
   FractionBox[
    RowBox[{"13312", " ", 
     SuperscriptBox["q", "21"]}], "3"], "-", 
   FractionBox[
    RowBox[{"11264", " ", 
     SuperscriptBox["q", "22"]}], "3"], "+", 
   FractionBox[
    RowBox[{"17408", " ", 
     SuperscriptBox["q", "23"]}], "3"], "-", 
   FractionBox[
    RowBox[{"1024", " ", 
     SuperscriptBox["q", "24"]}], "3"], "+", 
   FractionBox[
    RowBox[{"20480", " ", 
     SuperscriptBox["q", "25"]}], "3"], "+", 
   FractionBox[
    RowBox[{"8192", " ", 
     SuperscriptBox["q", "26"]}], "3"], "-", 
   RowBox[{"1024", " ", 
    SuperscriptBox["q", "27"]}], "-", 
   FractionBox[
    RowBox[{"10240", " ", 
     SuperscriptBox["q", "28"]}], "3"], "-", 
   FractionBox[
    RowBox[{"25600", " ", 
     SuperscriptBox["q", "29"]}], "3"], "-", 
   FractionBox[
    RowBox[{"8192", " ", 
     SuperscriptBox["q", "30"]}], "3"], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", "q", "]"}], "31"],
    SeriesData[$CellContext`q, 0, {}, 8, 31, 1],
    Editable->False]}],
  SeriesData[$CellContext`q, 0, {
    Rational[-512, 3], 
    Rational[2560, 3], 
    Rational[-2048, 3], 
    Rational[-5120, 3], 2048, 
    Rational[-2560, 3], 
    Rational[5120, 3], 2048, -2048, 
    Rational[-9728, 3], 
    Rational[4096, 3], 
    Rational[1024, 3], 
    Rational[-16384, 3], 
    Rational[13312, 3], 
    Rational[-11264, 3], 
    Rational[17408, 3], 
    Rational[-1024, 3], 
    Rational[20480, 3], 
    Rational[8192, 3], -1024, 
    Rational[-10240, 3], 
    Rational[-25600, 3], 
    Rational[-8192, 3]}, 8, 31, 1],
  Editable->False]], "Output",
 CellChangeTimes->{3.91684608780525*^9, 3.91684618727525*^9, 
  3.91684680291525*^9, 3.91684735470725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"PP9", "=", 
   RowBox[{
    RowBox[{"P1", "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "\[Rule]", 
        RowBox[{
         RowBox[{
          RowBox[{"(", "0", ")"}], " ", "x"}], "-", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"1", "/", "3"}], ")"}], " ", "y"}]}]}], ",", 
       RowBox[{"y", "\[Rule]", " ", 
        RowBox[{
         RowBox[{
          RowBox[{"(", "3", ")"}], " ", "x"}], "+", " ", 
         RowBox[{
          RowBox[{"(", "2", ")"}], " ", "y"}]}]}], ",", 
       RowBox[{"z", "\[Rule]", "z"}], ",", 
       RowBox[{"t", "\[Rule]", 
        RowBox[{"-", "t"}]}]}], "}"}]}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"PP9", "/.", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"x", "\[Rule]", "h1"}], ",", 
    RowBox[{"y", "\[Rule]", "h2"}], ",", 
    RowBox[{"z", "\[Rule]", "h3"}], ",", 
    RowBox[{"t", "\[Rule]", "h4"}]}], "}"}]}]}], "Input",
 CellChangeTimes->{{3.91684623452525*^9, 3.91684626437525*^9}, 
   3.9168463199752502`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   RowBox[{"-", 
    FractionBox[
     RowBox[{"512", " ", 
      SuperscriptBox["q", "8"]}], "3"]}], "+", 
   FractionBox[
    RowBox[{"2560", " ", 
     SuperscriptBox["q", "9"]}], "3"], "-", 
   FractionBox[
    RowBox[{"2048", " ", 
     SuperscriptBox["q", "10"]}], "3"], "-", 
   FractionBox[
    RowBox[{"5120", " ", 
     SuperscriptBox["q", "11"]}], "3"], "+", 
   RowBox[{"2048", " ", 
    SuperscriptBox["q", "12"]}], "-", 
   FractionBox[
    RowBox[{"2560", " ", 
     SuperscriptBox["q", "13"]}], "3"], "+", 
   FractionBox[
    RowBox[{"5120", " ", 
     SuperscriptBox["q", "14"]}], "3"], "+", 
   RowBox[{"2048", " ", 
    SuperscriptBox["q", "15"]}], "-", 
   RowBox[{"2048", " ", 
    SuperscriptBox["q", "16"]}], "-", 
   FractionBox[
    RowBox[{"9728", " ", 
     SuperscriptBox["q", "17"]}], "3"], "+", 
   FractionBox[
    RowBox[{"4096", " ", 
     SuperscriptBox["q", "18"]}], "3"], "+", 
   FractionBox[
    RowBox[{"1024", " ", 
     SuperscriptBox["q", "19"]}], "3"], "-", 
   FractionBox[
    RowBox[{"16384", " ", 
     SuperscriptBox["q", "20"]}], "3"], "+", 
   FractionBox[
    RowBox[{"13312", " ", 
     SuperscriptBox["q", "21"]}], "3"], "-", 
   FractionBox[
    RowBox[{"11264", " ", 
     SuperscriptBox["q", "22"]}], "3"], "+", 
   FractionBox[
    RowBox[{"17408", " ", 
     SuperscriptBox["q", "23"]}], "3"], "-", 
   FractionBox[
    RowBox[{"1024", " ", 
     SuperscriptBox["q", "24"]}], "3"], "+", 
   FractionBox[
    RowBox[{"20480", " ", 
     SuperscriptBox["q", "25"]}], "3"], "+", 
   FractionBox[
    RowBox[{"8192", " ", 
     SuperscriptBox["q", "26"]}], "3"], "-", 
   RowBox[{"1024", " ", 
    SuperscriptBox["q", "27"]}], "-", 
   FractionBox[
    RowBox[{"10240", " ", 
     SuperscriptBox["q", "28"]}], "3"], "-", 
   FractionBox[
    RowBox[{"25600", " ", 
     SuperscriptBox["q", "29"]}], "3"], "-", 
   FractionBox[
    RowBox[{"8192", " ", 
     SuperscriptBox["q", "30"]}], "3"], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", "q", "]"}], "31"],
    SeriesData[$CellContext`q, 0, {}, 8, 31, 1],
    Editable->False]}],
  SeriesData[$CellContext`q, 0, {
    Rational[-512, 3], 
    Rational[2560, 3], 
    Rational[-2048, 3], 
    Rational[-5120, 3], 2048, 
    Rational[-2560, 3], 
    Rational[5120, 3], 2048, -2048, 
    Rational[-9728, 3], 
    Rational[4096, 3], 
    Rational[1024, 3], 
    Rational[-16384, 3], 
    Rational[13312, 3], 
    Rational[-11264, 3], 
    Rational[17408, 3], 
    Rational[-1024, 3], 
    Rational[20480, 3], 
    Rational[8192, 3], -1024, 
    Rational[-10240, 3], 
    Rational[-25600, 3], 
    Rational[-8192, 3]}, 8, 31, 1],
  Editable->False]], "Output",
 CellChangeTimes->{3.9168462726852503`*^9, 3.9168463262852497`*^9, 
  3.91684680295525*^9, 3.91684735474725*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"PP10", "=", 
   RowBox[{
    RowBox[{"P1", "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "\[Rule]", 
        RowBox[{
         RowBox[{
          RowBox[{"-", 
           RowBox[{"(", 
            RowBox[{"1", "/", "2"}], ")"}]}], " ", "x"}], "-", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"1", "/", "2"}], ")"}], " ", "y"}]}]}], ",", 
       RowBox[{"y", "\[Rule]", " ", 
        RowBox[{
         RowBox[{
          RowBox[{"-", 
           RowBox[{"(", 
            RowBox[{"3", "/", "2"}], ")"}]}], " ", "x"}], "+", " ", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"1", "/", "2"}], ")"}], " ", "y"}]}]}], ",", 
       RowBox[{"z", "\[Rule]", 
        RowBox[{"-", "z"}]}], ",", 
       RowBox[{"t", "\[Rule]", "t"}]}], "}"}]}], "//", "Expand"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"PP10", "/.", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"x", "\[Rule]", "h1"}], ",", 
    RowBox[{"y", "\[Rule]", "h2"}], ",", 
    RowBox[{"z", "\[Rule]", "h3"}], ",", 
    RowBox[{"t", "\[Rule]", "h4"}]}], "}"}]}]}], "Input",
 CellChangeTimes->{{3.91684638104525*^9, 3.91684645359525*^9}}],

Cell[BoxData[
 InterpretationBox[
  SuperscriptBox[
   RowBox[{"O", "[", "q", "]"}], "31"],
  SeriesData[$CellContext`q, 0, {}, 31, 31, 1],
  Editable->False]], "Output",
 CellChangeTimes->{3.91684646005525*^9, 3.9168468029952497`*^9, 
  3.9168473547772503`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"%", "Involution", " ", "S2w32S2", " ", "correesponds", " ", "to", 
  " ", "such", " ", 
  RowBox[{"action", "."}]}]], "Input",
 CellChangeTimes->{{3.91684737282725*^9, 3.91684740289725*^9}}]
},
WindowSize->{1920, 998},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
FrontEndVersion->"11.0 for Microsoft Windows (64-bit) (July 28, 2016)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 2683, 84, 69, "Input"],
Cell[3244, 106, 208, 3, 30, "Input"],
Cell[3455, 111, 461, 13, 30, "Input"],
Cell[3919, 126, 691, 22, 50, "Input"],
Cell[4613, 150, 608, 21, 50, "Input"],
Cell[5224, 173, 330, 11, 30, "Input"],
Cell[5557, 186, 949, 30, 88, "Input"],
Cell[CellGroupData[{
Cell[6531, 220, 638, 18, 30, "Input"],
Cell[7172, 240, 1727, 34, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8936, 279, 1209, 36, 50, "Input"],
Cell[10148, 317, 2004, 45, 46, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12189, 367, 235, 8, 30, "Input"],
Cell[12427, 377, 1679, 33, 30, "Output"]
}, Open  ]],
Cell[14121, 413, 134, 3, 30, "Input"],
Cell[CellGroupData[{
Cell[14280, 420, 246, 7, 30, "Input"],
Cell[14529, 429, 1400, 24, 30, "Output"]
}, Open  ]],
Cell[15944, 456, 86, 1, 30, "Input"],
Cell[CellGroupData[{
Cell[16055, 461, 658, 23, 30, "Input"],
Cell[16716, 486, 1767, 39, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18520, 530, 726, 24, 30, "Input"],
Cell[19249, 556, 1666, 33, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[20952, 594, 683, 23, 30, "Input"],
Cell[21638, 619, 1611, 32, 30, "Output"]
}, Open  ]],
Cell[23264, 654, 112, 2, 30, "Input"],
Cell[CellGroupData[{
Cell[23401, 660, 657, 18, 30, "Input"],
Cell[24061, 680, 2483, 51, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26581, 736, 349, 12, 30, "Input"],
Cell[26933, 750, 1596, 28, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28566, 783, 386, 11, 30, "Input"],
Cell[28955, 796, 1645, 35, 30, "Output"]
}, Open  ]],
Cell[30615, 834, 457, 12, 30, "Input"],
Cell[CellGroupData[{
Cell[31097, 850, 563, 16, 30, "Input"],
Cell[31663, 868, 2236, 47, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33936, 920, 349, 12, 30, "Input"],
Cell[34288, 934, 1070, 22, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[35395, 961, 386, 11, 30, "Input"],
Cell[35784, 974, 1385, 31, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[37206, 1010, 409, 13, 30, "Input"],
Cell[37618, 1025, 1348, 25, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[39003, 1055, 1029, 31, 30, "Input"],
Cell[40035, 1088, 1570, 35, 33, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[41642, 1128, 674, 18, 30, "Input"],
Cell[42319, 1148, 1949, 43, 30, "Output"]
}, Open  ]],
Cell[44283, 1194, 157, 2, 50, "Input"],
Cell[44443, 1198, 105, 1, 30, "Input"],
Cell[44551, 1201, 313, 5, 30, "Input"],
Cell[44867, 1208, 111, 2, 30, "Input"],
Cell[44981, 1212, 146, 2, 30, "Input"],
Cell[45130, 1216, 67, 1, 30, "Input"],
Cell[45200, 1219, 85, 1, 30, "Input"],
Cell[45288, 1222, 67, 1, 30, "Input"],
Cell[CellGroupData[{
Cell[45380, 1227, 1739, 51, 88, "Input"],
Cell[47122, 1280, 618, 16, 46, "Output"]
}, Open  ]],
Cell[47755, 1299, 948, 29, 30, "Input"],
Cell[CellGroupData[{
Cell[48728, 1332, 314, 8, 30, "Input"],
Cell[49045, 1342, 2936, 97, 49, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[52018, 1444, 1612, 48, 88, "Input"],
Cell[53633, 1494, 441, 11, 46, "Output"]
}, Open  ]],
Cell[54089, 1508, 1093, 32, 50, "Input"],
Cell[CellGroupData[{
Cell[55207, 1544, 568, 18, 49, "Input"],
Cell[55778, 1564, 703, 20, 49, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[56518, 1589, 1643, 49, 88, "Input"],
Cell[58164, 1640, 513, 15, 46, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[58714, 1660, 1636, 49, 88, "Input"],
Cell[60353, 1711, 444, 14, 46, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[60834, 1730, 1174, 36, 50, "Input"],
Cell[62011, 1768, 2837, 96, 49, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[64885, 1869, 1048, 32, 50, "Input"],
Cell[65936, 1903, 2843, 96, 49, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[68816, 2004, 1174, 36, 50, "Input"],
Cell[69993, 2042, 262, 7, 30, "Output"]
}, Open  ]],
Cell[70270, 2052, 214, 4, 30, "Input"]
}
]
*)