(* 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[ 23795, 722] NotebookOptionsPosition[ 19408, 588] NotebookOutlinePosition[ 22191, 668] CellTagsIndexPosition[ 22148, 665] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Green's function solution of Heat Equation 01-11-11", "Title", CellChangeTimes->{{3.4176063396294003`*^9, 3.4176063432174*^9}, { 3.432688554589577*^9, 3.4326885915025773`*^9}, {3.432988387002*^9, 3.432988392038*^9}, {3.4415437730018*^9, 3.4415437945922003`*^9}, { 3.4415542865534*^9, 3.4415542919384003`*^9}, {3.4680660053921576`*^9, 3.4680660321773577`*^9}, {3.468066662083358*^9, 3.4680666665449576`*^9}, { 3.468540050681733*^9, 3.4685400628837337`*^9}, {3.4685407328117332`*^9, 3.4685407374897337`*^9}, {3.47243380830192*^9, 3.47243381470492*^9}, { 3.4726800180721188`*^9, 3.4726800215509186`*^9}, {3.4864996693814096`*^9, 3.4864996740926094`*^9}, {3.486552929915*^9, 3.486552932888*^9}, { 3.5006763285975285`*^9, 3.5006763323259354`*^9}, {3.500749088940343*^9, 3.500749092723559*^9}, {3.5007675028238792`*^9, 3.5007675160366344`*^9}, { 3.5007681974356084`*^9, 3.5007681991157045`*^9}, {3.5035088220813575`*^9, 3.5035088355285807`*^9}, {3.5035813196375813`*^9, 3.503581322296734*^9}, { 3.503780031894681*^9, 3.5037800537659197`*^9}, 3.690566942172567*^9}], Cell["N. T. Gladd", "Subtitle", CellChangeTimes->{{3.6914351493330326`*^9, 3.691435164695605*^9}}], Cell[TextData[{ StyleBox["Initialization: ", FontWeight->"Bold"], "Be sure the files ", StyleBox["NTGStylesheet2.nb", FontSlant->"Italic"], " and ", StyleBox["NTGUtilityFunctions.m", FontSlant->"Italic"], " is are in the same directory as that from which this notebook was loaded. \ Then execute the cell immediately below by mousing left on the cell bar to \ the right of that cell and then typing \[OpenCurlyDoubleQuote]shift\ \[CloseCurlyDoubleQuote] + \[OpenCurlyDoubleQuote]enter\ \[CloseCurlyDoubleQuote]. Respond \[OpenCurlyDoubleQuote]Yes\ \[CloseCurlyDoubleQuote] in response to the query to evaluate initialization \ cells." }], "Text", CellChangeTimes->{{3.69023148371922*^9, 3.6902314907825966`*^9}, { 3.690231568228395*^9, 3.690231573307065*^9}, {3.690231629378787*^9, 3.6902317034648433`*^9}, {3.6902318285880823`*^9, 3.6902318790622425`*^9}, {3.690561803160874*^9, 3.690561837304328*^9}, 3.691435129792831*^9, 3.6914351763077483`*^9}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"SetDirectory", "[", RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}], " ", RowBox[{"(*", " ", RowBox[{ "set", " ", "directory", " ", "where", " ", "source", " ", "files", " ", "are", " ", "located"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"SetOptions", "[", RowBox[{ RowBox[{"EvaluationNotebook", "[", "]"}], ",", " ", RowBox[{"(*", " ", RowBox[{"load", " ", "the", " ", "StyleSheet"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"StyleDefinitions", " ", "\[Rule]", " ", RowBox[{"Get", "[", "\"\\"", "]"}]}]}], "]"}], ";", RowBox[{"Get", "[", "\"\\"", "]"}], ";", " ", RowBox[{"(*", " ", RowBox[{"Load", " ", "utilities", " ", "package"}], " ", "*)"}]}]}], "Input", InitializationCell->True, CellChangeTimes->{{3.690231304262511*^9, 3.6902313300154333`*^9}, { 3.6902313821303463`*^9, 3.6902313878808374`*^9}, {3.690231423337756*^9, 3.6902314287291145`*^9}, {3.6902315145507984`*^9, 3.6902315306304874`*^9}, {3.6905614204501295`*^9, 3.6905616087617702`*^9}, {3.690561669084352*^9, 3.6905616745378094`*^9}, { 3.690561776716816*^9, 3.6905617961927767`*^9}, 3.6911976801547356`*^9}], Cell[CellGroupData[{ Cell["Solving heat equation using Green\[CloseCurlyQuote]s function \ approach", "Section", CellChangeTimes->{{3.6905669741488705`*^9, 3.6905669863611994`*^9}, { 3.6905672127726564`*^9, 3.690567215535365*^9}}], Cell[TextData[{ "This notebook is a revised version of early work \[LongDash] ", StyleBox["Green' s function - Heat equation 02 - 05 - 04. nb", FontSlant->"Italic"] }], "Text", CellChangeTimes->{{3.5037955482276835`*^9, 3.5037955669277534`*^9}}], Cell["Define the basic PDE.", "Text", CellChangeTimes->{{3.5037955848047757`*^9, 3.5037955935462756`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"w", "[", "1", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"G", "[", RowBox[{"x", ",", "t"}], "]"}], ",", " ", "t"}], "]"}], " ", "-", " ", RowBox[{"\[Kappa]", " ", RowBox[{"D", "[", RowBox[{ RowBox[{"G", "[", RowBox[{"x", ",", "t"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"x", ",", " ", "2"}], "}"}]}], "]"}]}]}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"DiracDelta", "[", "t", "]"}], RowBox[{"DiracDelta", "[", RowBox[{"x", " ", "-", " ", "x0"}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.503779612066744*^9, 3.5037796333295813`*^9}, { 3.50379589259838*^9, 3.5037959106484127`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["G", TagBox[ RowBox[{"(", RowBox[{"0", ",", "1"}], ")"}], Derivative], MultilineFunction->None], "[", RowBox[{"x", ",", "t"}], "]"}], "-", RowBox[{"\[Kappa]", " ", RowBox[{ SuperscriptBox["G", TagBox[ RowBox[{"(", RowBox[{"2", ",", "0"}], ")"}], Derivative], MultilineFunction->None], "[", RowBox[{"x", ",", "t"}], "]"}]}]}], "\[Equal]", RowBox[{ RowBox[{"DiracDelta", "[", "t", "]"}], " ", RowBox[{"DiracDelta", "[", RowBox[{"x", "-", "x0"}], "]"}]}]}]], "Output", CellChangeTimes->{ 3.5037796471356053`*^9, 3.5037955974835005`*^9, {3.503795897510661*^9, 3.5037959125425215`*^9}, 3.6905670057736287`*^9}] }, Open ]], Cell["Easier to work with an expression than with an equation", "Text", CellChangeTimes->{{3.5037956047339153`*^9, 3.503795618517704*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"w", "[", "2", "]"}], " ", "=", " ", RowBox[{ RowBox[{"w", "[", "1", "]"}], " ", "/.", " ", RowBox[{"Equal", " ", "\[Rule]", " ", "Subtract"}]}]}]], "Input", CellChangeTimes->{{3.503779686104474*^9, 3.503779702219302*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"DiracDelta", "[", "t", "]"}]}], " ", RowBox[{"DiracDelta", "[", RowBox[{"x", "-", "x0"}], "]"}]}], "+", RowBox[{ SuperscriptBox["G", TagBox[ RowBox[{"(", RowBox[{"0", ",", "1"}], ")"}], Derivative], MultilineFunction->None], "[", RowBox[{"x", ",", "t"}], "]"}], "-", RowBox[{"\[Kappa]", " ", RowBox[{ SuperscriptBox["G", TagBox[ RowBox[{"(", RowBox[{"2", ",", "0"}], ")"}], Derivative], MultilineFunction->None], "[", RowBox[{"x", ",", "t"}], "]"}]}]}]], "Output", CellChangeTimes->{3.5037797029369035`*^9, 3.5037956215428767`*^9, 3.503795914829652*^9, 3.6905670123211927`*^9}] }, Open ]], Cell["Laplace transform with respect to time", "Text", CellChangeTimes->{{3.5037956310514207`*^9, 3.503795640559965*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"w", "[", "3", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"LaplaceTransform", "[", RowBox[{ RowBox[{"w", "[", "2", "]"}], ",", " ", "t", ",", " ", "\[Omega]"}], "]"}], " ", "/.", " ", RowBox[{ RowBox[{"LaplaceTransform", "[", RowBox[{ RowBox[{"G", "[", RowBox[{"x", ",", "t"}], "]"}], ",", "t", ",", "\[Omega]"}], "]"}], " ", "\[Rule]", " ", RowBox[{ OverscriptBox["G", "^"], "[", RowBox[{"x", ",", " ", "\[Omega]"}], "]"}]}]}], " ", "/.", " ", RowBox[{ RowBox[{"LaplaceTransform", "[", RowBox[{ RowBox[{ SuperscriptBox["G", TagBox[ RowBox[{"(", RowBox[{"2", ",", "0"}], ")"}], Derivative], MultilineFunction->None], "[", RowBox[{"x", ",", "t"}], "]"}], ",", "t", ",", "\[Omega]"}], "]"}], " ", "\[Rule]", " ", RowBox[{"D", "[", RowBox[{ RowBox[{ OverscriptBox["G", "^"], "[", RowBox[{"x", ",", " ", "\[Omega]"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"x", ",", " ", "2"}], "}"}]}], "]"}]}]}], " ", "/.", " ", RowBox[{ RowBox[{"G", "[", RowBox[{"x", ",", "0"}], "]"}], " ", "\[Rule]", " ", "0"}]}]}]], "Input", CellChangeTimes->{{3.5037797071645107`*^9, 3.503779753246992*^9}, { 3.5037798447567525`*^9, 3.5037798767212086`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"-", RowBox[{"DiracDelta", "[", RowBox[{"x", "-", "x0"}], "]"}]}], "+", RowBox[{"\[Omega]", " ", RowBox[{ OverscriptBox["G", "^"], "[", RowBox[{"x", ",", "\[Omega]"}], "]"}]}], "-", RowBox[{"\[Kappa]", " ", RowBox[{ SuperscriptBox[ OverscriptBox["G", "^"], TagBox[ RowBox[{"(", RowBox[{"2", ",", "0"}], ")"}], Derivative], MultilineFunction->None], "[", RowBox[{"x", ",", "\[Omega]"}], "]"}]}]}]], "Output", CellChangeTimes->{ 3.5037797183497305`*^9, 3.503779754416994*^9, {3.503779852385166*^9, 3.50377987732961*^9}, 3.5037956452402325`*^9, 3.5037959172297893`*^9, 3.690567019540705*^9}] }, Open ]], Cell["Fourier transform with respect to space", "Text", CellChangeTimes->{{3.5037956540917387`*^9, 3.50379566198219*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"w", "[", "4", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"FourierTransform", "[", RowBox[{"#", ",", " ", "x", ",", " ", "k"}], "]"}], ")"}], "&"}], " ", "/@", " ", RowBox[{"w", "[", "3", "]"}], " "}]}]], "Input", CellChangeTimes->{{3.5037797632310095`*^9, 3.5037798007022753`*^9}, 3.50377990019925*^9}], Cell[BoxData[ RowBox[{ RowBox[{"-", FractionBox[ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "k", " ", "x0"}]], SqrtBox[ RowBox[{"2", " ", "\[Pi]"}]]]}], "+", RowBox[{ SuperscriptBox["k", "2"], " ", "\[Kappa]", " ", RowBox[{"FourierTransform", "[", RowBox[{ RowBox[{ OverscriptBox["G", "^"], "[", RowBox[{"x", ",", "\[Omega]"}], "]"}], ",", "x", ",", "k"}], "]"}]}], "+", RowBox[{"FourierTransform", "[", RowBox[{ RowBox[{"\[Omega]", " ", RowBox[{ OverscriptBox["G", "^"], "[", RowBox[{"x", ",", "\[Omega]"}], "]"}]}], ",", "x", ",", "k"}], "]"}]}]], "Output", CellChangeTimes->{{3.503779777551834*^9, 3.503779801716277*^9}, { 3.5037798864400253`*^9, 3.503779901462852*^9}, 3.503795665971418*^9, 3.5037959229581165`*^9, 3.690567030151198*^9}] }, Open ]], Cell["\<\ This has to be manipulated to remove the \[Omega] from FourierTransform. I \ use a rule\ \>", "Text", CellChangeTimes->{{3.5037801596121054`*^9, 3.5037801636369123`*^9}, { 3.50378023602104*^9, 3.503780242978652*^9}, {3.503795677994106*^9, 3.5037957106609745`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"w", "[", "5", "]"}], " ", "=", " ", RowBox[{ RowBox[{"w", "[", "4", "]"}], " ", "/.", " ", RowBox[{ RowBox[{ RowBox[{"FourierTransform", "[", RowBox[{ RowBox[{"a_.", RowBox[{ OverscriptBox["G", "^"], "[", RowBox[{"x", ",", "\[Omega]"}], "]"}]}], ",", "x", ",", "k"}], "]"}], " ", "/;", " ", RowBox[{"FreeQ", "[", RowBox[{"a", ",", "x"}], "]"}]}], " ", "\[Rule]", " ", RowBox[{"a", " ", RowBox[{"FourierTransform", "[", RowBox[{ RowBox[{ OverscriptBox["G", "^"], "[", RowBox[{"x", ",", "\[Omega]"}], "]"}], ",", "x", ",", "k"}], "]"}]}]}]}]}]], "Input", CellChangeTimes->{{3.5037801649785147`*^9, 3.503780224087019*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"-", FractionBox[ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "k", " ", "x0"}]], SqrtBox[ RowBox[{"2", " ", "\[Pi]"}]]]}], "+", RowBox[{ SuperscriptBox["k", "2"], " ", "\[Kappa]", " ", RowBox[{"FourierTransform", "[", RowBox[{ RowBox[{ OverscriptBox["G", "^"], "[", RowBox[{"x", ",", "\[Omega]"}], "]"}], ",", "x", ",", "k"}], "]"}]}], "+", RowBox[{"\[Omega]", " ", RowBox[{"FourierTransform", "[", RowBox[{ RowBox[{ OverscriptBox["G", "^"], "[", RowBox[{"x", ",", "\[Omega]"}], "]"}], ",", "x", ",", "k"}], "]"}]}]}]], "Output", CellChangeTimes->{3.503780226863824*^9, 3.5037959307125607`*^9, 3.6905670435588627`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"w", "[", "6", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"w", "[", "5", "]"}], " ", "\[Equal]", " ", "0"}], " ", "/.", " ", RowBox[{ RowBox[{"FourierTransform", "[", RowBox[{ RowBox[{ OverscriptBox["G", "^"], "[", RowBox[{"x", ",", "\[Omega]"}], "]"}], ",", "x", ",", "k"}], "]"}], " ", "\[Rule]", " ", RowBox[{ OverscriptBox["G", "~"], "[", RowBox[{"k", ",", "\[Omega]"}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.503779935080911*^9, 3.5037799940178146`*^9}, { 3.503780249905064*^9, 3.503780265193091*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"-", FractionBox[ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "k", " ", "x0"}]], SqrtBox[ RowBox[{"2", " ", "\[Pi]"}]]]}], "+", RowBox[{ SuperscriptBox["k", "2"], " ", "\[Kappa]", " ", RowBox[{ OverscriptBox["G", "~"], "[", RowBox[{"k", ",", "\[Omega]"}], "]"}]}], "+", RowBox[{"\[Omega]", " ", RowBox[{ OverscriptBox["G", "~"], "[", RowBox[{"k", ",", "\[Omega]"}], "]"}]}]}], "\[Equal]", "0"}]], "Output", CellChangeTimes->{ 3.5037799956246176`*^9, {3.503780256316675*^9, 3.503780274085107*^9}, 3.503795933486719*^9, 3.6905670569369435`*^9}] }, Open ]], Cell["Solve the transformed dependent variable.", "Text", CellChangeTimes->{{3.5037957465470266`*^9, 3.5037957602838125`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"w", "[", "7", "]"}], " ", "=", " ", RowBox[{ RowBox[{"Solve", "[", RowBox[{ RowBox[{"w", "[", "6", "]"}], " ", ",", " ", RowBox[{ OverscriptBox["G", "~"], "[", RowBox[{"k", ",", "\[Omega]"}], "]"}]}], "]"}], "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}]}]], "Input", CellChangeTimes->{{3.503780096135594*^9, 3.5037801265244474`*^9}, { 3.5037802766903114`*^9, 3.503780295316744*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ OverscriptBox["G", "~"], "[", RowBox[{"k", ",", "\[Omega]"}], "]"}], "\[Rule]", FractionBox[ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "k", " ", "x0"}]], RowBox[{ SqrtBox[ RowBox[{"2", " ", "\[Pi]"}]], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["k", "2"], " ", "\[Kappa]"}], "+", "\[Omega]"}], ")"}]}]]}]], "Output", CellChangeTimes->{{3.503780104169608*^9, 3.5037801268052483`*^9}, { 3.5037802600918818`*^9, 3.5037802960967455`*^9}, 3.503795939117041*^9, 3.6905670685631313`*^9}] }, Open ]], Cell["First, invert the Laplace transform", "Text", CellChangeTimes->{{3.503780314769978*^9, 3.5037803228039923`*^9}, { 3.503795852150067*^9, 3.5037958525450897`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"w", "[", "8", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"InverseLaplaceTransform", "[", RowBox[{"#", ",", " ", "\[Omega]", ",", " ", "t"}], "]"}], ")"}], "&"}], " ", "/@", " ", RowBox[{"w", "[", "7", "]"}]}], " ", "/.", " ", RowBox[{ RowBox[{"InverseLaplaceTransform", "[", RowBox[{ RowBox[{ OverscriptBox["G", "~"], "[", RowBox[{"k", ",", "\[Omega]"}], "]"}], ",", "\[Omega]", ",", "t"}], "]"}], " ", "\[Rule]", " ", RowBox[{ OverscriptBox["G", "_"], "[", RowBox[{"k", ",", "t"}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.5037803480292363`*^9, 3.503780420538164*^9}, { 3.5037804893654847`*^9, 3.5037804921734896`*^9}, {3.503795800543115*^9, 3.50379580411932*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ OverscriptBox["G", "_"], "[", RowBox[{"k", ",", "t"}], "]"}], "\[Rule]", FractionBox[ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"\[ImaginaryI]", " ", "k", " ", "x0"}], "-", RowBox[{ SuperscriptBox["k", "2"], " ", "t", " ", "\[Kappa]"}]}]], SqrtBox[ RowBox[{"2", " ", "\[Pi]"}]]]}]], "Output", CellChangeTimes->{3.5037803412588243`*^9, 3.5037804936554923`*^9, 3.5037959487825937`*^9, 3.690567075641757*^9}] }, Open ]], Cell["where I'm sort of running short of overscripts for notation.", "Text", CellChangeTimes->{{3.5037958249875135`*^9, 3.503795839023316*^9}}], Cell["Then, invert the Fourier transform", "Text", CellChangeTimes->{{3.5037805086471186`*^9, 3.5037805175235343`*^9}, { 3.5037958439235964`*^9, 3.5037958453816795`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{" ", RowBox[{ RowBox[{"w", "[", "9", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"InverseFourierTransform", "[", RowBox[{"#", ",", " ", "k", ",", " ", "x"}], "]"}], ")"}], "&"}], " ", "/@", " ", RowBox[{"w", "[", "8", "]"}]}], " ", "/.", " ", RowBox[{ RowBox[{"InverseFourierTransform", "[", RowBox[{ RowBox[{ OverscriptBox["G", "_"], "[", RowBox[{"k", ",", "t"}], "]"}], ",", "k", ",", "x"}], "]"}], " ", "\[Rule]", " ", RowBox[{"G", "[", RowBox[{"x", ",", "t"}], "]"}]}]}], " ", "/.", " ", RowBox[{"Rule", " ", "\[Rule]", " ", "Equal", " "}]}]}]}]], "Input", CellChangeTimes->{{3.5037805379907703`*^9, 3.50378057795804*^9}, { 3.5037959649745197`*^9, 3.503795987073784*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"G", "[", RowBox[{"x", ",", "t"}], "]"}], "\[Equal]", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ SuperscriptBox[ RowBox[{"(", RowBox[{"x", "-", "x0"}], ")"}], "2"], RowBox[{"4", " ", "t", " ", "\[Kappa]"}]]}]], " ", SqrtBox[ FractionBox["1", RowBox[{"t", " ", "\[Kappa]"}]]]}], RowBox[{"2", " ", SqrtBox["\[Pi]"]}]]}]], "Output", CellChangeTimes->{{3.5037805329519615`*^9, 3.503780578301241*^9}, { 3.503795968729735*^9, 3.5037959877168207`*^9}, 3.6905670885821466`*^9}] }, Open ]], Cell["\<\ and we have the Green's function solution of the original pde\ \>", "Text", CellChangeTimes->{{3.503795861238587*^9, 3.503795876308449*^9}, { 3.503795995168247*^9, 3.5037960374156637`*^9}, {3.6905671167259035`*^9, 3.6905671183196907`*^9}, {3.690567156126048*^9, 3.6905671681741514`*^9}}] }, Open ]] }, Open ]] }, WindowSize->{1085, 657}, WindowMargins->{{37, Automatic}, {Automatic, 0}}, PrintingCopies->1, PrintingPageRange->{32000, 32000}, PrintingOptions->{"Magnification"->1., "PaperOrientation"->"Portrait", "PaperSize"->{612, 792}}, ShowSelection->True, Magnification:>1.5 Inherited, FrontEndVersion->"11.0 for Microsoft Windows (64-bit) (July 28, 2016)", StyleDefinitions->Notebook[{ Cell[ StyleData[ StyleDefinitions -> FrontEnd`FileName[{"Report"}, "StandardReport.nb", CharacterEncoding -> "WindowsANSI"]]], Cell[ StyleData["Notebook"], MenuSortingValue -> 10000, Background -> RGBColor[1., 1., 0.878431]], Cell[ StyleData["Title"], FontColor -> GrayLevel[0.3]], Cell[ StyleData["Subsubtitle"]], Cell[ StyleData["Subsubtitle"]], Cell[ StyleData["Section"], CellFrameColor -> RGBColor[0.266667, 0.403922, 0.631373], FontColor -> RGBColor[0.266667, 0.403922, 0.631373]], Cell[ StyleData["Subsection"], FontColor -> RGBColor[0.266667, 0.403922, 0.631373]], Cell[ StyleData["Subsubsection"], CellMargins -> {{48, Inherited}, {Inherited, Inherited}}, FontFamily -> "Gill Sans MT", FontSize -> 16, FontWeight -> "Plain", FontSlant -> "Plain", FontVariations -> { "StrikeThrough" -> False, "Underline" -> False}, FontColor -> RGBColor[0.266667, 0.403922, 0.631373]], Cell[ StyleData["Input"], Background -> GrayLevel[1]], Cell[ StyleData["Output"], Background -> GrayLevel[1]], Cell[ StyleData["Graphics"], MenuSortingValue -> 10000, Background -> GrayLevel[1]], Cell[ StyleData["Figure"], CellFrame -> {{1, 1}, {0, 1}}, CellMargins -> {{60, 10}, {0, 0}}, MenuSortingValue -> 10000, Background -> RGBColor[1., 1., 1]], Cell[ StyleData["FigureCaption"], CellFrame -> {{1, 1}, {1, 0}}, CellMargins -> {{60, 10}, {0, 0}}, MenuSortingValue -> 10000, FontFamily -> "Arial Narrow", FontSize -> 14, FontWeight -> "Plain", FontSlant -> "Plain", FontVariations -> { "StrikeThrough" -> False, "Underline" -> False}, Background -> RGBColor[1, 1, 1]], Cell[ CellGroupData[{ Cell[ StyleData["Text"]], Cell[ StyleData["Notebook", "Printout"], PageFooters -> {{"", "", "copyright \[Copyright] N T Gladd 2016"}, { "", "", "copyright \[Copyright] N T Gladd 2016"}}, PrintingOptions -> {"FirstPageHeader" -> True}]}, Open]]}, Visible -> False, FrontEndVersion -> "11.0 for Microsoft Windows (64-bit) (July 28, 2016)", StyleDefinitions -> "PrivateStylesheetFormatting.nb"] ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 1110, 14, 243, "Title"], Cell[1692, 38, 99, 1, 54, "Subtitle"], Cell[1794, 41, 977, 21, 131, "Text"], Cell[2774, 64, 1287, 28, 205, "Input", InitializationCell->True], Cell[CellGroupData[{ Cell[4086, 96, 212, 3, 109, "Section"], Cell[4301, 101, 250, 5, 74, "Text"], Cell[4554, 108, 107, 1, 45, "Text"], Cell[CellGroupData[{ Cell[4686, 113, 758, 22, 59, "Input"], Cell[5447, 137, 785, 26, 61, "Output"] }, Open ]], Cell[6247, 166, 139, 1, 45, "Text"], Cell[CellGroupData[{ Cell[6411, 171, 267, 6, 59, "Input"], Cell[6681, 179, 738, 25, 189, "Output"] }, Open ]], Cell[7434, 207, 122, 1, 142, "Text"], Cell[CellGroupData[{ Cell[7581, 212, 1448, 42, 585, "Input"], Cell[9032, 256, 710, 22, 196, "Output"] }, Open ]], Cell[9757, 281, 122, 1, 142, "Text"], Cell[CellGroupData[{ Cell[9904, 286, 391, 11, 189, "Input"], Cell[10298, 299, 867, 25, 405, "Output"] }, Open ]], Cell[11180, 327, 280, 6, 142, "Text"], Cell[CellGroupData[{ Cell[11485, 337, 784, 23, 415, "Input"], Cell[12272, 362, 771, 24, 405, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13080, 391, 619, 17, 210, "Input"], Cell[13702, 410, 688, 20, 288, "Output"] }, Open ]], Cell[14405, 433, 127, 1, 142, "Text"], Cell[CellGroupData[{ Cell[14557, 438, 475, 13, 210, "Input"], Cell[15035, 453, 614, 18, 301, "Output"] }, Open ]], Cell[15664, 474, 170, 2, 142, "Text"], Cell[CellGroupData[{ Cell[15859, 480, 834, 23, 301, "Input"], Cell[16696, 505, 496, 14, 301, "Output"] }, Open ]], Cell[17207, 522, 144, 1, 142, "Text"], Cell[17354, 525, 173, 2, 142, "Text"], Cell[CellGroupData[{ Cell[17552, 531, 875, 24, 483, "Input"], Cell[18430, 557, 632, 19, 351, "Output"] }, Open ]], Cell[19077, 579, 303, 5, 142, "Text"] }, Open ]] }, Open ]] } ] *)