(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 10.1' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 21263, 569] NotebookOptionsPosition[ 20187, 531] NotebookOutlinePosition[ 20531, 546] CellTagsIndexPosition[ 20488, 543] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ RowBox[{ RowBox[{"h", "=", "0.01"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"c", "=", "1"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Tau]", "=", FractionBox["h", RowBox[{"2", " ", "c"}]]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"T", "=", "3"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"n", "=", RowBox[{"Ceiling", "[", RowBox[{"10", "/", "h"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"m", "=", RowBox[{"Ceiling", "[", RowBox[{"T", "/", "\[Tau]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"u0", "[", "x_", "]"}], ":=", TagBox[GridBox[{ {"\[Piecewise]", GridBox[{ {"1", RowBox[{ RowBox[{"1", "/", "4"}], "\[LessEqual]", "x", "\[LessEqual]", RowBox[{"3", "/", "4"}]}]}, { RowBox[{"1", "-", RowBox[{"Abs", "[", RowBox[{ RowBox[{"4", "x"}], "-", "6"}], "]"}]}], RowBox[{ RowBox[{"5", "/", "4"}], "\[LessEqual]", "x", "\[LessEqual]", RowBox[{"7", "/", "4"}]}]}, { RowBox[{"Cos", "[", RowBox[{"Pi", RowBox[{"(", RowBox[{ RowBox[{"2", "x"}], "-", "5"}], ")"}]}], "]"}], RowBox[{ RowBox[{"9", "/", "4"}], "\[LessEqual]", " ", "x", "\[LessEqual]", " ", RowBox[{"11", "/", "4"}]}]}, {"0", "True"} }, AllowedDimensions->{2, Automatic}, Editable->True, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxItemSize->{ "Columns" -> {{Automatic}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.84]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}, Selectable->True]} }, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxItemSize->{ "Columns" -> {{Automatic}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.35]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "Piecewise", DeleteWithContents->True, Editable->False, SelectWithContents->True, Selectable->False]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", "+", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"m", "+", "1"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"i", " ", "h"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", "\[LessEqual]", RowBox[{"n", "+", "1"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "=", RowBox[{"u0", "[", RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", " ", RowBox[{"j", "\[LessEqual]", "m"}], ",", " ", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "2"}], ",", RowBox[{"i", "\[LessEqual]", "n"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", ",", RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "=", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1", "-", SuperscriptBox[ RowBox[{"(", FractionBox[ RowBox[{"c", " ", "\[Tau]"}], "h"], ")"}], "2"]}], ")"}], RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ FractionBox["1", "2"], SuperscriptBox[ RowBox[{"(", FractionBox[ RowBox[{"c", " ", "\[Tau]"}], "h"], ")"}], "2"]}], "+", RowBox[{ FractionBox["1", "2"], RowBox[{"(", FractionBox[ RowBox[{"c", " ", "\[Tau]"}], "h"], ")"}]}]}], ")"}], RowBox[{"y", "[", RowBox[{"[", RowBox[{ RowBox[{"i", "-", "1"}], ",", "j"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ FractionBox["1", "2"], SuperscriptBox[ RowBox[{"(", FractionBox[ RowBox[{"c", " ", "\[Tau]"}], "h"], ")"}], "2"]}], "-", RowBox[{ FractionBox["1", "2"], RowBox[{"(", FractionBox[ RowBox[{"c", " ", "\[Tau]"}], "h"], ")"}]}]}], ")"}], RowBox[{"y", "[", RowBox[{"[", RowBox[{ RowBox[{"i", "+", "1"}], ",", "j"}], "]"}], "]"}]}]}]}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"1", ",", RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "=", "0"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{ RowBox[{"n", "+", "1"}], ",", RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "=", "0"}], ";"}]}], "\[IndentingNewLine]", "]"}]}], "Input", CellChangeTimes->{{3.6729945350826855`*^9, 3.6729946326762676`*^9}, { 3.6729947141439276`*^9, 3.6729948038280573`*^9}, {3.6729948894419537`*^9, 3.6729949365006456`*^9}, {3.672996709863076*^9, 3.672996738751728*^9}, { 3.7048619439003305`*^9, 3.7048620543004875`*^9}, {3.7048622664197874`*^9, 3.7048622870438166`*^9}, {3.704869258709442*^9, 3.704869297079497*^9}, { 3.7048693306395445`*^9, 3.7048693832096186`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Manipulate", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"plot1", "=", RowBox[{"ListLinePlot", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"n", "+", "1"}]}], "}"}]}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "6"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.5"}], ",", "1.5"}], "}"}]}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"plot2", "=", RowBox[{"Plot", "[", RowBox[{ RowBox[{"u0", "[", RowBox[{"x", "-", RowBox[{"c", " ", "j", " ", "\[Tau]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "10"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Red"}], ",", RowBox[{"ExclusionsStyle", "\[Rule]", "Red"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot1", ",", "plot2"}], "]"}]}], "\[IndentingNewLine]", ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"m", "+", "1"}], ",", "1"}], "}"}]}], "\[IndentingNewLine]", "]"}]], "Input", CellChangeTimes->{{3.6729948080743*^9, 3.672994872975012*^9}, { 3.704862068401508*^9, 3.7048620973125486`*^9}, {3.7048621415756116`*^9, 3.7048622447297573`*^9}, 3.7048623027378387`*^9, {3.704862335840886*^9, 3.704862340830893*^9}, {3.7048623826629524`*^9, 3.7048623989839754`*^9}, { 3.7048695206398163`*^9, 3.7048695206798162`*^9}, {3.7048696899460573`*^9, 3.7048696906160583`*^9}}], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`j$$ = 601, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`j$$], 1, 601, 1}}, Typeset`size$$ = { 360., {108., 111.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`j$8800$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`j$$ = 1}, "ControllerVariables" :> { Hold[$CellContext`j$$, $CellContext`j$8800$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> ($CellContext`plot1 = ListLinePlot[ Table[{ Part[$CellContext`x, $CellContext`i], Part[$CellContext`y, $CellContext`i, $CellContext`j$$]}, \ {$CellContext`i, 1, $CellContext`n + 1}], PlotRange -> {{0, 6}, {-0.5, 1.5}}]; $CellContext`plot2 = Plot[ $CellContext`u0[$CellContext`x - $CellContext`c $CellContext`j$$ \ $CellContext`\[Tau]], {$CellContext`x, 0, 10}, PlotStyle -> Red, ExclusionsStyle -> Red]; Show[$CellContext`plot1, $CellContext`plot2]), "Specifications" :> {{$CellContext`j$$, 1, 601, 1}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{409., {166., 171.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{ 3.704861849199904*^9, {3.704862049990481*^9, 3.70486209797255*^9}, 3.704862141955612*^9, {3.7048622039856997`*^9, 3.704862245199758*^9}, { 3.7048622918338237`*^9, 3.7048623034878397`*^9}, {3.704862389474962*^9, 3.704862399893977*^9}, {3.7048692207573886`*^9, 3.7048692326494055`*^9}, 3.7048693927396326`*^9, 3.7048695214198174`*^9, {3.7048696912460585`*^9, 3.704869707926082*^9}, 3.704869761266158*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`j$$ = 100, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`j$$], 1, 101, 1}}, Typeset`size$$ = { 360., {103., 107.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`j$969$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`j$$ = 1}, "ControllerVariables" :> { Hold[$CellContext`j$$, $CellContext`j$969$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> ListPlot[ Table[{ Part[$CellContext`x, $CellContext`i], Part[$CellContext`y, $CellContext`i, $CellContext`j$$]}, \ {$CellContext`i, 1, $CellContext`n + 1}]], "Specifications" :> {{$CellContext`j$$, 1, 101, 1}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{409., {149., 154.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Input"], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Table", "::", "iterb"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Iterator \[NoBreak]\\!\\({i, 1, \\(\\(1 + n\\)\\)}\\)\ \[NoBreak] does not have appropriate bounds. \\!\\(\\*ButtonBox[\\\"\ \[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/iterb\\\", ButtonNote -> \ \\\"Table::iterb\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.704861830610876*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Table", "::", "iterb"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Iterator \[NoBreak]\\!\\({i, 1, \\(\\(1 + n\\)\\)}\\)\ \[NoBreak] does not have appropriate bounds. \\!\\(\\*ButtonBox[\\\"\ \[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/iterb\\\", ButtonNote -> \ \\\"Table::iterb\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.7048618306308765`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Table", "::", "iterb"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Iterator \[NoBreak]\\!\\({i, 1.`16., \\(\\(1.`16. + \ n\\)\\)}\\)\[NoBreak] does not have appropriate bounds. \\!\\(\\*ButtonBox[\\\ \"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/iterb\\\", ButtonNote -> \ \\\"Table::iterb\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.7048618313688784`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"General", "::", "stop"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(Table :: \ iterb\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \ calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \ \\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.7048618314288783`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ListPlot", "::", "lpn"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\(Table[\\(\\(\\(\\({\\(\\(x \ \[LeftDoubleBracket] i \[RightDoubleBracket]\\)\\), \\(\\(y \ \[LeftDoubleBracket] \\(\\(i, FE`j$$16\\)\\) \ \[RightDoubleBracket]\\)\\)}\\)\\), \\(\\({i, 1.`, \\(\\(\\(\\(1.` \ \[VeryThinSpace]\\)\\) + n\\)\\)}\\)\\)\\)\\)]\\)\[NoBreak] is not a list of \ numbers or pairs of numbers. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/ListPlot\\\", ButtonNote -> \\\"ListPlot::lpn\\\"]\ \\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.704861831648879*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ListPlot", "::", "lpn"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\(Table[\\(\\(\\(\\({\\(\\(x \ \[LeftDoubleBracket] i \[RightDoubleBracket]\\)\\), \\(\\(y \ \[LeftDoubleBracket] \\(\\(i, FE`j$$16\\)\\) \ \[RightDoubleBracket]\\)\\)}\\)\\), \\(\\({i, 1.`, \\(\\(\\(\\(1.` \ \[VeryThinSpace]\\)\\) + n\\)\\)}\\)\\)\\)\\)]\\)\[NoBreak] is not a list of \ numbers or pairs of numbers. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/ListPlot\\\", ButtonNote -> \\\"ListPlot::lpn\\\"]\ \\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.7048618317488785`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Part", "::", "partw"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Part \[NoBreak]\\!\\(100\\)\[NoBreak] of \ \[NoBreak]\\!\\({0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \ 0, 0, 0, 0, \\(\\(\[LeftSkeleton] 11 \[RightSkeleton]\\)\\)}\\)\[NoBreak] \ does not exist. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/partw\\\", ButtonNote -> \ \\\"Part::partw\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.704862046330476*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Part", "::", "partw"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Part \[NoBreak]\\!\\(100\\)\[NoBreak] of \ \[NoBreak]\\!\\({0, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, \ 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, \ 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, \ 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, \\(\\(\[LeftSkeleton] 11 \ \[RightSkeleton]\\)\\)}\\)\[NoBreak] does not exist. \\!\\(\\*ButtonBox[\\\"\ \[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/partw\\\", ButtonNote -> \ \\\"Part::partw\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.704862046370476*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Part", "::", "partw"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Part \[NoBreak]\\!\\(100\\)\[NoBreak] of \ \[NoBreak]\\!\\({0, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, \ 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, \ 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, \ 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, 0.`, \\(\\(\[LeftSkeleton] 11 \ \[RightSkeleton]\\)\\)}\\)\[NoBreak] does not exist. \\!\\(\\*ButtonBox[\\\"\ \[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/partw\\\", ButtonNote -> \ \\\"Part::partw\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.7048620464004755`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"General", "::", "stop"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(Part :: \ partw\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \ calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \ \\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.7048620464304757`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"y", "[", RowBox[{"[", "2", "]"}], "]"}]], "Input", CellChangeTimes->{{3.672994882975584*^9, 3.67299488395264*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"u0", "[", "0.01`", "]"}], ",", RowBox[{ RowBox[{"5000.`", " ", RowBox[{"u0", "[", "0.`", "]"}]}], "-", RowBox[{"4999.`", " ", RowBox[{"u0", "[", "0.01`", "]"}]}]}]}], "}"}]], "Output", CellChangeTimes->{3.672994884215655*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData["\[Tau]"], "Input", CellChangeTimes->{{3.672999850566706*^9, 3.67299985080472*^9}}], Cell[BoxData["0.05`"], "Output", CellChangeTimes->{3.6729998510257325`*^9}] }, Open ]] }, WindowSize->{1904, 998}, WindowMargins->{{-8, Automatic}, {Automatic, 0}}, FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (November 7, 2010)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[558, 20, 6693, 195, 502, "Input"], Cell[CellGroupData[{ Cell[7276, 219, 1980, 51, 152, "Input"], Cell[9259, 272, 2408, 47, 354, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11704, 324, 1620, 35, 320, "Input"], Cell[13327, 361, 478, 10, 21, "Message"], Cell[13808, 373, 480, 10, 21, "Message"], Cell[14291, 385, 490, 10, 21, "Message"], Cell[14784, 397, 531, 11, 21, "Message"], Cell[15318, 410, 698, 14, 21, "Message"], Cell[16019, 426, 700, 14, 21, "Message"], Cell[16722, 442, 670, 13, 21, "Message"], Cell[17395, 457, 770, 14, 21, "Message"], Cell[18168, 473, 772, 14, 21, "Message"], Cell[18943, 489, 530, 11, 21, "Message"] }, Open ]], Cell[CellGroupData[{ Cell[19510, 505, 142, 3, 31, "Input"], Cell[19655, 510, 303, 9, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[19995, 524, 97, 1, 31, "Input"], Cell[20095, 527, 76, 1, 30, "Output"] }, Open ]] } ] *) (* End of internal cache information *)