(* Content-type: application/mathematica *)

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

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     28383,        801]
NotebookOptionsPosition[     26048,        723]
NotebookOutlinePosition[     27056,        758]
CellTagsIndexPosition[     26930,        752]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"Example", " ", "file", " ", "for", " ", 
      RowBox[{"GREATER2", ".", "m"}]}], ";", " ", 
     RowBox[{"run", " ", "this", " ", "cell", " ", "first"}]}], ",", " ", 
    RowBox[{"then", " ", "the", " ", 
     RowBox[{"rest", ".", " ", "If"}], " ", "this", " ", 
     RowBox[{"doesn", "'"}], "t", " ", "work", " ", "then", " ", "you", " ", 
     "put", " ", 
     RowBox[{"the", " ", ".", "m"}], " ", "files", " ", "in", " ", "the", " ",
      "wrong", " ", 
     RowBox[{
     "directory", ".", " ", "\[IndentingNewLine]", "\[IndentingNewLine]", 
      "Read"}], " ", "the", " ", 
     RowBox[{"GREATER2", ".", "m"}], " ", "file", " ", "for", " ", "more", 
     " ", 
     RowBox[{
     "commands", ".", "\[IndentingNewLine]", "\[IndentingNewLine]", "If"}], 
     " ", "you", " ", "cannot", " ", "put", " ", "things", " ", "in", " ", 
     "the", " ", "Mathematica", " ", "apps", " ", "directory", " ", "for", 
     " ", "some", " ", "reason"}], ",", " ", 
    RowBox[{
    "uncomment", " ", "the", " ", "next", " ", "line", " ", "and", " ", 
     "point", " ", "it", " ", "to", " ", "the", " ", "location", " ", "of", 
     " ", "your", " ", 
     RowBox[{"GREATER2", ".", "m"}], " ", 
     RowBox[{"file", "."}]}]}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"PrependTo", "[", 
     RowBox[{"$Path", ",", " ", 
      RowBox[{"ToFileName", "[", 
       RowBox[{"{", 
        RowBox[{
        "$HomeDirectory", ",", " ", "\"\<physics\>\"", ",", " ", 
         "\"\<my-math-apps\>\""}], "}"}], "]"}]}], "]"}], ";"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Needs", "[", "\"\<GREATER2`\>\"", "]"}], ";"}]}]], "Input",
 CellChangeTimes->{{3.43523344003125*^9, 3.43523351621875*^9}, {
  3.435233935625*^9, 3.43523394528125*^9}, {3.549981558480381*^9, 
  3.5499815670208387`*^9}, {3.549981764036931*^9, 3.549981809920496*^9}, {
  3.54998187744914*^9, 3.5499818934180737`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"define", " ", "our", " ", "coordinates"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"X", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{"t", ",", "r", ",", "\[Phi]"}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "define", " ", "a", " ", "line", " ", "element", " ", "for", " ", 
     "AdS3"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ds2", " ", "=", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"-", 
        RowBox[{"r", "^", "2"}]}], " ", 
       RowBox[{"dt", "^", "2"}]}], " ", "+", " ", 
      RowBox[{
       RowBox[{"dr", "^", "2"}], "/", 
       RowBox[{"r", "^", "2"}]}], " ", "+", " ", 
      RowBox[{
       RowBox[{"r", "^", "2"}], " ", 
       RowBox[{"d\[Phi]", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"define", " ", "the", " ", "metric"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{"Note", ":", " ", 
      RowBox[{
      "I", " ", "like", " ", "to", " ", "use", " ", "the", " ", "tensor", " ",
        "name", " ", "to", " ", "keep", " ", "track", " ", "of", " ", "where",
        " ", "the", " ", "indices", " ", 
       RowBox[{"are", ".", "  ", "Since"}], " ", "the", " ", "metric", " ", 
       "has", " ", "two", " ", "down", " ", "indices"}]}], ",", " ", 
     RowBox[{"I", " ", "end", " ", "the", " ", "name", " ", "with", " ", 
      RowBox[{"\"\<dd\>\"", ".", "  ", "If"}], " ", "I", " ", "defined", " ", 
      "a", " ", "tensor", " ", "with", " ", "two", " ", "up", " ", "indices", 
      " ", "and", " ", "two", " ", "down"}], ",", " ", 
     RowBox[{
      RowBox[{"I", "'"}], "d", " ", "call", " ", "it", " ", "Xuudd"}]}], " ", 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Gdd", " ", "=", " ", 
     RowBox[{"Metric", "[", 
      RowBox[{"ds2", ",", " ", "X"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"inverse", " ", "metric"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Guu", " ", "=", " ", 
     RowBox[{"IMetric", "[", "Gdd", "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"Gdd", "//", "SMF"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.43523344003125*^9, 3.43523366875*^9}, {
  3.435233963765625*^9, 3.4352340548125*^9}, {3.549981596432922*^9, 
  3.549981644711142*^9}, {3.549981830313285*^9, 3.549981831695394*^9}, {
  3.549982135069282*^9, 3.549982139076964*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", 
       SuperscriptBox["r", "2"]}], "0", "0"},
     {"0", 
      FractionBox["1", 
       SuperscriptBox["r", "2"]], "0"},
     {"0", "0", 
      SuperscriptBox["r", "2"]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{
  3.549981961737955*^9, {3.549982140563778*^9, 3.549982143671213*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Compute", " ", "the", " ", "Ricci", " ", "scalar"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"R", " ", "=", " ", 
    RowBox[{"SCurvature", "[", 
     RowBox[{"Gdd", ",", " ", "X"}], "]"}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"Define", " ", "a", " ", "vector", " ", "field"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Zeta]u", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"e", "[", 
        RowBox[{"t", "-", "\[Phi]"}], "]"}], ",", " ", "0", ",", "0"}], 
      "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"Norm", " ", "of", " ", "\[Zeta]"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Zeta]u", ".", "Gdd", ".", "\[Zeta]u"}], "//", "S"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "Take", " ", "lie", " ", "derivative", " ", "of", " ", "G", " ", "along",
       " ", "\[Zeta]"}], ";", " ", 
     RowBox[{"the", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "1"}], ",", 
        RowBox[{"-", "1"}]}], "}"}], " ", "indicates", " ", "we", " ", "are", 
      " ", "taking", " ", "the", " ", "Lie", " ", "derivative", " ", "of", 
      " ", "a", " ", "tensor", " ", "with", " ", "two", " ", "lower", " ", 
      "indices"}]}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"dGdd", " ", "=", " ", 
     RowBox[{"LieD", "[", 
      RowBox[{"\[Zeta]u", ",", " ", "Gdd", ",", " ", "X", ",", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "1"}], ",", 
         RowBox[{"-", "1"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "Take", " ", "the", " ", "covariant", " ", "derivative", " ", "of", " ", 
      "\[Zeta]"}], ";", " ", 
     RowBox[{"the", " ", 
      RowBox[{"{", "1", "}"}], " ", "indicates", " ", "the", " ", "\[Zeta]", 
      " ", "has", " ", "one", " ", "upper", " ", "index"}]}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Dd\[Zeta]u", " ", "=", " ", 
     RowBox[{"CoD", "[", 
      RowBox[{"\[Zeta]u", ",", " ", "Gdd", ",", " ", "X", ",", " ", 
       RowBox[{"{", "1", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"Display", " ", "results", " ", "nicely"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{"dGdd", "//", "SMF"}], "\[IndentingNewLine]", 
   RowBox[{"Dd\[Zeta]u", "//", "SMF"}], "\[IndentingNewLine]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"-", "6"}]], "Output",
 CellChangeTimes->{3.549981826316345*^9, 3.5499819617842827`*^9, 
  3.5499821437094727`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   SuperscriptBox["r", "2"]}], " ", 
  SuperscriptBox[
   RowBox[{"e", "[", 
    RowBox[{"t", "-", "\[Phi]"}], "]"}], "2"]}]], "Output",
 CellChangeTimes->{3.549981826316345*^9, 3.5499819617842827`*^9, 
  3.5499821437130527`*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{
       RowBox[{"-", "2"}], " ", 
       SuperscriptBox["r", "2"], " ", 
       RowBox[{
        SuperscriptBox["e", "\[Prime]",
         MultilineFunction->None], "[", 
        RowBox[{"t", "-", "\[Phi]"}], "]"}]}], "0", 
      RowBox[{
       SuperscriptBox["r", "2"], " ", 
       RowBox[{
        SuperscriptBox["e", "\[Prime]",
         MultilineFunction->None], "[", 
        RowBox[{"t", "-", "\[Phi]"}], "]"}]}]},
     {"0", "0", "0"},
     {
      RowBox[{
       SuperscriptBox["r", "2"], " ", 
       RowBox[{
        SuperscriptBox["e", "\[Prime]",
         MultilineFunction->None], "[", 
        RowBox[{"t", "-", "\[Phi]"}], "]"}]}], "0", "0"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.549981826316345*^9, 3.5499819617842827`*^9, 
  3.5499821437153597`*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{
       SuperscriptBox["e", "\[Prime]",
        MultilineFunction->None], "[", 
       RowBox[{"t", "-", "\[Phi]"}], "]"}], 
      RowBox[{
       SuperscriptBox["r", "3"], " ", 
       RowBox[{"e", "[", 
        RowBox[{"t", "-", "\[Phi]"}], "]"}]}], "0"},
     {
      FractionBox[
       RowBox[{"e", "[", 
        RowBox[{"t", "-", "\[Phi]"}], "]"}], "r"], "0", "0"},
     {
      RowBox[{"-", 
       RowBox[{
        SuperscriptBox["e", "\[Prime]",
         MultilineFunction->None], "[", 
        RowBox[{"t", "-", "\[Phi]"}], "]"}]}], "0", "0"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.549981826316345*^9, 3.5499819617842827`*^9, 
  3.549982143720375*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"Change", " ", "coordinate", " ", "systems"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"X2", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{"xplus", ",", " ", "xminus", ",", " ", "y"}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"trans12", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"xplus", "\[Equal]", " ", 
        RowBox[{"t", "+", "\[Phi]"}]}], ",", " ", 
       RowBox[{"xminus", " ", "\[Equal]", " ", 
        RowBox[{"t", " ", "-", " ", "\[Phi]"}]}], ",", " ", 
       RowBox[{"y", " ", "\[Equal]", " ", 
        RowBox[{"1", "/", "r"}]}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "last", " ", "parameter", " ", "specifies", " ", "the", " ", "index", " ",
      "positions", " ", "on", " ", "Gdd"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Gdd2", " ", "=", " ", 
     RowBox[{"ChangeCoords", "[", 
      RowBox[{
      "Gdd", ",", " ", "X", ",", " ", "X2", ",", " ", "trans12", ",", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "1"}], ",", 
         RowBox[{"-", "1"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Zeta]u2", " ", "=", " ", 
     RowBox[{"ChangeCoords", "[", 
      RowBox[{
      "\[Zeta]u", ",", " ", "X", ",", " ", "X2", ",", " ", "trans12", ",", 
       " ", 
       RowBox[{"{", "1", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"Gdd2", "//", "SMF"}], "\[IndentingNewLine]", 
   "\[Zeta]u2"}]}]], "Input",
 CellChangeTimes->{{3.549981648688156*^9, 3.549981738442952*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", 
      RowBox[{"-", 
       FractionBox["1", 
        RowBox[{"2", " ", 
         SuperscriptBox["y", "2"]}]]}], "0"},
     {
      RowBox[{"-", 
       FractionBox["1", 
        RowBox[{"2", " ", 
         SuperscriptBox["y", "2"]}]]}], "0", "0"},
     {"0", "0", 
      FractionBox["1", 
       SuperscriptBox["y", "2"]]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.54998169108602*^9, 3.549981739722728*^9}, 
   3.5499818263657913`*^9, 3.549981961867343*^9, 3.549982143776432*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"e", "[", "xminus", "]"}], ",", 
   RowBox[{"e", "[", "xminus", "]"}], ",", "0"}], "}"}]], "Output",
 CellChangeTimes->{{3.54998169108602*^9, 3.549981739722728*^9}, 
   3.5499818263657913`*^9, 3.549981961867343*^9, 3.549982143779293*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"Compute", " ", "some", " ", "tensor", " ", "expressions"}], " ", 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Vu", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{"v1", ",", "v2", ",", "v3"}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Kdd", " ", "=", " ", 
    RowBox[{
     RowBox[{"ExtrinsicCurvature", "[", 
      RowBox[{"Gdd", ",", " ", "2", ",", " ", "1", ",", " ", "X"}], "]"}], "//", 
     RowBox[{
      RowBox[{"Simplify", "[", 
       RowBox[{"#", ",", " ", 
        RowBox[{"r", ">", "0"}]}], "]"}], "&"}]}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "TensorMath", " ", "computes", " ", "arbitrary", " ", "tensor", " ", 
     "expressions", " ", "written", " ", "in", " ", "Einstein", " ", 
     "summation", " ", "convention"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{"sometensor", " ", "=", " ", 
    RowBox[{"TensorMath", "[", 
     RowBox[{
      RowBox[{"Vu", "[", 
       RowBox[{"[", "\[Mu]", "]"}], "]"}], 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"Kdd", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"-", "\[Mu]"}], ",", 
           RowBox[{"-", "\[Nu]"}]}], "]"}], "]"}], " ", "+", " ", 
        RowBox[{"3", " ", 
         RowBox[{"f", "[", "x", "]"}], " ", 
         RowBox[{"Gdd", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"-", "\[Mu]"}], ",", 
            RowBox[{"-", "\[Nu]"}]}], "]"}], "]"}]}]}], ")"}]}], " ", "]"}]}],
    "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{
     "The", " ", "next", " ", "statement", " ", "is", " ", "equivalent"}], 
     ",", " ", 
     RowBox[{
     "but", " ", "writing", " ", "things", " ", "this", " ", "way", " ", 
      "gets", " ", "tricky", " ", "when", " ", "tensors", " ", "have", " ", 
      "more", " ", "than", " ", "2", " ", "indices"}]}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{"sometensor2", " ", "=", " ", 
    RowBox[{
     RowBox[{"Vu", ".", "Kdd"}], " ", "+", " ", 
     RowBox[{"3", " ", 
      RowBox[{"f", "[", "x", "]"}], " ", 
      RowBox[{"Vu", ".", "Gdd"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sometensor", " ", "-", " ", "sometensor2"}], "//", 
    "S"}]}]}]], "Input",
 CellChangeTimes->{{3.549981896959918*^9, 3.5499819387335033`*^9}, {
  3.549981969358603*^9, 3.549982085987924*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     SuperscriptBox["r", "2"], ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", 
     RowBox[{"-", 
      SuperscriptBox["r", "2"]}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.5499819408710546`*^9, 3.549981973016027*^9}, {
   3.549982015579833*^9, 3.5499820891622553`*^9}, 3.549982143822054*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{
     SuperscriptBox["r", "2"], " ", "v1"}], "-", 
    RowBox[{"3", " ", 
     SuperscriptBox["r", "2"], " ", "v1", " ", 
     RowBox[{"f", "[", "x", "]"}]}]}], ",", 
   FractionBox[
    RowBox[{"3", " ", "v2", " ", 
     RowBox[{"f", "[", "x", "]"}]}], 
    SuperscriptBox["r", "2"]], ",", 
   RowBox[{
    RowBox[{
     RowBox[{"-", 
      SuperscriptBox["r", "2"]}], " ", "v3"}], "+", 
    RowBox[{"3", " ", 
     SuperscriptBox["r", "2"], " ", "v3", " ", 
     RowBox[{"f", "[", "x", "]"}]}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.5499819408710546`*^9, 3.549981973016027*^9}, {
   3.549982015579833*^9, 3.5499820891622553`*^9}, 3.549982143824443*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{
     SuperscriptBox["r", "2"], " ", "v1"}], "-", 
    RowBox[{"3", " ", 
     SuperscriptBox["r", "2"], " ", "v1", " ", 
     RowBox[{"f", "[", "x", "]"}]}]}], ",", 
   FractionBox[
    RowBox[{"3", " ", "v2", " ", 
     RowBox[{"f", "[", "x", "]"}]}], 
    SuperscriptBox["r", "2"]], ",", 
   RowBox[{
    RowBox[{
     RowBox[{"-", 
      SuperscriptBox["r", "2"]}], " ", "v3"}], "+", 
    RowBox[{"3", " ", 
     SuperscriptBox["r", "2"], " ", "v3", " ", 
     RowBox[{"f", "[", "x", "]"}]}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.5499819408710546`*^9, 3.549981973016027*^9}, {
   3.549982015579833*^9, 3.5499820891622553`*^9}, 3.5499821438266687`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0", ",", "0", ",", "0"}], "}"}]], "Output",
 CellChangeTimes->{{3.5499819408710546`*^9, 3.549981973016027*^9}, {
   3.549982015579833*^9, 3.5499820891622553`*^9}, 3.549982143828452*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"Use", " ", "?", "functionname"}], " ", "to", " ", "see", " ", 
    "the", " ", "definition", " ", "of", " ", "a", " ", "function", " ", "so",
     " ", "you", " ", "can", " ", "see", " ", "the", " ", "required", " ", 
    "syntax"}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"?", "Metric"}], "\[IndentingNewLine]", 
   RowBox[{"?", "Riemann"}], "\[IndentingNewLine]", 
   RowBox[{"?", "ExtrinsicCurvature"}]}]}]], "Input",
 CellChangeTimes->{{3.549981928993081*^9, 3.549981931941813*^9}, {
  3.54998209904242*^9, 3.549982187122085*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 StyleBox["\<\"\!\(\\\"Metric[\\\\!\\\\(\\\\*\\\\nStyleBox[\\\\\\\"ds2\\\\\\\"\
,\\\\nFontSlant->\\\\\\\"Italic\\\\\\\"]\\\\), \
{\\\\!\\\\(\\\\*SubscriptBox[\\\\nStyleBox[\\\\\\\"x\\\\\\\",\\\\nFontSlant->\
\\\\\\\"Italic\\\\\\\"], \\\\(1\\\\)]\\\\), \
\\\\!\\\\(\\\\*SubscriptBox[\\\\nStyleBox[\\\\\\\"x\\\\\\\",\\\\nFontSlant->\\\
\\\\\"Italic\\\\\\\"], \\\\(2\\\\)]\\\\), ...}] returns the metric tensor \
specified by the line element \
\\\\!\\\\(\\\\*\\\\nStyleBox[\\\\\\\"ds2\\\\\\\",\\\\nFontSlant->\\\\\\\"\
Italic\\\\\\\"]\\\\) with coordinates \
\\\\!\\\\(\\\\*SubscriptBox[\\\\nStyleBox[\\\\\\\"x\\\\\\\",\\\\nFontSlant->\\\
\\\\\"Italic\\\\\\\"], \\\\(i\\\\)]\\\\). (Just like ds2met but does not \
require you to define differentials.)\\\"\)\"\>", "MSG"]], "Print", \
"PrintUsage",
 CellChangeTimes->{3.549982187571146*^9},
 CellTags->"Info3549967787-9563680"],

Cell[BoxData[
 StyleBox["\<\"Riemann[g,x], with g a n.n-matrix and x n-vector of\\n  \
coordinates, gives the Riemann tensor (1st upper, three lower\\n  indices).\"\
\>", "MSG"]], "Print", "PrintUsage",
 CellChangeTimes->{3.549982187820879*^9},
 CellTags->"Info3549967787-9563680"],

Cell["GREATER2`ExtrinsicCurvature", "Print", "PrintUsage",
 CellChangeTimes->{3.549982188055702*^9},
 CellTags->"Info3549967787-9563680"],

Cell[BoxData[
 InterpretationBox[GridBox[{
    {GridBox[{
       {
        RowBox[{
         RowBox[{"ExtrinsicCurvature", "[", 
          RowBox[{"Gdd_", ",", "coord_", ",", "csign_", ",", "X_"}], "]"}], ":=", 
         RowBox[{"Module", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
            "nd", ",", "\[Gamma]dd", ",", "delN", ",", "delSym", ",", "Kdd", 
             ",", "\[Gamma]ud"}], "}"}], ",", 
           RowBox[{
            RowBox[{"nd", "=", 
             RowBox[{"UnitNormal", "[", 
              RowBox[{"Gdd", ",", "coord", ",", "csign"}], "]"}]}], ";", 
            RowBox[{"\[Gamma]dd", "=", 
             RowBox[{"Gdd", "-", 
              RowBox[{"Outer", "[", 
               RowBox[{"Times", ",", "nd", ",", "nd"}], "]"}]}]}], ";", 
            RowBox[{"\[Gamma]ud", "=", 
             RowBox[{
              RowBox[{"IMetric", "[", "Gdd", "]"}], ".", "\[Gamma]dd"}]}], 
            ";", 
            RowBox[{"delN", "=", 
             RowBox[{"CoD", "[", 
              RowBox[{"nd", ",", "Gdd", ",", "X", ",", 
               RowBox[{"{", 
                RowBox[{"-", "1"}], "}"}]}], "]"}]}], ";", 
            RowBox[{"delN", "=", 
             RowBox[{"Simplify", "[", 
              RowBox[{"delN", ",", 
               RowBox[{"TimeConstraint", "\[Rule]", "0.1`"}]}], "]"}]}], ";", 
            RowBox[{"delSym", "=", 
             RowBox[{
              FractionBox["1", "2"], " ", 
              RowBox[{"(", 
               RowBox[{"delN", "+", 
                RowBox[{"Transpose", "[", "delN", "]"}]}], ")"}]}]}], ";", 
            RowBox[{"Kdd", "=", 
             RowBox[{"-", 
              RowBox[{
               RowBox[{"Transpose", "[", "\[Gamma]ud", "]"}], ".", "delSym", 
               ".", "\[Gamma]ud"}]}]}], ";", "Kdd"}]}], "]"}]}]}
      },
      BaselinePosition->{Baseline, {1, 1}},
      GridBoxAlignment->{
       "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
        "RowsIndexed" -> {}},
      GridBoxItemSize->{"Columns" -> {{
           Scaled[0.999]}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, 
        "RowsIndexed" -> {}}]}
   },
   BaselinePosition->{Baseline, {1, 1}},
   GridBoxAlignment->{
    "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
     "RowsIndexed" -> {}}],
  Definition[GREATER2`ExtrinsicCurvature],
  Editable->False]], "Print",
 CellChangeTimes->{3.549982188105831*^9},
 CellTags->"Info3549967787-9563680"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"??", "Metric"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 StyleBox["\<\"\\!\\(\\\"Metric[\\\\!\\\\(\\\\*\\\\nStyleBox[\\\\\\\"ds2\\\\\\\
\",\\\\nFontSlant->\\\\\\\"Italic\\\\\\\"]\\\\), \
{\\\\!\\\\(\\\\*SubscriptBox[\\\\nStyleBox[\\\\\\\"x\\\\\\\",\\\\nFontSlant->\
\\\\\\\"Italic\\\\\\\"], \\\\(1\\\\)]\\\\), \
\\\\!\\\\(\\\\*SubscriptBox[\\\\nStyleBox[\\\\\\\"x\\\\\\\",\\\\nFontSlant->\\\
\\\\\"Italic\\\\\\\"], \\\\(2\\\\)]\\\\), ...}] returns the metric tensor \
specified by the line element \
\\\\!\\\\(\\\\*\\\\nStyleBox[\\\\\\\"ds2\\\\\\\",\\\\nFontSlant->\\\\\\\"\
Italic\\\\\\\"]\\\\) with coordinates \
\\\\!\\\\(\\\\*SubscriptBox[\\\\nStyleBox[\\\\\\\"x\\\\\\\",\\\\nFontSlant->\\\
\\\\\"Italic\\\\\\\"], \\\\(i\\\\)]\\\\). (Just like ds2met but does not \
require you to define differentials.)\\\"\\)\"\>", "MSG"]], "Print", \
"PrintUsage",
 CellChangeTimes->{3.549982184714356*^9},
 CellTags->"Info3549967784-9563680"],

Cell[BoxData[
 InterpretationBox[GridBox[{
    {GridBox[{
       {
        RowBox[{
         RowBox[{"Metric", "[", 
          RowBox[{"ds2_", ",", "X_"}], "]"}], ":=", 
         RowBox[{"Module", "[", 
          RowBox[{
           RowBox[{"{", "dX", "}"}], ",", 
           RowBox[{
            RowBox[{"dX", "=", 
             RowBox[{"Differential", "[", "X", "]"}]}], ";", 
            RowBox[{"ds2met", "[", 
             RowBox[{"ds2", ",", "dX"}], "]"}]}]}], "]"}]}]}
      },
      BaselinePosition->{Baseline, {1, 1}},
      GridBoxAlignment->{
       "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
        "RowsIndexed" -> {}},
      GridBoxItemSize->{"Columns" -> {{
           Scaled[0.999]}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, 
        "RowsIndexed" -> {}}]}
   },
   BaselinePosition->{Baseline, {1, 1}},
   GridBoxAlignment->{
    "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
     "RowsIndexed" -> {}}],
  Definition[GREATER2`Metric],
  Editable->False]], "Print",
 CellChangeTimes->{3.549982184763112*^9},
 CellTags->"Info3549967784-9563680"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{966, 1268},
WindowMargins->{{177, Automatic}, {Automatic, 10}},
ShowSelection->True,
FrontEndVersion->"8.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (October 5, \
2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{
 "Info3549967787-9563680"->{
  Cell[20102, 571, 888, 15, 60, "Print",
   CellTags->"Info3549967787-9563680"],
  Cell[20993, 588, 281, 5, 77, "Print",
   CellTags->"Info3549967787-9563680"],
  Cell[21277, 595, 137, 2, 39, "Print",
   CellTags->"Info3549967787-9563680"],
  Cell[21417, 599, 2466, 61, 89, "Print",
   CellTags->"Info3549967787-9563680"]},
 "Info3549967784-9563680"->{
  Cell[24007, 671, 891, 15, 42, "Print",
   CellTags->"Info3549967784-9563680"],
  Cell[24901, 688, 1119, 31, 20, "Print",
   CellTags->"Info3549967784-9563680"]}
 }
*)
(*CellTagsIndex
CellTagsIndex->{
 {"Info3549967787-9563680", 26390, 735},
 {"Info3549967784-9563680", 26743, 744}
 }
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 2112, 46, 208, "Input"],
Cell[CellGroupData[{
Cell[2682, 70, 2596, 61, 238, "Input"],
Cell[5281, 133, 841, 25, 89, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6159, 163, 2584, 66, 268, "Input"],
Cell[8746, 231, 138, 3, 27, "Output"],
Cell[8887, 236, 269, 8, 30, "Output"],
Cell[9159, 246, 1325, 39, 79, "Output"],
Cell[10487, 287, 1217, 36, 87, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11741, 328, 1677, 44, 133, "Input"],
Cell[13421, 374, 1024, 30, 109, "Output"],
Cell[14448, 406, 288, 6, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14773, 417, 2454, 65, 148, "Input"],
Cell[17230, 484, 478, 13, 33, "Output"],
Cell[17711, 499, 724, 21, 47, "Output"],
Cell[18438, 522, 726, 21, 47, "Output"],
Cell[19167, 545, 225, 4, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19429, 554, 648, 13, 88, "Input"],
Cell[CellGroupData[{
Cell[20102, 571, 888, 15, 60, "Print",
 CellTags->"Info3549967787-9563680"],
Cell[20993, 588, 281, 5, 77, "Print",
 CellTags->"Info3549967787-9563680"],
Cell[21277, 595, 137, 2, 39, "Print",
 CellTags->"Info3549967787-9563680"],
Cell[21417, 599, 2466, 61, 89, "Print",
 CellTags->"Info3549967787-9563680"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[23932, 666, 50, 1, 27, "Input"],
Cell[CellGroupData[{
Cell[24007, 671, 891, 15, 42, "Print",
 CellTags->"Info3549967784-9563680"],
Cell[24901, 688, 1119, 31, 20, "Print",
 CellTags->"Info3549967784-9563680"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)