(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 10391, 365] NotebookOptionsPosition[ 9289, 323] NotebookOutlinePosition[ 9627, 338] CellTagsIndexPosition[ 9584, 335] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"f", "[", RowBox[{"x_", ",", "y__"}], "]"}], ":=", " ", RowBox[{"{", RowBox[{"x", ",", "y", ",", RowBox[{ SuperscriptBox["x", "2"], "-", SuperscriptBox["y", "2"]}]}], "}"}]}]], "Input", CellChangeTimes->{{3.4469777153151093`*^9, 3.44697773307019*^9}, { 3.446978084917082*^9, 3.4469780890015993`*^9}, {3.446978295993376*^9, 3.446978296717963*^9}, {3.446978354012147*^9, 3.44697835437987*^9}, { 3.446982267118168*^9, 3.4469822808327*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["r", "x"], "[", RowBox[{"x_", ",", "y_"}], "]"}], ":=", " ", RowBox[{"D", "[", RowBox[{ RowBox[{"f", "[", RowBox[{"x", ",", "y"}], "]"}], ",", "x"}], "]"}]}]], "Input", CellChangeTimes->{{3.446977738496561*^9, 3.4469777666487083`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["r", "y"], "[", RowBox[{"x_", ",", "y_"}], "]"}], ":=", " ", RowBox[{"D", "[", RowBox[{ RowBox[{"f", "[", RowBox[{"x", ",", "y"}], "]"}], ",", "y"}], "]"}]}]], "Input", CellChangeTimes->{{3.4469777760875797`*^9, 3.4469777803104267`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"g", "=", RowBox[{"(", GridBox[{ { RowBox[{ RowBox[{ SubscriptBox["r", "x"], "[", RowBox[{"x", ",", "y"}], "]"}], " ", ".", " ", RowBox[{ SubscriptBox["r", "x"], "[", RowBox[{"x", ",", "y"}], "]"}]}], RowBox[{ RowBox[{ SubscriptBox["r", "x"], "[", RowBox[{"x", ",", "y"}], "]"}], " ", ".", " ", RowBox[{ SubscriptBox["r", "y"], "[", RowBox[{"x", ",", "y"}], "]"}]}]}, { RowBox[{ RowBox[{ SubscriptBox["r", "x"], "[", RowBox[{"x", ",", "y"}], "]"}], " ", ".", RowBox[{ SubscriptBox["r", "y"], "[", RowBox[{"x", ",", "y"}], "]"}]}], RowBox[{ RowBox[{ SubscriptBox["r", "y"], "[", RowBox[{"x", ",", "y"}], "]"}], ".", " ", RowBox[{ SubscriptBox["r", "y"], "[", RowBox[{"x", ",", "y"}], "]"}]}]} }], ")"}]}]], "Input", CellChangeTimes->{{3.446978422508575*^9, 3.446978484770862*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"1", "+", RowBox[{"4", " ", SuperscriptBox["x", "2"]}]}], ",", RowBox[{ RowBox[{"-", "4"}], " ", "x", " ", "y"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "4"}], " ", "x", " ", "y"}], ",", RowBox[{"1", "+", RowBox[{"4", " ", SuperscriptBox["y", "2"]}]}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.446978485872553*^9, 3.446982285835558*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["r", "xx"], "[", RowBox[{"x_", ",", "y_"}], "]"}], ":=", " ", RowBox[{"D", "[", RowBox[{ RowBox[{"f", "[", RowBox[{"x", ",", "y"}], "]"}], ",", "x", ",", "x"}], "]"}]}]], "Input", CellChangeTimes->{{3.446977946130082*^9, 3.446977958325079*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["r", "xy"], "[", RowBox[{"x_", ",", "y_"}], "]"}], ":=", " ", RowBox[{"D", "[", RowBox[{ RowBox[{"f", "[", RowBox[{"x", ",", "y"}], "]"}], ",", "x", ",", "y"}], "]"}]}]], "Input", CellChangeTimes->{{3.446977967753179*^9, 3.44697797088339*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["r", "yy"], "[", RowBox[{"x_", ",", "y_"}], "]"}], ":=", " ", RowBox[{"D", "[", RowBox[{ RowBox[{"f", "[", RowBox[{"x", ",", "y"}], "]"}], ",", "y", ",", "y"}], "]"}]}]], "Input", CellChangeTimes->{{3.446977976536498*^9, 3.446977981391097*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"norm", "=", " ", RowBox[{"Cross", "[", RowBox[{ RowBox[{ SubscriptBox["r", "x"], "[", RowBox[{"x", ",", "y"}], "]"}], ",", RowBox[{ SubscriptBox["r", "y"], "[", RowBox[{"x", ",", "y"}], "]"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.446977782100934*^9, 3.4469778111640787`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "x"}], ",", RowBox[{"2", " ", "y"}], ",", "1"}], "}"}]], "Output", CellChangeTimes->{{3.4469778060671043`*^9, 3.446977813046978*^9}, 3.446978098096925*^9, 3.446978304605147*^9, 3.446978357582222*^9, 3.446982291140275*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"unitnorm", "=", FractionBox["norm", SqrtBox[ RowBox[{"norm", ".", "norm"}]]]}]], "Input", CellChangeTimes->{{3.446977843710034*^9, 3.446977856832185*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"2", " ", "x"}], SqrtBox[ RowBox[{"1", "+", RowBox[{"4", " ", SuperscriptBox["x", "2"]}], "+", RowBox[{"4", " ", SuperscriptBox["y", "2"]}]}]]]}], ",", FractionBox[ RowBox[{"2", " ", "y"}], SqrtBox[ RowBox[{"1", "+", RowBox[{"4", " ", SuperscriptBox["x", "2"]}], "+", RowBox[{"4", " ", SuperscriptBox["y", "2"]}]}]]], ",", FractionBox["1", SqrtBox[ RowBox[{"1", "+", RowBox[{"4", " ", SuperscriptBox["x", "2"]}], "+", RowBox[{"4", " ", SuperscriptBox["y", "2"]}]}]]]}], "}"}]], "Output", CellChangeTimes->{3.446977857739973*^9, 3.446978099374817*^9, 3.446978307357561*^9, 3.446978359397627*^9, 3.4469822930706472`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ SubscriptBox["r", "xx"], "[", RowBox[{"x", ",", "y"}], "]"}]], "Input", CellChangeTimes->{{3.4469779324871683`*^9, 3.446977933097499*^9}, { 3.446977991988421*^9, 3.446977998943015*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"0", ",", "0", ",", "2"}], "}"}]], "Output", CellChangeTimes->{3.446977999544606*^9, 3.446978101317737*^9, 3.446978310780607*^9, 3.446978361876317*^9, 3.446982294914431*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"L", "=", " ", RowBox[{"(", GridBox[{ { RowBox[{ RowBox[{ SubscriptBox["r", "xx"], "[", RowBox[{"x", ",", "y"}], "]"}], " ", ".", " ", "unitnorm"}], RowBox[{ RowBox[{ SubscriptBox["r", "xy"], "[", RowBox[{"x", ",", "y"}], "]"}], " ", ".", " ", "unitnorm"}]}, { RowBox[{ RowBox[{ SubscriptBox["r", "xy"], "[", RowBox[{"x", ",", "y"}], "]"}], " ", ".", " ", "unitnorm"}], RowBox[{ RowBox[{ SubscriptBox["r", "yy"], "[", RowBox[{"x", ",", "y"}], "]"}], " ", ".", " ", "unitnorm"}]} }], ")"}]}]], "Input", CellChangeTimes->{{3.446978004897355*^9, 3.446978048200841*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ FractionBox["2", SqrtBox[ RowBox[{"1", "+", RowBox[{"4", " ", SuperscriptBox["x", "2"]}], "+", RowBox[{"4", " ", SuperscriptBox["y", "2"]}]}]]], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", FractionBox["2", SqrtBox[ RowBox[{"1", "+", RowBox[{"4", " ", SuperscriptBox["x", "2"]}], "+", RowBox[{"4", " ", SuperscriptBox["y", "2"]}]}]]]}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.446978049030675*^9, 3.4469781033188553`*^9, 3.446978313128008*^9, 3.446978364798625*^9, 3.446982299074934*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"MatrixForm", "[", "L", "]"}]], "Input", CellChangeTimes->{{3.4469780530794764`*^9, 3.446978056122285*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { FractionBox["2", SqrtBox[ RowBox[{"1", "+", RowBox[{"4", " ", SuperscriptBox["x", "2"]}], "+", RowBox[{"4", " ", SuperscriptBox["y", "2"]}]}]]], "0"}, {"0", RowBox[{"-", FractionBox["2", SqrtBox[ RowBox[{"1", "+", RowBox[{"4", " ", SuperscriptBox["x", "2"]}], "+", RowBox[{"4", " ", 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.4469780569062634`*^9, 3.446978118351741*^9, 3.446978367145875*^9, 3.446982300871826*^9}] }, Open ]] }, WindowSize->{640, 656}, WindowMargins->{{4, Automatic}, {Automatic, 4}}, FrontEndVersion->"7.0 for Mac OS X x86 (32-bit) (November 10, 2008)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[545, 20, 505, 12, 33, "Input"], Cell[1053, 34, 308, 9, 27, "Input"], Cell[1364, 45, 310, 9, 30, "Input"], Cell[CellGroupData[{ Cell[1699, 58, 1042, 34, 49, "Input"], Cell[2744, 94, 515, 17, 33, "Output"] }, Open ]], Cell[3274, 114, 317, 9, 27, "Input"], Cell[3594, 125, 316, 9, 30, "Input"], Cell[3913, 136, 317, 9, 30, "Input"], Cell[CellGroupData[{ Cell[4255, 149, 343, 10, 33, "Input"], Cell[4601, 161, 315, 8, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4953, 174, 190, 5, 50, "Input"], Cell[5146, 181, 847, 28, 56, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6030, 214, 217, 5, 27, "Input"], Cell[6250, 221, 216, 4, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6503, 230, 729, 22, 51, "Input"], Cell[7235, 254, 730, 23, 56, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8002, 282, 131, 2, 27, "Input"], Cell[8136, 286, 1137, 34, 89, "Output"] }, Open ]] } ] *) (* End of internal cache information *)