(* 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[ 15341, 447] NotebookOptionsPosition[ 15066, 433] NotebookOutlinePosition[ 15428, 449] CellTagsIndexPosition[ 15385, 446] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ RowBox[{ RowBox[{"Needs", "[", "\"\\"", "]"}], ";"}], "\n", RowBox[{"\"\\"", ";", RowBox[{ RowBox[{"F", "[", RowBox[{"{", RowBox[{"x_", ",", "y_", ",", "z_"}], "}"}], "]"}], ":=", RowBox[{"{", RowBox[{"z", ",", "x", ",", "y"}], "}"}]}], ";", " ", RowBox[{ RowBox[{"CRL", "[", RowBox[{"{", RowBox[{"x_", ",", "y_", ",", "z_"}], "}"}], "]"}], "=", RowBox[{"Curl", "[", RowBox[{ RowBox[{"F", "[", RowBox[{"{", RowBox[{"x", ",", "y", ",", "z"}], "}"}], "]"}], ",", RowBox[{"Cartesian", "[", RowBox[{"x", ",", "y", ",", "z"}], "]"}]}], "]"}]}], ";"}], "\n", RowBox[{"\"\\"", ";", RowBox[{ RowBox[{"f", "[", "v_", "]"}], ":=", RowBox[{"2", " ", RowBox[{"Cos", "[", "v", "]"}]}]}], ";", RowBox[{ RowBox[{"g", "[", "v_", "]"}], ":=", RowBox[{"2", " ", RowBox[{"Sin", "[", "v", "]"}]}]}], ";", RowBox[{ RowBox[{"h", "[", "v_", "]"}], ":=", "0"}], ";", " ", RowBox[{ RowBox[{"r", "[", "v_", "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"f", "[", "v", "]"}], ",", RowBox[{"g", "[", "v", "]"}], ",", RowBox[{"h", "[", "v", "]"}]}], "}"}]}], ";", " ", RowBox[{ RowBox[{"Tn", "[", "v_", "]"}], ":=", RowBox[{"Normalize", "[", RowBox[{ SuperscriptBox["r", "\[Prime]", MultilineFunction->None], "[", "v", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"\"\\"", ";", RowBox[{"a", "=", "0"}], ";", RowBox[{"b", "=", RowBox[{"2", " ", "\[Pi]"}]}], ";", "\"\\"", ";", RowBox[{"n", "=", "35"}], ";", RowBox[{"s", "=", FractionBox[ RowBox[{"b", "-", "a"}], "n"]}], ";", RowBox[{"SC", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"f", "[", "v", "]"}], ",", RowBox[{"g", "[", "v", "]"}], ",", RowBox[{"h", "[", "v", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"v", ",", "a", ",", "b"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{"\"\\"", ";", RowBox[{ RowBox[{"Surf", "[", RowBox[{"u_", ",", "v_"}], "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"u", " ", RowBox[{"f", "[", "v", "]"}]}], ",", RowBox[{"u", " ", RowBox[{"g", "[", "v", "]"}]}], ",", RowBox[{"4", " ", RowBox[{"(", RowBox[{"1", "-", SuperscriptBox["u", "2"]}], ")"}]}]}], "}"}]}], ";", "\"\\"", ";", RowBox[{"Au", "=", "0"}], ";", RowBox[{"Bu", "=", "1"}], ";", RowBox[{"Av", "=", "a"}], ";", RowBox[{"Bv", "=", "b"}], ";", "\"\\"", ";", RowBox[{"mu", "=", "10"}], ";", RowBox[{"mv", "=", "10"}], ";", RowBox[{"Su", "=", FractionBox[ RowBox[{"Bu", "-", "Au"}], "mu"]}], ";", RowBox[{"Sv", "=", FractionBox[ RowBox[{"Bv", "-", "Av"}], "mv"]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"SD", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"Surf", "[", RowBox[{"u", ",", "v"}], "]"}], ",", RowBox[{"{", RowBox[{"u", ",", "Au", ",", "Bu"}], "}"}], ",", RowBox[{"{", RowBox[{"v", ",", "Av", ",", "Bv"}], "}"}]}], "]"}]}], ";", RowBox[{ RowBox[{"Nr", "[", RowBox[{"u_", ",", "v_"}], "]"}], "=", RowBox[{"-", RowBox[{"Normalize", "[", RowBox[{ RowBox[{ SubscriptBox["\[PartialD]", "v"], RowBox[{"Surf", "[", RowBox[{"u", ",", "v"}], "]"}]}], "\[Cross]", RowBox[{ SubscriptBox["\[PartialD]", "u"], RowBox[{"Surf", "[", RowBox[{"u", ",", "v"}], "]"}]}]}], "]"}]}]}], ";", "\"\\"", ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Integrand", "=", RowBox[{"Simplify", "[", RowBox[{ RowBox[{ RowBox[{"F", "[", RowBox[{"r", "[", "v", "]"}], "]"}], ".", RowBox[{ SuperscriptBox["r", "\[Prime]", MultilineFunction->None], "[", "v", "]"}]}], ",", RowBox[{"Assumptions", "\[Rule]", RowBox[{"a", "\[LessEqual]", "v", "\[LessEqual]", "b"}]}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"Integrand1", "=", RowBox[{"Simplify", "[", RowBox[{ RowBox[{"-", RowBox[{ RowBox[{"CRL", "[", RowBox[{"Surf", "[", RowBox[{"u", ",", "v"}], "]"}], "]"}], ".", " ", RowBox[{ RowBox[{ SubscriptBox["\[PartialD]", "v"], RowBox[{"Surf", "[", RowBox[{"u", ",", "v"}], "]"}]}], "\[Cross]", RowBox[{ SubscriptBox["\[PartialD]", "u"], RowBox[{"Surf", "[", RowBox[{"u", ",", "v"}], "]"}]}]}]}]}], ",", RowBox[{"Assumptions", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"Av", "\[LessEqual]", "v", "\[LessEqual]", "Bv"}], "&&", RowBox[{"Au", "\[LessEqual]", "u", "\[LessEqual]", "Bu"}]}], "}"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"LI", "=", RowBox[{ SubsuperscriptBox["\[Integral]", "a", "b"], RowBox[{ RowBox[{ RowBox[{"F", "[", RowBox[{"r", "[", "v", "]"}], "]"}], ".", RowBox[{ SuperscriptBox["r", "\[Prime]", MultilineFunction->None], "[", "v", "]"}]}], RowBox[{"\[DifferentialD]", "v"}]}]}]}], ";", " ", RowBox[{"CFI", "=", RowBox[{ SubsuperscriptBox["\[Integral]", "Au", "Bu"], RowBox[{ SubsuperscriptBox["\[Integral]", "Av", "Bv"], RowBox[{ RowBox[{"-", RowBox[{ RowBox[{"CRL", "[", RowBox[{"Surf", "[", RowBox[{"u", ",", "v"}], "]"}], "]"}], ".", RowBox[{ RowBox[{ SubscriptBox["\[PartialD]", "v"], RowBox[{"Surf", "[", RowBox[{"u", ",", "v"}], "]"}]}], "\[Cross]", RowBox[{ SubscriptBox["\[PartialD]", "u"], RowBox[{"Surf", "[", RowBox[{"u", ",", "v"}], "]"}]}]}]}]}], RowBox[{"\[DifferentialD]", "v"}], RowBox[{"\[DifferentialD]", "u"}]}]}]}]}], ";"}], "\n", RowBox[{ RowBox[{"P1", "=", RowBox[{"Flatten", "[", RowBox[{ RowBox[{"Table", "[", " ", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{"Green", ",", RowBox[{"Arrowheads", "\[Rule]", "Small"}], ",", RowBox[{"Arrow", "[", RowBox[{"{", " ", RowBox[{ RowBox[{"r", "[", RowBox[{"a", "+", RowBox[{"s", " ", "i"}]}], "]"}], ",", RowBox[{ RowBox[{"r", "[", RowBox[{"a", "+", RowBox[{"s", " ", "i"}]}], "]"}], "+", RowBox[{"Tn", "[", RowBox[{"a", "+", RowBox[{"s", " ", "i"}]}], "]"}]}]}], "}"}], "]"}]}], " ", "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}], ",", "2"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"P2", "=", RowBox[{"Flatten", "[", RowBox[{ RowBox[{"Table", "[", " ", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{"Blue", ",", RowBox[{"Arrowheads", "\[Rule]", "Small"}], ",", RowBox[{"Arrow", "[", RowBox[{"{", " ", RowBox[{ RowBox[{"r", "[", RowBox[{"a", "+", RowBox[{"s", " ", "i"}]}], "]"}], ",", RowBox[{ RowBox[{"r", "[", RowBox[{"a", "+", RowBox[{"s", " ", "i"}]}], "]"}], "+", " ", RowBox[{"0.2", "*", RowBox[{"F", "[", RowBox[{"r", "[", RowBox[{"a", "+", RowBox[{"s", " ", "i"}]}], "]"}], "]"}]}]}]}], "}"}], "]"}]}], " ", "}"}], " ", "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}], ",", "2"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"P3", "=", RowBox[{"Flatten", "[", RowBox[{ RowBox[{"Table", "[", " ", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{"Magenta", ",", RowBox[{"Arrowheads", "\[Rule]", "Small"}], ",", "\[IndentingNewLine]", RowBox[{"Arrow", "[", RowBox[{"{", " ", RowBox[{ RowBox[{"Surf", "[", RowBox[{ RowBox[{"Au", "+", RowBox[{"Su", " ", "i"}]}], ",", RowBox[{"Av", "+", RowBox[{"Sv", " ", "j"}]}]}], "]"}], ",", RowBox[{ RowBox[{"Surf", "[", RowBox[{ RowBox[{"Au", "+", RowBox[{"Su", " ", "i"}]}], ",", RowBox[{"Av", "+", RowBox[{"Sv", " ", "j"}]}]}], "]"}], "+", RowBox[{"CRL", "[", RowBox[{"Surf", "[", RowBox[{ RowBox[{"Au", "+", RowBox[{"Su", " ", "i"}]}], ",", RowBox[{"Av", "+", RowBox[{"Sv", " ", "j"}]}]}], "]"}], "]"}]}]}], "}"}], "]"}]}], " ", "}"}], " ", "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "mu"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "mv"}], "}"}]}], "]"}], ",", "2"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"P4", "=", RowBox[{"Flatten", "[", RowBox[{ RowBox[{"Table", "[", " ", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{"Red", ",", RowBox[{"Arrowheads", "\[Rule]", "Small"}], ",", "\[IndentingNewLine]", RowBox[{"Arrow", "[", RowBox[{"{", " ", RowBox[{ RowBox[{"Surf", "[", RowBox[{ RowBox[{"Au", "+", RowBox[{"Su", " ", "i"}]}], ",", RowBox[{"Av", "+", RowBox[{"Sv", " ", "j"}]}]}], "]"}], ",", RowBox[{ RowBox[{"Surf", "[", RowBox[{ RowBox[{"Au", "+", RowBox[{"Su", " ", "i"}]}], ",", RowBox[{"Av", "+", RowBox[{"Sv", " ", "j"}]}]}], "]"}], "+", " ", RowBox[{"Nr", "[", RowBox[{ RowBox[{"Au", "+", RowBox[{"Su", " ", "i"}]}], ",", RowBox[{"Av", "+", RowBox[{"Sv", " ", "j"}]}]}], "]"}]}]}], "}"}], "]"}]}], " ", "}"}], " ", "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "mu"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "mv"}], "}"}]}], "]"}], ",", "2"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"g1", "=", RowBox[{"Show", "[", RowBox[{"P1", ",", "P2", ",", "SC", ",", RowBox[{"ViewPoint", "\[Rule]", RowBox[{"{", RowBox[{"1.3`", ",", RowBox[{"-", "2.4`"}], ",", "2.`"}], "}"}]}]}], "]"}]}], ";", RowBox[{"g2", "=", RowBox[{"Show", "[", RowBox[{"P1", ",", "P2", ",", "P3", ",", "P4", ",", "SC", ",", "SD", ",", RowBox[{"ViewPoint", "\[Rule]", RowBox[{"{", RowBox[{"1.3`", ",", RowBox[{"-", "2.4`"}], ",", "2.`"}], "}"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"Style", "[", "\"\<\!\(\* StyleBox[\"Curve\", FontColor->GrayLevel[0]]\) \!\(\* StyleBox[\"Unit\", FontColor->RGBColor[0, 1, 0]]\)\!\(\* StyleBox[\" \", FontColor->RGBColor[0, 1, 0]]\)\!\(\* StyleBox[\"Tangent\", FontColor->RGBColor[0, 1, 0]]\) \!\(\* StyleBox[\"Vector\", FontColor->RGBColor[0, 0, 1]]\)\!\(\* StyleBox[\" \", FontColor->RGBColor[0, 0, 1]]\)\!\(\* StyleBox[\"Field\", FontColor->RGBColor[0, 0, 1]]\) \!\(\* StyleBox[\"Unit\", FontColor->RGBColor[1, 0, 0]]\)\!\(\* StyleBox[\" \", FontColor->RGBColor[1, 0, 0]]\)\!\(\* StyleBox[\"Normal\", FontColor->RGBColor[1, 0, 0]]\) \!\(\* StyleBox[\"Curl\", FontColor->RGBColor[1, 0, 1]]\)\>\"", "]"}], "]"}], ";", RowBox[{"Print", "[", RowBox[{ RowBox[{ "Style", "[", "\"\< \ Work=\[Integral]F\[CenterDot]Tds=\[Integral]F\[CenterDot]dr=\[Integral]\>\"", "]"}], ",", "Integrand", ",", "\"\<= \>\"", ",", "LI"}], "]"}], ";"}], "\n", RowBox[{ RowBox[{"Print", "[", RowBox[{ RowBox[{"Style", "[", "\"\\"", "]"}], ",", RowBox[{"CRL", "[", RowBox[{"{", RowBox[{"x", ",", "y", ",", "z"}], "}"}], "]"}], ",", "\"\<; \>\"", ",", RowBox[{ "Style", "[", "\"\\"", "]"}], ",", "Integrand1", ",", "\"\<= \>\"", ",", "CFI"}], "]"}], ";"}], "\n", RowBox[{"Show", "[", RowBox[{ RowBox[{"GraphicsRow", "[", RowBox[{"{", RowBox[{"g1", ",", "g2"}], "}"}], "]"}], ",", RowBox[{"ImageSize", "\[Rule]", RowBox[{"{", RowBox[{"1000", ",", "400"}], "}"}]}]}], "]"}], "\n"}], "Input", CellChangeTimes->{{3.4118623466875*^9, 3.41186247540625*^9}, 3.41229727090625*^9, 3.47248997459375*^9, {3.472490011203125*^9, 3.472490252390625*^9}, {3.472490300234375*^9, 3.472490508875*^9}, { 3.47249060459375*^9, 3.47249073703125*^9}, {3.472490841703125*^9, 3.47249088765625*^9}, 3.472490927875*^9, {3.47249097315625*^9, 3.472491081921875*^9}, 3.4724911194375*^9, {3.472491152859375*^9, 3.472491163171875*^9}, {3.472491277359375*^9, 3.472491355234375*^9}, { 3.47249200146875*^9, 3.472492072921875*^9}, {3.472492109765625*^9, 3.47249211521875*^9}, {3.472492328203125*^9, 3.472492332359375*^9}, { 3.472492433765625*^9, 3.472492476984375*^9}, {3.472492527875*^9, 3.47249257596875*^9}, 3.47249294475*^9, 3.472492980015625*^9, { 3.472493060140625*^9, 3.47249306890625*^9}, {3.472493105921875*^9, 3.472493302703125*^9}, {3.472493346609375*^9, 3.472493347453125*^9}, { 3.472493394953125*^9, 3.47249359940625*^9}, {3.47249364134375*^9, 3.4724937189375*^9}, {3.47249375975*^9, 3.472493762046875*^9}, { 3.472493809234375*^9, 3.4724938101875*^9}, 3.472514348453125*^9, { 3.472586791375*^9, 3.47258681909375*^9}}] }, WindowSize->{1272, 922}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, Magnification->1, FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (February 18, 2009)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[545, 20, 14517, 411, 492, "Input"] } ] *) (* End of internal cache information *)