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

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

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     27566,        765]
NotebookOptionsPosition[     26932,        739]
NotebookOutlinePosition[     27274,        754]
CellTagsIndexPosition[     27231,        751]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"ClearAll", "[", "inverseAbsoluteThirdOrder", "]"}], 
  "\[IndentingNewLine]", "\n", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"1", "/"}], "|", 
     RowBox[{"x", "+", "r"}], "|"}], "\[IndentingNewLine]", "=", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"x", "^", "2"}], "+", 
        RowBox[{"r", "^", "2"}], "+", 
        RowBox[{"2", " ", 
         RowBox[{"x", ".", "r"}]}]}], ")"}], "^", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "/", "2"}], ")"}]}], "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"1", "/"}], "|", "r", "|"}], ")"}], 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "+", 
          RowBox[{
           RowBox[{"x", "^", "2"}], "/", 
           RowBox[{"r", "^", "2"}]}], "+", 
          RowBox[{"2", " ", 
           RowBox[{"x", ".", "r"}]}]}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "/", "2"}], ")"}]}], "\\", "approx"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"1", "/"}], "|", "r", "|"}], ")"}], 
      RowBox[{"(", 
       RowBox[{"1", "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"x", "^", "2"}], "/", 
            RowBox[{"r", "^", "2"}]}], "+", 
           RowBox[{"2", " ", 
            RowBox[{"x", ".", "r"}]}]}], ")"}], "/", "2"}], "+", 
        RowBox[{"3", " ", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"x", "^", "2"}], "/", 
              RowBox[{"r", "^", "2"}]}], "+", 
             RowBox[{"2", " ", 
              RowBox[{"x", ".", "r"}]}]}], ")"}], "^", "2"}], "/", "8"}]}]}], 
       ")"}]}]}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"inverseAbsoluteThirdOrder", "[", 
     RowBox[{"x_", ",", "r_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "/", 
       RowBox[{"Norm", "[", "r", "]"}]}], ")"}], 
     RowBox[{"(", 
      RowBox[{"1", "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"x", ".", "x"}], ")"}], "/", 
           RowBox[{"(", 
            RowBox[{"r", ".", "r"}], ")"}]}], "+", 
          RowBox[{"2", " ", 
           RowBox[{
            RowBox[{"x", ".", "r"}], "/", 
            RowBox[{"(", 
             RowBox[{"r", ".", "r"}], ")"}]}]}]}], ")"}], "/", "2"}], "+", 
       RowBox[{"3", " ", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"x", ".", "x"}], ")"}], "/", 
             RowBox[{"(", 
              RowBox[{"r", ".", "r"}], ")"}]}], "+", 
            RowBox[{"2", " ", 
             RowBox[{
              RowBox[{"x", ".", "r"}], "/", 
              RowBox[{"(", 
               RowBox[{"r", ".", "r"}], ")"}]}]}]}], ")"}], "^", "2"}], "/", 
         "8"}]}]}], ")"}]}]}], " ", ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"x1", "=", 
   RowBox[{"{", 
    RowBox[{
    "\"\<\!\(\*SuperscriptBox[SubscriptBox[\(x\), \(1\)], \((1)\)]\)\>\"", 
     ",", "\"\<\!\(\*SuperscriptBox[SubscriptBox[\(x\), \(2\)], \
\((1)\)]\)\>\"", ",", 
     "\"\<\!\(\*SuperscriptBox[SubscriptBox[\(x\), \(3\)], \((1)\)]\)\>\""}], 
    "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"x2", "=", 
   RowBox[{"{", 
    RowBox[{
    "\"\<\!\(\*SuperscriptBox[SubscriptBox[\(x\), \(1\)], \((2)\)]\)\>\"", 
     ",", "\"\<\!\(\*SuperscriptBox[SubscriptBox[\(x\), \(2\)], \
\((2)\)]\)\>\"", ",", 
     "\"\<\!\(\*SuperscriptBox[SubscriptBox[\(x\), \(3\)], \((2)\)]\)\>\""}], 
    "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"r", "=", 
    RowBox[{"{", 
     RowBox[{
     "\"\<\!\(\*SubscriptBox[\(R\), \(1\)]\)\>\"", ",", 
      "\"\<\!\(\*SubscriptBox[\(R\), \(2\)]\)\>\"", ",", 
      "\"\<\!\(\*SubscriptBox[\(R\), \(3\)]\)\>\""}], "}"}]}], ";"}], 
  "\n"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"taylor", " ", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"1", "/", 
        RowBox[{"Norm", "[", "r", "]"}]}], "+", 
       RowBox[{"inverseAbsoluteThirdOrder", "[", 
        RowBox[{
         RowBox[{"x1", "-", "x2"}], ",", "r"}], "]"}], "-", 
       RowBox[{"inverseAbsoluteThirdOrder", "[", 
        RowBox[{"x1", ",", "r"}], "]"}], "-", 
       RowBox[{"inverseAbsoluteThirdOrder", "[", 
        RowBox[{"x2", ",", 
         RowBox[{"-", "r"}]}], "]"}]}], ")"}], ")"}], " ", "//", " ", 
    "FullSimplify"}]}], "\n"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.658579430875585*^9, 3.658579456042781*^9}, {
  3.658579495014436*^9, 3.658579507832864*^9}, {3.658579541786708*^9, 
  3.658579573625909*^9}, {3.658579632756649*^9, 3.6585798026817517`*^9}, {
  3.658579835851212*^9, 3.658579838100144*^9}, {3.6585798964111605`*^9, 
  3.658579901317869*^9}, {3.6585799606704836`*^9, 3.658579977665346*^9}, {
  3.6586108061725836`*^9, 3.658610807721204*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"4", " ", 
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "2"], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "2"], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "2"]}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
       "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((1)\\)]\\)\"\>"}], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\
\\)], \\((1)\\)]\\)\"\>", "2"], "+", 
       RowBox[{
       "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((1)\\)]\\)\"\>"}], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\
\\)], \\((1)\\)]\\)\"\>", "2"], "+", 
       RowBox[{
       "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((1)\\)]\\)\"\>"}], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\
\\)], \\((1)\\)]\\)\"\>", "2"]}], ")"}]}], "-", 
    RowBox[{"3", " ", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{
        RowBox[{
        "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
         " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((1)\\)]\\)\"\>"}], "+", 
        SuperscriptBox["\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \
\\(1\\)], \\((1)\\)]\\)\"\>", "2"], "+", 
        RowBox[{
        "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
         " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((1)\\)]\\)\"\>"}], "+", 
        SuperscriptBox["\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \
\\(2\\)], \\((1)\\)]\\)\"\>", "2"], "+", 
        RowBox[{
        "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
         " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((1)\\)]\\)\"\>"}], "+", 
        SuperscriptBox["\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \
\\(3\\)], \\((1)\\)]\\)\"\>", "2"]}], ")"}], "2"]}], "-", 
    RowBox[{"4", " ", 
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "2"], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "2"], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "2"]}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((1)\\)]\\)\"\>", 
          "-", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((2)\\)]\\)\"\>"}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
          "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>"}], 
          "+", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((1)\\)]\\)\"\>", 
          "-", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((2)\\)]\\)\"\>"}], ")"}]}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((1)\\)]\\)\"\>", 
          "-", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((2)\\)]\\)\"\>"}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
          "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>"}], 
          "+", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((1)\\)]\\)\"\>", 
          "-", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((2)\\)]\\)\"\>"}], ")"}]}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((1)\\)]\\)\"\>", 
          "-", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((2)\\)]\\)\"\>"}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
          "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>"}], 
          "+", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((1)\\)]\\)\"\>", 
          "-", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((2)\\)]\\)\"\>"}], ")"}]}]}], ")"}]}], "+", 
    RowBox[{"3", " ", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          
          RowBox[{"\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \
\\(1\\)], \\((1)\\)]\\)\"\>", 
           "-", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((2)\\)]\\)\"\>"}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
           "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>"}], 
           "+", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((1)\\)]\\)\"\>", 
           "-", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((2)\\)]\\)\"\>"}], ")"}]}], "+", 
        RowBox[{
         RowBox[{"(", 
          
          RowBox[{"\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \
\\(2\\)], \\((1)\\)]\\)\"\>", 
           "-", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((2)\\)]\\)\"\>"}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
           "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>"}], 
           "+", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((1)\\)]\\)\"\>", 
           "-", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((2)\\)]\\)\"\>"}], ")"}]}], "+", 
        RowBox[{
         RowBox[{"(", 
          
          RowBox[{"\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \
\\(3\\)], \\((1)\\)]\\)\"\>", 
           "-", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((2)\\)]\\)\"\>"}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
           "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>"}], 
           "+", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((1)\\)]\\)\"\>", 
           "-", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((2)\\)]\\)\"\>"}], ")"}]}]}], ")"}], "2"]}], "+", 
    RowBox[{"4", " ", 
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "2"], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "2"], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "2"]}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "2"}], 
        " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((2)\\)]\\)\"\>"}], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\
\\)], \\((2)\\)]\\)\"\>", "2"], "-", 
       RowBox[{
       "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((2)\\)]\\)\"\>"}], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\
\\)], \\((2)\\)]\\)\"\>", "2"], "-", 
       RowBox[{
       "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((2)\\)]\\)\"\>"}], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\
\\)], \\((2)\\)]\\)\"\>", "2"]}], ")"}]}], "-", 
    RowBox[{"3", " ", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "2"}], 
         " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
         " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((2)\\)]\\)\"\>"}], "+", 
        SuperscriptBox["\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \
\\(1\\)], \\((2)\\)]\\)\"\>", "2"], "-", 
        RowBox[{
        "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
         " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((2)\\)]\\)\"\>"}], "+", 
        SuperscriptBox["\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \
\\(2\\)], \\((2)\\)]\\)\"\>", "2"], "-", 
        RowBox[{
        "2", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
         " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((2)\\)]\\)\"\>"}], "+", 
        SuperscriptBox["\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \
\\(3\\)], \\((2)\\)]\\)\"\>", "2"]}], ")"}], "2"]}]}], ")"}], "/", 
  RowBox[{"(", 
   RowBox[{"8", " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "2"], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "2"], "+", 
       SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "2"]}], ")"}], "2"], " ", 
    SqrtBox[
     RowBox[{
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "]"}], "2"]}]]}], ")"}]}]], "Output",
 CellChangeTimes->{{3.658579432844396*^9, 3.6585794566209235`*^9}, 
   3.6585795091154547`*^9, {3.6585795545610914`*^9, 3.658579574297849*^9}, 
   3.658579695226824*^9, 3.658579775469773*^9, 3.6585798089744344`*^9, 
   3.658579892180535*^9, 3.6585799499680505`*^9, 3.6585799999341197`*^9, 
   3.658610817794407*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.658579519584515*^9, 3.6585795215863485`*^9}, 
   3.658579907623641*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"explicit", " ", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"x1", ".", "x2"}], "/", 
     RowBox[{
      RowBox[{"Norm", "[", "r", "]"}], "^", "3"}]}], "-", 
    RowBox[{"3", " ", 
     RowBox[{"(", 
      RowBox[{"x1", ".", "r"}], ")"}], 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"x2", ".", "r"}], ")"}], "/", 
      RowBox[{
       RowBox[{"Norm", "[", "r", "]"}], "^", "5"}]}]}]}], " ", "//", " ", 
   "Expand"}]}], "\n", 
 RowBox[{" "}]}], "Input",
 CellChangeTimes->{{3.658611988046239*^9, 3.658612009174883*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   FractionBox[
    RowBox[{"3", " ", 
     SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
      "2"], " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((1)\\)]\\)\"\>", 
     " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \\((2)\
\\)]\\)\"\>"}], 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox[
        RowBox[{
        "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
         "]"}], "2"], "+", 
       SuperscriptBox[
        RowBox[{
        "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
         "]"}], "2"], "+", 
       SuperscriptBox[
        RowBox[{
        "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
         "]"}], "2"]}], ")"}], 
     RowBox[{"5", "/", "2"}]]]}], "-", 
  FractionBox[
   RowBox[{
   "3", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((2)\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((1)\\)]\\)\"\>"}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "]"}], "2"]}], ")"}], 
    RowBox[{"5", "/", "2"}]]], "-", 
  FractionBox[
   RowBox[{
   "3", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((1)\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((2)\\)]\\)\"\>"}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "]"}], "2"]}], ")"}], 
    RowBox[{"5", "/", "2"}]]], "-", 
  FractionBox[
   RowBox[{"3", " ", 
    SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", "2"],
     " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \\((1)\
\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((2)\\)]\\)\"\>"}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "]"}], "2"]}], ")"}], 
    RowBox[{"5", "/", "2"}]]], "-", 
  FractionBox[
   RowBox[{
   "3", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((2)\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((1)\\)]\\)\"\>"}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "]"}], "2"]}], ")"}], 
    RowBox[{"5", "/", "2"}]]], "-", 
  FractionBox[
   RowBox[{
   "3", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((2)\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((1)\\)]\\)\"\>"}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "]"}], "2"]}], ")"}], 
    RowBox[{"5", "/", "2"}]]], "-", 
  FractionBox[
   RowBox[{
   "3", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((1)\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((2)\\)]\\)\"\>"}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "]"}], "2"]}], ")"}], 
    RowBox[{"5", "/", "2"}]]], "-", 
  FractionBox[
   RowBox[{
   "3", " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((1)\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((2)\\)]\\)\"\>"}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "]"}], "2"]}], ")"}], 
    RowBox[{"5", "/", "2"}]]], "-", 
  FractionBox[
   RowBox[{"3", " ", 
    SuperscriptBox["\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", "2"],
     " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \\((1)\
\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((2)\\)]\\)\"\>"}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "]"}], "2"]}], ")"}], 
    RowBox[{"5", "/", "2"}]]], "+", 
  FractionBox[
   RowBox[{"\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((1)\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(1\\)], \
\\((2)\\)]\\)\"\>"}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "]"}], "2"]}], ")"}], 
    RowBox[{"3", "/", "2"}]]], "+", 
  FractionBox[
   RowBox[{"\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((1)\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(2\\)], \
\\((2)\\)]\\)\"\>"}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "]"}], "2"]}], ")"}], 
    RowBox[{"3", "/", "2"}]]], "+", 
  FractionBox[
   RowBox[{"\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((1)\\)]\\)\"\>", 
    " ", "\<\"\\!\\(\\*SuperscriptBox[SubscriptBox[\\(x\\), \\(3\\)], \
\\((2)\\)]\\)\"\>"}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(1\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(2\\)]\\)\"\>", 
        "]"}], "2"], "+", 
      SuperscriptBox[
       RowBox[{
       "Abs", "[", "\<\"\\!\\(\\*SubscriptBox[\\(R\\), \\(3\\)]\\)\"\>", 
        "]"}], "2"]}], ")"}], 
    RowBox[{"3", "/", "2"}]]]}]], "Output",
 CellChangeTimes->{{3.6586119937489047`*^9, 3.658612009799945*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"taylor", " ", "-", " ", "explicit"}], " ", "//", " ", 
  "FullSimplify"}]], "Input",
 CellChangeTimes->{{3.6586120495031414`*^9, 3.6586120564620595`*^9}}],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{3.658612657152606*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.658579534889491*^9, 3.6585795349519415`*^9}}]
},
WindowSize->{707, 636},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 5151, 150, 554, "Input"],
Cell[5733, 174, 10192, 239, 417, "Output"]
}, Open  ]],
Cell[15940, 416, 120, 2, 39, "Input"],
Cell[CellGroupData[{
Cell[16085, 422, 567, 18, 91, "Input"],
Cell[16655, 442, 9858, 279, 754, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26550, 726, 189, 4, 39, "Input"],
Cell[26742, 732, 77, 1, 40, "Output"]
}, Open  ]],
Cell[26834, 736, 94, 1, 39, InheritFromParent]
}
]
*)

(* End of internal cache information *)