(* Content-type: application/vnd.wolfram.cdf.text *) (*** Wolfram CDF File ***) (* http://www.wolfram.com/cdf *) (* CreatedBy='Mathematica 10.1' *) (*************************************************************************) (* *) (* The Mathematica License under which this file was created prohibits *) (* restricting third parties in receipt of this file from republishing *) (* or redistributing it by any means, including but not limited to *) (* rights management or terms of use, without the express consent of *) (* Wolfram Research, Inc. For additional information concerning CDF *) (* licensing and redistribution see: *) (* *) (* www.wolfram.com/cdf/adopting-cdf/licensing-options.html *) (* *) (*************************************************************************) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 1064, 20] NotebookDataLength[ 22497, 615] NotebookOptionsPosition[ 23054, 610] NotebookOutlinePosition[ 23408, 626] CellTagsIndexPosition[ 23365, 623] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ StyleBox["Manipulate", FontSize->12], StyleBox["[", FontSize->12], RowBox[{ RowBox[{ StyleBox[ RowBox[{"Tmax", ":=", "100"}], FontSize->12], StyleBox[";", FontSize->12], StyleBox["\[IndentingNewLine]", FontSize->12], RowBox[{ StyleBox["Column", FontSize->12], StyleBox["[", FontSize->12], RowBox[{ RowBox[{ StyleBox["{", FontSize->12], StyleBox["\[IndentingNewLine]", FontSize->12], RowBox[{ StyleBox[ RowBox[{"Show", "[", RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "+", RowBox[{"r", "*", "x", "*", RowBox[{"(", RowBox[{"1", "-", "x"}], ")"}]}]}], ",", "x"}], "}"}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "1.6"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"Thick", ",", "Black"}], "}"}], ",", RowBox[{"{", RowBox[{"Dashed", ",", "Thick", ",", "Black"}], "}"}]}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1.6"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1.5"}], "}"}]}], "}"}]}], ",", RowBox[{"Frame", "\[Rule]", RowBox[{"{", RowBox[{"True", ",", "True", ",", "False", ",", "False"}], "}"}]}], ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"Style", "[", RowBox[{"\"\\"", ",", "18", ",", RowBox[{"FontFamily", "\[Rule]", "\"\\""}]}], "]"}], ",", RowBox[{"Style", "[", RowBox[{"\"\\"", ",", "18", ",", RowBox[{"FontFamily", "\[Rule]", "\"\\""}]}], "]"}]}], "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Graphics", "[", RowBox[{"Flatten", "[", RowBox[{"MapIndexed", "[", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Hue", "[", RowBox[{"0.8", "-", RowBox[{ RowBox[{"(", RowBox[{"0.8", "*", RowBox[{"#2", "[", RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "/", "Tmax"}]}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"#", ",", RowBox[{"If", "[", RowBox[{ RowBox[{"#", "\[Equal]", "x0"}], ",", "0", ",", "#"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"#", ",", RowBox[{"#", "+", RowBox[{"r", "#", " ", RowBox[{"(", RowBox[{"1", "-", "#"}], ")"}]}]}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"#", ",", " ", RowBox[{"#", "+", RowBox[{"r", "#", " ", RowBox[{"(", RowBox[{"1", "-", "#"}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"#", "+", RowBox[{"r", "#", " ", RowBox[{"(", RowBox[{"1", "-", "#"}], ")"}]}]}], ",", RowBox[{"#", "+", RowBox[{"r", "#", " ", RowBox[{"(", RowBox[{"1", "-", "#"}], ")"}]}]}]}], "}"}]}], "}"}], "]"}]}], "}"}], " ", "&"}], " ", ",", RowBox[{"NestWhileList", "[", RowBox[{ RowBox[{ RowBox[{"#", "+", RowBox[{"r", "#", " ", RowBox[{"(", RowBox[{"1", "-", "#"}], ")"}]}]}], "&"}], ",", " ", "x0", ",", " ", RowBox[{ RowBox[{"#1", " ", "\[NotEqual]", " ", RowBox[{"#", "+", RowBox[{"r", "#", " ", RowBox[{"(", RowBox[{"1", "-", "#"}], ")"}]}]}]}], " ", "&"}], ",", " ", "2", ",", " ", RowBox[{"Tnow", "-", "1"}]}], "]"}]}], "]"}], "]"}], "]"}], ",", RowBox[{"ImageSize", "\[Rule]", RowBox[{"{", RowBox[{"500", ",", "350"}], "}"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", RowBox[{"3.5", "/", "5"}]}]}], "\[IndentingNewLine]", "]"}], FontSize->12], StyleBox[",", FontSize->12], StyleBox["\[IndentingNewLine]", FontSize->12], RowBox[{ StyleBox[ RowBox[{"t", "=", RowBox[{"Range", "[", RowBox[{"1", ",", "Tnow"}], "]"}]}], FontSize->12], StyleBox[";", FontSize->12], StyleBox["\[IndentingNewLine]", FontSize->12], StyleBox[ RowBox[{"NN", "=", RowBox[{"NestList", "[", RowBox[{ RowBox[{ RowBox[{"#", "+", RowBox[{"r", " ", "#", " ", RowBox[{"(", RowBox[{"1", "-", "#"}], ")"}]}]}], "&"}], ",", "x0", ",", "Tnow"}], "]"}]}], FontSize->12], StyleBox[";", FontSize->12], StyleBox["\[IndentingNewLine]", FontSize->12], RowBox[{ StyleBox["Show", FontSize->12], StyleBox["[", FontSize->12], RowBox[{ StyleBox["{", FontSize->12], RowBox[{ RowBox[{ StyleBox["ListPlot", FontSize->12], StyleBox["[", FontSize->12], RowBox[{ StyleBox["NN", FontSize->12], StyleBox[",", FontSize->12], StyleBox["\[IndentingNewLine]", FontSize->12], StyleBox[" ", FontSize->12], StyleBox[ RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"Max", "[", RowBox[{"20", ",", "Tnow"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1.5"}], "}"}]}], "}"}]}], FontSize->12], StyleBox[",", FontSize->12], StyleBox["\[IndentingNewLine]", FontSize->12], StyleBox[" ", FontSize->12], StyleBox[ RowBox[{"Joined", "\[Rule]", "True"}], FontSize->12], StyleBox[",", FontSize->12], StyleBox[ RowBox[{"PlotStyle", "\[Rule]", "Thick"}], FontSize->12], StyleBox[",", FontSize->12], StyleBox[ RowBox[{"Frame", "\[Rule]", RowBox[{"{", RowBox[{"True", ",", "True", ",", "False", ",", "False"}], "}"}]}], FontSize->12], StyleBox[",", FontSize->12], StyleBox[ RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"Style", "[", RowBox[{"\"\