(* 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[ 13607, 424] NotebookOptionsPosition[ 12913, 399] NotebookOutlinePosition[ 13271, 415] CellTagsIndexPosition[ 13228, 412] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"f", "[", RowBox[{"x_", ",", "y_"}], "]"}], "=", RowBox[{ RowBox[{"x", "^", "2"}], "-", RowBox[{"y", "^", "2"}]}]}]], "Input", CellChangeTimes->{{3.474281245329436*^9, 3.474281251309776*^9}}], Cell[BoxData[ RowBox[{"surf", "=", RowBox[{"Plot3D", "[", RowBox[{ RowBox[{"f", "[", RowBox[{"x", ",", "y"}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "->", RowBox[{"Opacity", "[", "0.5", "]"}]}], ",", RowBox[{"Mesh", "\[Rule]", "False"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.474281258609886*^9, 3.4742812717599*^9}, { 3.47428302390049*^9, 3.474283054587433*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"p", "[", RowBox[{"x_", ",", "y_"}], "]"}], "=", RowBox[{"{", RowBox[{"x", ",", "y", ",", RowBox[{ RowBox[{"x", "^", "2"}], "-", RowBox[{"y", "^", "2"}]}]}], "}"}]}]], "Input", CellChangeTimes->{{3.474281301435609*^9, 3.474281371916804*^9}, { 3.474281444692228*^9, 3.474281445004281*^9}}], Cell[BoxData[ RowBox[{"p", "[", RowBox[{"t", ",", "1"}], "]"}]], "Input", CellChangeTimes->{{3.474281770527296*^9, 3.474281773436019*^9}}], Cell[BoxData[ RowBox[{"Clear", "[", RowBox[{ "x", ",", "y", ",", "t", ",", "y0", ",", "y1", ",", "y2", ",", "y3", ",", "y4", ",", "p", ",", "P"}], "]"}]], "Input", CellChangeTimes->{{3.474281826565684*^9, 3.474281841281253*^9}, { 3.474281892649397*^9, 3.474281894505433*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"p", "[", RowBox[{"x_", ",", "y_"}], "]"}], "=", RowBox[{"{", RowBox[{"x", ",", "y", ",", RowBox[{ RowBox[{"x", "^", "2"}], "-", RowBox[{"y", "^", "2"}]}]}], "}"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y3", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"p", "[", RowBox[{"t", ",", RowBox[{"-", "2"}]}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y4", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"p", "[", RowBox[{"t", ",", RowBox[{"-", "1"}]}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y0", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"p", "[", RowBox[{"t", ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y5", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"p", "[", RowBox[{"t", ",", "1"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y2", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"p", "[", RowBox[{"t", ",", "2"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"ycross", "=", RowBox[{"Show", "[", RowBox[{"y3", ",", "y4", ",", "y2", ",", "y0", ",", "y5", ",", "surf", ",", RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.474281280257199*^9, 3.47428129481238*^9}, { 3.474281374306808*^9, 3.474281399735661*^9}, {3.474281453602179*^9, 3.474281554679594*^9}, {3.474281606093403*^9, 3.474281607159151*^9}, { 3.474281670681953*^9, 3.474281688669382*^9}, {3.474281789142511*^9, 3.474281822193476*^9}, {3.474281900738437*^9, 3.474281953179332*^9}, { 3.474282023176175*^9, 3.474282079518057*^9}, {3.474283064745966*^9, 3.474283094897212*^9}, {3.474283159029897*^9, 3.474283161608619*^9}, { 3.474283236014305*^9, 3.474283241641521*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"x3", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"p", "[", RowBox[{ RowBox[{"-", "2"}], ",", "t"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x4", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"p", "[", RowBox[{ RowBox[{"-", "1"}], ",", "t"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x0", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"p", "[", RowBox[{"0", ",", "t"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"p", "[", RowBox[{"1", ",", "t"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"p", "[", RowBox[{"2", ",", "t"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"xcross", "=", RowBox[{"Show", "[", RowBox[{"x0", ",", "x1", ",", "x2", ",", "x3", ",", "x4", ",", "surf", ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"x", ",", "y", ",", "z"}], "}"}]}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.474282110934413*^9, 3.474282164627339*^9}, { 3.474282272682835*^9, 3.474282292135804*^9}, {3.474283100968727*^9, 3.474283106053118*^9}, {3.474283149120794*^9, 3.474283153511454*^9}, { 3.474283199090425*^9, 3.474283207731473*^9}}], Cell[BoxData[{ RowBox[{"Clear", "[", RowBox[{"x", ",", "y"}], "]"}], "\[IndentingNewLine]", RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{"f", "[", RowBox[{"x", ",", "y"}], "]"}], "\[Equal]", "1"}], ",", "x"}], "]"}]}], "Input", CellChangeTimes->{{3.474282172643755*^9, 3.474282176541153*^9}, { 3.474282206902854*^9, 3.474282264015686*^9}, {3.47428234782468*^9, 3.474282388680599*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"q1", "[", RowBox[{"y_", ",", "z_"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"y", "^", "2"}], "+", "z"}], "]"}]}], ",", "y", ",", "z"}], "}"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"q2", "[", RowBox[{"y_", ",", "z_"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"y", "^", "2"}], "+", "z"}], "]"}], ",", "y", ",", "z"}], "}"}]}]}], "Input", CellChangeTimes->{{3.47428240752329*^9, 3.474282479357873*^9}, { 3.474282642275406*^9, 3.474282755273246*^9}, 3.47428291230563*^9}], Cell[BoxData[{ RowBox[{ RowBox[{"z1", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"q1", "[", RowBox[{"t", ",", RowBox[{"-", "2"}]}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"z2", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"q2", "[", RowBox[{"t", ",", RowBox[{"-", "2"}]}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"z3", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"q1", "[", RowBox[{"t", ",", RowBox[{"-", "1"}]}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"z4", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"q2", "[", RowBox[{"t", ",", RowBox[{"-", "1"}]}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"z5", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"q1", "[", RowBox[{"t", ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"z6", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"q2", "[", RowBox[{"t", ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"z7", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"q1", "[", RowBox[{"t", ",", "1"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"z8", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"q2", "[", RowBox[{"t", ",", "1"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"z9", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"q1", "[", RowBox[{"t", ",", "2"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"z10", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"q2", "[", RowBox[{"t", ",", "2"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Thick"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.474282760246307*^9, 3.474282877162232*^9}, { 3.474282907290495*^9, 3.47428293275113*^9}, {3.474283112782466*^9, 3.474283125311975*^9}}], Cell[BoxData[ RowBox[{"levelsets", "=", RowBox[{"Show", "[", RowBox[{ "surf", ",", "z1", ",", "z2", ",", "z3", ",", "z4", ",", "z5", ",", "z6", ",", "z7", ",", "z8", ",", "z9", ",", "z10", ",", RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.474282881843254*^9, 3.474282898083459*^9}, { 3.474283136652474*^9, 3.474283139059214*^9}, {3.474283221088901*^9, 3.474283221740019*^9}}], Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"surf", ",", "xcross", ",", "ycross", ",", "levelsets"}], "]"}]], "Input", CellChangeTimes->{{3.474283130911656*^9, 3.47428313289455*^9}, { 3.474283173686767*^9, 3.474283180637496*^9}}] }, WindowSize->{640, 623}, WindowMargins->{{24, Automatic}, {Automatic, 36}}, ShowSelection->True, FrontEndVersion->"7.0 for Linux x86 (32-bit) (February 25, 2009)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[545, 20, 240, 7, 32, "Input"], Cell[788, 29, 584, 16, 55, "Input"], Cell[1375, 47, 352, 10, 32, "Input"], Cell[1730, 59, 143, 3, 32, "Input"], Cell[1876, 64, 288, 6, 32, "Input"], Cell[2167, 72, 2735, 77, 165, "Input"], Cell[4905, 151, 2322, 68, 165, "Input"], Cell[7230, 221, 418, 11, 55, "Input"], Cell[7651, 234, 664, 21, 55, "Input"], Cell[8318, 257, 3907, 122, 363, "Input"], Cell[12228, 381, 440, 9, 55, "Input"], Cell[12671, 392, 238, 5, 32, "Input"] } ] *) (* End of internal cache information *)