(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 11.2' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 59853, 1380] NotebookOptionsPosition[ 57780, 1342] NotebookOutlinePosition[ 58126, 1357] CellTagsIndexPosition[ 58083, 1354] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"Quit", "[", "]"}], ";"}]], "Input", CellChangeTimes->{{3.8579466011214185`*^9, 3.8579466033732166`*^9}},ExpressionUUID->"9b358bc8-61b5-45f2-98ae-\ a2f0f3b6e844"], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "para", " ", "la", " ", "salida", " ", "de", " ", "gr\[AAcute]ficos"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"antialias", "[", RowBox[{"g_", ",", RowBox[{"n_:", " ", "1"}]}], "]"}], ":=", RowBox[{"ImageResize", "[", RowBox[{ RowBox[{"Rasterize", "[", RowBox[{"g", ",", "\"\\"", ",", RowBox[{"ImageResolution", "\[Rule]", RowBox[{"n", " ", "72"}]}], ",", " ", RowBox[{"Background", "\[Rule]", "None"}]}], "]"}], ",", RowBox[{"Scaled", "[", RowBox[{"1", "/", "n"}], "]"}]}], "]"}]}], ";"}], " "}]], "Input", CellChangeTimes->{{3.8580846283952775`*^9, 3.8580846503437457`*^9}},ExpressionUUID->"ef9c96d3-dc19-445e-a48a-\ bb9c27e27e0a"], Cell[CellGroupData[{ Cell[TextData[{ "P\[EAcute]ndulo doble: \[Lambda] = l1/l2, \[Mu] = ", Cell[BoxData[ FormBox[ FractionBox["m2", RowBox[{"m1", " ", "+", " ", "m2"}]], TraditionalForm]], FormatType->"TraditionalForm",ExpressionUUID-> "717194e7-71e6-4d85-be6a-0184a588d3ce"], "." }], "Section", CellChangeTimes->{{3.858084547081337*^9, 3.858084567483689*^9}},ExpressionUUID->"3e90e1f8-09cb-4dac-a915-\ db3d64c0c680"], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"par\[AAcute]metros", " ", "del", " ", RowBox[{"p\[EAcute]ndulo", ".", " ", "\[Mu]"}], " ", "tiene", " ", "que", " ", "ser", " ", "menor", " ", "que", " ", "1"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Lambda]", " ", "=", " ", RowBox[{"1", "/", "2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Mu]", " ", "=", " ", RowBox[{"2", "/", "3"}]}], ";"}], " ", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"condiciones", " ", "iniciales"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Phi]10", " ", "=", " ", "3"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Phi]20", " ", "=", " ", "2"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Phi]1p0", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Phi]2p0", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "Este", " ", "valor", " ", "es", " ", "para", " ", "integrar", " ", "las", " ", "ecuaciones", " ", "un", " ", "poco", " ", "en", " ", "el", " ", "pasado"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tpast", " ", "=", " ", RowBox[{"-", "10"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "Integrador", " ", "num\[EAcute]rico", " ", "de", " ", "las", " ", "ecuaciones"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"sol", " ", "=", " ", RowBox[{ RowBox[{"NDSolve", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"ecuaciones", " ", "diferenciales"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Lambda]", " ", RowBox[{ RowBox[{"\[Phi]1", "''"}], "[", "t", "]"}]}], " ", "+", " ", RowBox[{"\[Mu]", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Phi]2", "''"}], "[", "t", "]"}], " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{"\[Phi]2", "[", "t", "]"}], " ", "-", " ", RowBox[{"\[Phi]1", "[", "t", "]"}]}], "]"}]}], " ", "-", " ", RowBox[{ SuperscriptBox[ RowBox[{ RowBox[{"\[Phi]2", "'"}], "[", "t", "]"}], "2"], RowBox[{"Sin", "[", RowBox[{ RowBox[{"\[Phi]2", "[", "t", "]"}], " ", "-", " ", RowBox[{"\[Phi]1", "[", "t", "]"}]}], "]"}]}]}], ")"}]}], " ", "+", RowBox[{"Sin", "[", RowBox[{"\[Phi]1", "[", "t", "]"}], "]"}]}], " ", "\[Equal]", " ", "0"}], ",", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Phi]2", "''"}], "[", "t", "]"}], " ", "+", " ", RowBox[{"\[Lambda]", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Phi]1", "''"}], "[", "t", "]"}], " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{"\[Phi]2", "[", "t", "]"}], " ", "-", " ", RowBox[{"\[Phi]1", "[", "t", "]"}]}], "]"}]}], " ", "+", " ", RowBox[{ SuperscriptBox[ RowBox[{ RowBox[{"\[Phi]1", "'"}], "[", "t", "]"}], "2"], RowBox[{"Sin", "[", RowBox[{ RowBox[{"\[Phi]2", "[", "t", "]"}], " ", "-", " ", RowBox[{"\[Phi]1", "[", "t", "]"}]}], "]"}]}]}], ")"}]}], " ", "+", RowBox[{"Sin", "[", RowBox[{"\[Phi]2", "[", "t", "]"}], "]"}]}], " ", "\[Equal]", " ", "0"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"condiciones", " ", "iniciales"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Phi]1", "[", "0", "]"}], " ", "\[Equal]", " ", "\[Phi]10"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Phi]2", "[", "0", "]"}], " ", "\[Equal]", " ", "\[Phi]20"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Phi]1", "'"}], "[", "0", "]"}], " ", "\[Equal]", "\[Phi]1p0"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Phi]2", "'"}], "[", "0", "]"}], " ", "\[Equal]", " ", "\[Phi]2p0"}]}], "\[IndentingNewLine]", "}"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"\[Phi]1", ",", " ", "\[Phi]2"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "tpast", ",", " ", RowBox[{"tmax", " ", "=", " ", "50"}]}], "}"}]}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "Se", " ", "puede", " ", "aumentar", " ", "la", " ", "precisi\[OAcute]n", " ", "del", " ", "c\[AAcute]lculo", " ", "descomentando", " ", "la", " ", "siguiente", " ", RowBox[{"l\[IAcute]nea", ".", " ", "Ajusten"}], " ", "los", " ", "valores", " ", "como", " ", "mejor", " ", "les", " ", RowBox[{"convenga", ".", " ", "Esto"}], " ", "se", " ", "nota", " ", "si", " ", "grafican", " ", "la", " ", "conservaci\[OAcute]n", " ", "de", " ", "la", " ", "energ\[IAcute]a"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"WorkingPrecision", "\[Rule]", " ", "20"}], ",", " ", RowBox[{"MaxSteps", "\[Rule]", " ", SuperscriptBox["10", "6"]}], ","}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "Extrae", " ", "cada", " ", "\[AAcute]ngulo", " ", "de", " ", "la", " ", "soluci\[OAcute]n"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"\[Phi]1sol", ",", " ", "\[Phi]2sol"}], "}"}], " ", "=", " ", RowBox[{"sol", "[", RowBox[{"[", RowBox[{"All", ",", "2"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "contruye", " ", "las", " ", "funciones", " ", "posici\[OAcute]n", " ", "de", " ", "cada", " ", "masa"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"r1sol", "[", "t_", "]"}], ":=", " ", RowBox[{"\[Lambda]", " ", RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", RowBox[{"\[Phi]1sol", "[", "t", "]"}], "]"}], ",", " ", RowBox[{"-", RowBox[{"Cos", "[", RowBox[{"\[Phi]1sol", "[", "t", "]"}], "]"}]}]}], "}"}]}]}], " ", ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"r2sol", "[", "t_", "]"}], ":=", " ", RowBox[{ RowBox[{"r1sol", "[", "t", "]"}], " ", "+", " ", RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", RowBox[{"\[Phi]2sol", "[", "t", "]"}], "]"}], ",", " ", RowBox[{"-", RowBox[{"Cos", "[", RowBox[{"\[Phi]2sol", "[", "t", "]"}], "]"}]}]}], "}"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"energ\[IAcute]a", " ", "mec\[AAcute]nica"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"e", "[", "t_", "]"}], ":=", " ", RowBox[{ RowBox[{ FractionBox["1", "2"], SuperscriptBox["\[Lambda]", "2"], SuperscriptBox[ RowBox[{ RowBox[{"\[Phi]1sol", "'"}], "[", "t", "]"}], "2"]}], " ", "+", " ", RowBox[{ FractionBox["1", "2"], "\[Mu]", RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{ RowBox[{"\[Phi]2sol", "'"}], "[", "t", "]"}], "2"], "+", " ", RowBox[{"2", " ", "\[Lambda]", " ", RowBox[{ RowBox[{"\[Phi]1sol", "'"}], "[", "t", "]"}], " ", RowBox[{ RowBox[{"\[Phi]2sol", "'"}], "[", "t", "]"}], " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{"\[Phi]2sol", "[", "t", "]"}], " ", "-", " ", RowBox[{"\[Phi]1sol", "[", "t", "]"}]}], "]"}]}]}], ")"}]}], " ", "-", " ", RowBox[{"\[Lambda]", " ", RowBox[{"Cos", "[", RowBox[{"\[Phi]1sol", "[", "t", "]"}], "]"}]}], " ", "-", " ", RowBox[{"\[Mu]", " ", RowBox[{"Cos", "[", RowBox[{"\[Phi]2sol", "[", "t", "]"}], "]"}]}]}]}]}]}]], "Input", CellChangeTimes->{{3.8579034695664053`*^9, 3.857903478505571*^9}, { 3.8579035185838118`*^9, 3.857903703945401*^9}, {3.857903741382861*^9, 3.8579039176066656`*^9}, {3.8579039707971287`*^9, 3.857904133691585*^9}, 3.857904980231278*^9, {3.8579050644653263`*^9, 3.8579050950901237`*^9}, 3.857905207088758*^9, {3.857905322898846*^9, 3.857905326227235*^9}, 3.857905398685835*^9, {3.85790546690773*^9, 3.8579054912531605`*^9}, 3.8579055438864384`*^9, {3.8579057786835265`*^9, 3.8579057787929096`*^9}, { 3.8579058184600525`*^9, 3.8579058203656635`*^9}, {3.8579098019411764`*^9, 3.8579098553962326`*^9}, {3.8579099289566627`*^9, 3.857909955985346*^9}, 3.8579100222158976`*^9, {3.857913256044149*^9, 3.857913256153489*^9}, 3.857913301146863*^9, {3.857942858973325*^9, 3.8579428658028316`*^9}, { 3.857942913212735*^9, 3.8579429382539062`*^9}, {3.8579430091014194`*^9, 3.857943009159853*^9}, {3.857943049092434*^9, 3.8579431032724667`*^9}, { 3.858003384657833*^9, 3.8580033914777822`*^9}, {3.8580034273676023`*^9, 3.858003436577983*^9}, {3.8580034793577967`*^9, 3.8580034816077785`*^9}, 3.858003759082328*^9, {3.8580051487277513`*^9, 3.858005148857129*^9}, { 3.858006132324269*^9, 3.8580061325275154`*^9}, {3.8580165431231537`*^9, 3.858016543212401*^9}, {3.8580180270040183`*^9, 3.8580180298228903`*^9}, { 3.858018257095784*^9, 3.8580182584626083`*^9}, {3.8580193696008997`*^9, 3.8580193755496144`*^9}, 3.8580194449914455`*^9, {3.8580194840906663`*^9, 3.8580194954038115`*^9}, {3.8580195641755695`*^9, 3.8580195677069936`*^9}, {3.8580196138550105`*^9, 3.858019627850872*^9}, { 3.85801987406936*^9, 3.8580198791551733`*^9}, {3.8580205186852875`*^9, 3.858020527438131*^9}, {3.858024417543952*^9, 3.858024418015189*^9}, { 3.8580846577710233`*^9, 3.8580848096602354`*^9}, {3.858084869058654*^9, 3.8580850000219564`*^9}, {3.858085616598116*^9, 3.8580857132761383`*^9}, { 3.8580857464131675`*^9, 3.8580857680969486`*^9}, {3.8580857986890783`*^9, 3.8580857988295164`*^9}, {3.8580927967839656`*^9, 3.858092796924554*^9}, { 3.858092899278448*^9, 3.8580929220457363`*^9}},ExpressionUUID->"8c33cc18-b350-403a-9e82-\ 7edf0bd84570"], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "un", " ", "gr\[AAcute]fico", " ", "de", " ", "la", " ", "posici\[OAcute]n", " ", "de", " ", "la", " ", "segunda", " ", "masa"}], "*)"}], "\[IndentingNewLine]", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"r2sol", "[", "tt", "]"}], ",", " ", RowBox[{"{", RowBox[{"tt", ",", " ", "0", ",", " ", "tmax"}], "}"}], ",", " ", RowBox[{"ColorFunction", "\[Rule]", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"x", ",", " ", "y", ",", " ", "u"}], "}"}], ",", RowBox[{"Directive", "[", RowBox[{"Opacity", "[", SqrtBox["u"], "]"}], "]"}]}], "]"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.857942883992222*^9, 3.85794289694549*^9}, { 3.8579430208959637`*^9, 3.8579430213802166`*^9}, {3.85800341035807*^9, 3.8580034133980837`*^9}, {3.858003447607875*^9, 3.8580034515275807`*^9}, { 3.8580158024184103`*^9, 3.8580158095952573`*^9}, {3.8580848159870596`*^9, 3.8580848262048197`*^9}},ExpressionUUID->"20627f30-a074-4ad5-a85b-\ de7b1268e934"], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ "verificaci\[OAcute]n", " ", "de", " ", "que", " ", "la", " ", "energ\[IAcute]a", " ", "se", " ", "conserva"}], " ", "-", " ", "bueno"}], ",", " ", RowBox[{"m\[AAcute]s", " ", "o", " ", RowBox[{"menos", ".", " ", "Descomenten"}], " ", "las", " ", "l\[IAcute]neas", " ", "del", " ", "WorkingPresicion", " ", "para", " ", "mejorar", " ", "la", " ", RowBox[{"precisi\[OAcute]n", "."}]}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"e0", " ", "=", " ", RowBox[{"e", "[", "0", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{ FractionBox[ RowBox[{"e", "[", "t", "]"}], "e0"], "-", "1"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "0", ",", " ", "tmax"}], "}"}], ",", " ", RowBox[{"PlotRange", "\[Rule]", " ", "All"}], ",", " ", RowBox[{"ImageSize", "\[Rule]", "600"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.8579040826300745`*^9, 3.857904100318185*^9}, { 3.857905337867498*^9, 3.857905375691428*^9}, {3.8580848484539547`*^9, 3.858084859763959*^9}, {3.8580850165963736`*^9, 3.8580850467774196`*^9}},ExpressionUUID->"a3bba777-776f-444c-837c-\ 8e5c07fa006c"], Cell[BoxData[ RowBox[{ RowBox[{"(*", "Animaci\[OAcute]n", "*)"}], "\[IndentingNewLine]", RowBox[{"DynamicModule", "[", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ "Para", " ", "que", " ", "no", " ", "se", " ", "mezclen", " ", "con", " ", "variables", " ", "din\[AAcute]micas", " ", "de", " ", "otros", " ", "Manipulate"}], ",", " ", RowBox[{ "las", " ", "siguientes", " ", "variables", " ", "se", " ", "definen", " ", "localmente", " ", "dentro", " ", "de", " ", "este", " ", "\"\\""}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{ "\[Phi]1pr", ",", " ", "\[Phi]2pr", ",", " ", "r0", ",", " ", "r1", ",", " ", "r2", ",", " ", "traza", ",", " ", "tpr"}], "}"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"r0", " ", "=", " ", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Manipulate", "[", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "La", " ", "variable", " ", "de", " ", "control", " ", "de", " ", "este", " ", "bloque", " ", "\"\\"", " ", "es", " ", "el", " ", "tiempo", " ", "tpr"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "Eval\[UAcute]a", " ", "los", " ", "\[AAcute]ngulos", " ", "y", " ", "las", " ", "posiciones", " ", "a", " ", "tiempo", " ", "tpr"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Phi]1pr", " ", "=", " ", RowBox[{"\[Phi]1sol", "[", "tpr", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"\[Phi]2pr", " ", "=", " ", RowBox[{"\[Phi]2sol", "[", "tpr", "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"r1", " ", "=", " ", RowBox[{"\[Lambda]", RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", "\[Phi]1pr", "]"}], ",", RowBox[{"-", RowBox[{"Cos", "[", "\[Phi]1pr", "]"}]}]}], " ", "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"r2", " ", "=", " ", RowBox[{"r1", " ", "+", " ", RowBox[{"{", " ", RowBox[{ RowBox[{"Sin", "[", "\[Phi]2pr", "]"}], ",", RowBox[{"-", RowBox[{"Cos", "[", "\[Phi]2pr", "]"}]}]}], "}"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ "Grafica", " ", "parte", " ", "de", " ", "la", " ", "trayectoria", " ", "de", " ", "la", " ", "segunda", " ", "masa"}], ",", " ", RowBox[{ "para", " ", "mostrar", " ", "la", " ", "traza", " ", "que", " ", "va", " ", "dejando"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"traza", " ", "=", " ", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"r2sol", "[", "tt", "]"}], ",", " ", RowBox[{"{", RowBox[{"tt", ",", " ", RowBox[{"tpr", "-", "4"}], ",", " ", "tpr"}], "}"}], ",", " ", RowBox[{"ColorFunction", "\[Rule]", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"x", ",", " ", "y", ",", " ", "u"}], "}"}], ",", RowBox[{"Directive", "[", RowBox[{"Opacity", "[", SqrtBox["u"], "]"}], "]"}]}], "]"}]}], ",", " ", RowBox[{"PlotPoints", "\[Rule]", " ", "50"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "\"\\"", " ", "re\[UAcute]ne", " ", "en", " ", "un", " ", "mismo", " ", "gr\[AAcute]fico", " ", "varios", " ", "elementos"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Show", "[", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "primero", " ", "muestra", " ", "la", " ", "traza", " ", "generada", " ", "m\[AAcute]s", " ", "arriba"}], "*)"}], "\[IndentingNewLine]", RowBox[{"traza", ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "despu\[EAcute]s", " ", "grafica", " ", "las", " ", "barras", " ", "y", " ", "las", " ", "masas"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Graphics", "[", "\[IndentingNewLine]", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"AbsoluteThickness", "[", "3", "]"}], ",", " ", "\[IndentingNewLine]", "Orange", ",", " ", RowBox[{"Line", "[", RowBox[{"{", RowBox[{"r0", ",", " ", "r1"}], "}"}], "]"}], ",", "\[IndentingNewLine]", "Green", ",", " ", RowBox[{"Line", "[", RowBox[{"{", RowBox[{"r1", ",", " ", "r2"}], "}"}], "]"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", "White", ",", " ", RowBox[{"Disk", "[", RowBox[{"r1", ",", " ", RowBox[{"radio", " ", "=", " ", "0.05"}]}], "]"}], ",", "\[IndentingNewLine]", "Red", ",", " ", RowBox[{"Disk", "[", RowBox[{"r0", ",", " ", RowBox[{ FractionBox["2", "3"], "radio"}]}], "]"}], ",", "\[IndentingNewLine]", "Cyan", ",", " ", RowBox[{"Disk", "[", RowBox[{"r2", ",", " ", "radio"}], "]"}]}], "\[IndentingNewLine]", "}"}], "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ "Se", " ", "da", " ", "espacio", " ", "suficiente", " ", "para", " ", "que", " ", "las", " ", "masas", " ", "no", " ", "se", " ", "salgan", " ", "de", " ", "la", " ", "zona", " ", RowBox[{"graficada", ".", " ", "En"}], " ", "estas", " ", "variables", " ", "adimensionales"}], ",", " ", RowBox[{ "la", " ", "longitud", " ", "de", " ", "la", " ", "primera", " ", "barra", " ", "es", " ", "\[Lambda]", " ", "y", " ", "la", " ", "de", " ", "la", " ", "segunda", " ", "es", " ", "1."}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"PlotRange", "\[Rule]", " ", RowBox[{"1.2", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "1"}], "-", "\[Lambda]"}], ",", " ", RowBox[{"1", "+", "\[Lambda]"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "1"}], "-", "\[Lambda]"}], ",", " ", RowBox[{"1", "+", "\[Lambda]"}]}], "}"}]}], "}"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"Axes", "\[Rule]", " ", "False"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Background", "\[Rule]", " ", "Black"}], ",", " ", "\[IndentingNewLine]", RowBox[{"ImageSize", "\[Rule]", " ", "500"}]}], "\[IndentingNewLine]", "]"}]}], ",", " ", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"Esto", " ", "controla", " ", "el", " ", "Manipulate"}], "*)"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{ "tpr", ",", " ", "0", ",", " ", "tmax", ",", " ", "0.01", ",", " ", RowBox[{"Appearance", "\[Rule]", " ", "\"\\""}], ",", " ", RowBox[{"AnimationRate", "\[Rule]", " ", "2"}]}], "}"}], ",", " ", "\[IndentingNewLine]", RowBox[{"TrackedSymbols", "\[RuleDelayed]", RowBox[{"{", "tpr", "}"}]}]}], "\[IndentingNewLine]", "]"}]}]}], "\[IndentingNewLine]", "]"}]}]], "Input", CellChangeTimes->{{3.857904270521346*^9, 3.8579042773010054`*^9}, { 3.857904311293007*^9, 3.857904497131488*^9}, {3.857904569366128*^9, 3.8579046586677847`*^9}, {3.8579046985975647`*^9, 3.85790477977835*^9}, { 3.8579051859843473`*^9, 3.857905222682843*^9}, {3.8579096794606433`*^9, 3.857909689239403*^9}, {3.8579097363537946`*^9, 3.857909797770282*^9}, { 3.857909877990615*^9, 3.857909924379796*^9}, {3.857909974546788*^9, 3.857909993361618*^9}, 3.8579100301361*^9, {3.8579100708450823`*^9, 3.8579101047234445`*^9}, {3.8579101358879976`*^9, 3.8579101371247044`*^9}, {3.8579103582979765`*^9, 3.8579104753965034`*^9}, {3.8579144807500305`*^9, 3.857914510346921*^9}, { 3.8579367960972204`*^9, 3.857936797552803*^9}, {3.8580035012682695`*^9, 3.8580035023778667`*^9}, {3.858003643777933*^9, 3.858003644488147*^9}, { 3.8580037120139627`*^9, 3.8580037188193865`*^9}, {3.858085054463415*^9, 3.858085292666889*^9}, {3.8580858184223566`*^9, 3.858085982800053*^9}, { 3.8580925285730124`*^9, 3.858092629268042*^9}},ExpressionUUID->"662f790d-04c3-4b02-8091-\ 44675ec4add0"], Cell[CellGroupData[{ Cell["\<\ Salida de gr\[AAcute]ficos: lo mismo de antes pero para construir una tabla \ con los gr\[AAcute]ficos rasterizados v\[IAcute]a la funci\[OAcute]n \ \[OpenCurlyDoubleQuote]antialias\[CloseCurlyDoubleQuote] definida al comienzo \ del notebook.\ \>", "Subsection", CellChangeTimes->{{3.8579130541853848`*^9, 3.8579130570128465`*^9}, { 3.8580853368183823`*^9, 3.858085381872067*^9}, {3.858092652893099*^9, 3.858092664734049*^9}},ExpressionUUID->"f6ad20cb-733c-40e5-a05f-\ 214a4b7b8fc4"], Cell[BoxData[ RowBox[{ RowBox[{"Monitor", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tgraf", " ", "=", " ", "\[IndentingNewLine]", RowBox[{"Table", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Phi]1pr", " ", "=", " ", RowBox[{"\[Phi]1sol", "[", "tpr", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"\[Phi]2pr", " ", "=", " ", RowBox[{"\[Phi]2sol", "[", "tpr", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"r0", " ", "=", " ", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"r1", " ", "=", " ", RowBox[{"\[Lambda]", RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", "\[Phi]1pr", "]"}], ",", RowBox[{"-", RowBox[{"Cos", "[", "\[Phi]1pr", "]"}]}]}], " ", "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"r2", " ", "=", " ", RowBox[{"r1", " ", "+", " ", RowBox[{"{", " ", RowBox[{ RowBox[{"Sin", "[", "\[Phi]2pr", "]"}], ",", RowBox[{"-", RowBox[{"Cos", "[", "\[Phi]2pr", "]"}]}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"traza", " ", "=", " ", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"r2sol", "[", "tt", "]"}], ",", " ", RowBox[{"{", RowBox[{"tt", ",", " ", RowBox[{"tpr", "-", "5"}], ",", " ", "tpr"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"AbsoluteThickness", "[", "2", "]"}]}], ",", " ", RowBox[{"ColorFunction", "\[Rule]", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"x", ",", " ", "y", ",", " ", "u"}], "}"}], ",", RowBox[{"Directive", "[", RowBox[{"Opacity", "[", SqrtBox["u"], "]"}], "]"}]}], "]"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"gtemp", " ", "=", " ", RowBox[{"antialias", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"Show", "[", "\[IndentingNewLine]", RowBox[{"traza", ",", "\[IndentingNewLine]", RowBox[{"Graphics", "[", "\[IndentingNewLine]", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"AbsoluteThickness", "[", "3", "]"}], ",", " ", "\[IndentingNewLine]", "Orange", ",", " ", RowBox[{"Line", "[", RowBox[{"{", RowBox[{"r0", ",", " ", "r1"}], "}"}], "]"}], ",", "\[IndentingNewLine]", "Green", ",", " ", RowBox[{"Line", "[", RowBox[{"{", RowBox[{"r1", ",", " ", "r2"}], "}"}], "]"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", "White", ",", " ", RowBox[{"Disk", "[", RowBox[{"r1", ",", " ", RowBox[{"radio", " ", "=", " ", "0.05"}]}], "]"}], ",", "\[IndentingNewLine]", "Red", ",", " ", RowBox[{"Disk", "[", RowBox[{"r0", ",", " ", RowBox[{ FractionBox["2", "3"], "radio"}]}], "]"}], ",", "\[IndentingNewLine]", "Cyan", ",", " ", RowBox[{"Disk", "[", RowBox[{"r2", ",", " ", "radio"}], "]"}]}], "\[IndentingNewLine]", "}"}], "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", RowBox[{"PlotRange", "\[Rule]", " ", RowBox[{"1.2", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "1"}], "-", "\[Lambda]"}], ",", " ", RowBox[{"1", "+", "\[Lambda]"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "1"}], "-", "\[Lambda]"}], ",", " ", RowBox[{"1", "+", "\[Lambda]"}]}], "}"}]}], "}"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"Axes", "\[Rule]", " ", "False"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Background", "\[Rule]", " ", "Black"}], ",", " ", "\[IndentingNewLine]", RowBox[{"ImageSize", "\[Rule]", " ", RowBox[{"(", RowBox[{"imsize", " ", "=", " ", "500"}], ")"}]}]}], "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", "2"}], " ", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ "este", " ", "n\[UAcute]mero", " ", "gobierna", " ", "el", " ", "rasterizado", " ", RowBox[{"(", "resoluci\[OAcute]n", ")"}]}], ";", " ", RowBox[{ "1", " ", "deja", " ", "el", " ", "gr\[AAcute]fico", " ", "tal", " ", "como", " ", "sale"}]}], ",", " ", RowBox[{"2", " ", "lo", " ", "suaviza", " ", "un", " ", "poco"}], ",", " ", RowBox[{"y", " ", "de", " ", "ah\[IAcute]", " ", "para", " ", RowBox[{"arriba", ".", " ", "Un"}], " ", "valor", " ", "muy", " ", "alto", " ", "hace", " ", "lenta", " ", "la", " ", "generaci\[OAcute]n", " ", "de", " ", "los", " ", "gr\[AAcute]ficos"}]}], "*)"}], "\[IndentingNewLine]", "]"}]}]}], ",", " ", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"Grad\[UAcute]en", " ", "el", " ", "paso", " ", RowBox[{"temporal", ".", " ", "Empiecen"}], " ", "con", " ", "un", " ", "\[CapitalDelta]t", " ", "grande", " ", "y", " ", "si", " ", "ven", " ", "que", " ", "todo", " ", "va", " ", "bien"}], ",", " ", "dismin\[UAcute]yanlo"}], "*)"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{"tpr", ",", " ", "0", ",", " ", "tmax", ",", " ", "0.1"}], "}"}]}], "]"}]}], " ", ";"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "Esta", " ", "l\[IAcute]nea", " ", "pertenece", " ", "al", " ", "bloque", " ", RowBox[{"\"\\"", ".", " ", "Va"}], " ", "mostrando", " ", "las", " ", "salidas", " ", "a", " ", "medida", " ", "que", " ", "se", " ", RowBox[{"calculan", ".", " ", "Sirve"}], " ", "para", " ", "ver", " ", "que", " ", "todo", " ", "marche", " ", RowBox[{"bien", "."}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Show", "[", RowBox[{"gtemp", ",", " ", RowBox[{"ImageSize", "\[Rule]", " ", "imsize"}]}], "]"}], ",", " ", "tpr"}], "}"}]}], "\[IndentingNewLine]", "]"}], ";"}]], "Input", CellChangeTimes->{{3.857911898122059*^9, 3.857911999654297*^9}, { 3.8579177423653507`*^9, 3.857917742515665*^9}, {3.858003866153293*^9, 3.858003917617735*^9}, {3.858006164367814*^9, 3.8580061647897687`*^9}, { 3.85801592135305*^9, 3.8580159604289813`*^9}, {3.8580165009234524`*^9, 3.858016579449478*^9}, {3.8580244963218174`*^9, 3.8580245261369047`*^9}, { 3.8580245566231146`*^9, 3.8580245602350206`*^9}, {3.8580316343070264`*^9, 3.8580316460819335`*^9}, {3.858031749594594*^9, 3.8580318321611924`*^9}, { 3.858032365484688*^9, 3.8580323940302925`*^9}, {3.858085357970574*^9, 3.8580854930016437`*^9}, {3.8580926876158657`*^9, 3.85809268781894*^9}, { 3.858092851932968*^9, 3.858092877940652*^9}, {3.858092985966375*^9, 3.8580930508180275`*^9}, {3.8580931095101247`*^9, 3.85809316840291*^9}, 3.8580932020398684`*^9},ExpressionUUID->"2411b3af-b027-482f-b721-\ 0a5c44b3d0ad"], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "Aqu\[IAcute]", " ", "pueden", " ", "animar", " ", "la", " ", "tabla", " ", "de", " ", "gr\[AAcute]ficos"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"tgraf", "[", RowBox[{"[", "i", "]"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", RowBox[{"Length", "@", "tgraf"}], ",", " ", "1", ",", " ", RowBox[{"Appearance", "\[Rule]", " ", "\"\\""}], ",", " ", RowBox[{"AnimationRate", "\[Rule]", " ", "25"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.857912106068201*^9, 3.8579121352211018`*^9}, { 3.858016070417719*^9, 3.8580160705329456`*^9}, {3.8580166586404333`*^9, 3.8580166588456287`*^9}, {3.858085498078576*^9, 3.858085526088711*^9}},ExpressionUUID->"146fa7ae-c10c-42ff-81ac-\ fd08b6480725"] }, Open ]], Cell[CellGroupData[{ Cell["Salida de archivos", "Subsection", CellChangeTimes->{{3.858085543069085*^9, 3.8580855457562833`*^9}},ExpressionUUID->"f747a3e6-cedc-4410-af73-\ e22310b3bec4"], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "Crea", " ", "y", " ", "selecciona", " ", "el", " ", "directorio", " ", "de", " ", "salida", " ", "dentro", " ", "del", " ", "directorio", " ", "en", " ", "donde", " ", "se", " ", "encuentra", " ", "guardado", " ", "el", " ", "notebook"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"SetDirectory", "[", RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}], "\n", RowBox[{ RowBox[{"CreateDirectory", "[", RowBox[{ "dirsal", " ", "=", " ", "\"\\""}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"SetDirectory", "[", "dirsal", "]"}]}]}]], "Input", CellChangeTimes->{{3.857913383121125*^9, 3.857913395370955*^9}, 3.858085575963005*^9, {3.8580927080562434`*^9, 3.858092744201603*^9}},ExpressionUUID->"7f711c5b-f110-477d-ace2-\ ed407a0389b2"], Cell[BoxData[ RowBox[{"Monitor", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", " ", "=", " ", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Length", "@", "tgraf"}]}], ",", " ", RowBox[{"++", "i"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"Export", "[", RowBox[{ RowBox[{"StringJoin", "[", RowBox[{"\"\\"", ",", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"i", " ", "<", "10"}], ",", " ", "\"\<0000\>\"", ",", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"i", " ", "<", " ", "100"}], ",", " ", "\"\<000\>\"", ",", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"i", " ", "<", " ", "1000"}], ",", " ", "\"\<00\>\"", ",", " ", "\"\<0\>\""}], "]"}]}], "]"}]}], "]"}], ",", " ", RowBox[{"ToString", "[", "i", "]"}], ",", " ", "\"\<.png\>\""}], "]"}], ",", " ", RowBox[{"tgraf", "[", RowBox[{"[", "i", "]"}], "]"}], ",", " ", "\"\\""}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], " ", ";"}], "\[IndentingNewLine]", ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"i", ",", " ", RowBox[{"Length", "@", "tgraf"}]}], "}"}]}], "\[IndentingNewLine]", "]"}]], "Input", CellChangeTimes->{{3.857913462523011*^9, 3.857913601079079*^9}},ExpressionUUID->"83146592-6f2d-4f24-98a9-\ a347e6c2a18c"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Dependencia en las condiciones iniciales. Construye una tabla de soluciones \ con diferentes condiciones iniciales. Se omiten los comentarios.\ \>", "Section", CellChangeTimes->{{3.857913910058522*^9, 3.8579139148074055`*^9}, { 3.8580860403958826`*^9, 3.8580860612050543`*^9}},ExpressionUUID->"976a21aa-2f5b-4d36-b242-\ 67ac31b5e568"], Cell[BoxData[{ RowBox[{"Monitor", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tablasol", " ", "=", " ", "\[IndentingNewLine]", RowBox[{"Table", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Lambda]", " ", "=", " ", "1"}], ";", "\[IndentingNewLine]", RowBox[{"\[Mu]", " ", "=", " ", RowBox[{"1", "/", "2"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"\[Phi]10", " ", "=", " ", "2"}], ";", "\[IndentingNewLine]", RowBox[{"\[Phi]20", " ", "=", RowBox[{"3", "+", RowBox[{"n", " ", SuperscriptBox["10", RowBox[{"-", "6"}]]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"\[Phi]1p0", " ", "=", " ", "0"}], ";", "\[IndentingNewLine]", RowBox[{"\[Phi]2p0", " ", "=", " ", "0"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"sol", " ", "=", " ", RowBox[{ RowBox[{"NDSolve", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Lambda]", " ", RowBox[{ RowBox[{"\[Phi]1", "''"}], "[", "t", "]"}]}], " ", "+", " ", RowBox[{"\[Mu]", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Phi]2", "''"}], "[", "t", "]"}], " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{"\[Phi]2", "[", "t", "]"}], " ", "-", " ", RowBox[{"\[Phi]1", "[", "t", "]"}]}], "]"}]}], " ", "-", " ", RowBox[{ SuperscriptBox[ RowBox[{ RowBox[{"\[Phi]2", "'"}], "[", "t", "]"}], "2"], RowBox[{"Sin", "[", RowBox[{ RowBox[{"\[Phi]2", "[", "t", "]"}], " ", "-", " ", RowBox[{"\[Phi]1", "[", "t", "]"}]}], "]"}]}]}], ")"}]}], " ", "+", " ", RowBox[{"Sin", "[", RowBox[{"\[Phi]1", "[", "t", "]"}], "]"}]}], " ", "\[Equal]", " ", "0"}], ",", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Phi]2", "''"}], "[", "t", "]"}], " ", "+", " ", RowBox[{"\[Lambda]", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Phi]1", "''"}], "[", "t", "]"}], " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{"\[Phi]2", "[", "t", "]"}], " ", "-", " ", RowBox[{"\[Phi]1", "[", "t", "]"}]}], "]"}]}], " ", "+", " ", RowBox[{ SuperscriptBox[ RowBox[{ RowBox[{"\[Phi]1", "'"}], "[", "t", "]"}], "2"], RowBox[{"Sin", "[", RowBox[{ RowBox[{"\[Phi]2", "[", "t", "]"}], " ", "-", " ", RowBox[{"\[Phi]1", "[", "t", "]"}]}], "]"}]}]}], ")"}]}], " ", "+", " ", RowBox[{"Sin", "[", RowBox[{"\[Phi]2", "[", "t", "]"}], "]"}]}], " ", "\[Equal]", " ", "0"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Phi]1", "[", "0", "]"}], " ", "\[Equal]", " ", "\[Phi]10"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Phi]2", "[", "0", "]"}], " ", "\[Equal]", " ", "\[Phi]20"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Phi]1", "'"}], "[", "0", "]"}], " ", "\[Equal]", "\[Phi]1p0"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Phi]2", "'"}], "[", "0", "]"}], " ", "\[Equal]", " ", "\[Phi]2p0"}]}], "\[IndentingNewLine]", "}"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"\[Phi]1", ",", " ", "\[Phi]2"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", RowBox[{"tpast", " ", "=", " ", RowBox[{"-", "10"}]}], ",", " ", RowBox[{"tmax", " ", "=", " ", "100"}]}], "}"}], ",", " ", RowBox[{"WorkingPrecision", "\[Rule]", "20"}], ",", " ", RowBox[{"MaxSteps", "\[Rule]", SuperscriptBox["10", "6"]}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"sol", "[", RowBox[{"[", RowBox[{"All", ",", "2"}], "]"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"n", ",", " ", "1", ",", " ", "5"}], "}"}]}], "\[IndentingNewLine]", "]"}]}], ";"}], ",", " ", "\[IndentingNewLine]", "n"}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"r2soltabla", "[", "t_", "]"}], ":=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"\[Lambda]", " ", RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", RowBox[{ RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", "i", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "[", "t", "]"}], "]"}], ",", " ", RowBox[{"-", RowBox[{"Cos", "[", RowBox[{ RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", "i", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "[", "t", "]"}], "]"}]}]}], "}"}]}], " ", "+", " ", RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", RowBox[{ RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", "i", "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "[", "t", "]"}], "]"}], ",", " ", RowBox[{"-", RowBox[{"Cos", "[", RowBox[{ RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", "i", "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "[", "t", "]"}], "]"}]}]}], "}"}]}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", RowBox[{"Length", "@", "tablasol"}]}], "}"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.8579034695664053`*^9, 3.857903478505571*^9}, { 3.8579035185838118`*^9, 3.857903703945401*^9}, {3.857903741382861*^9, 3.8579039176066656`*^9}, {3.8579039707971287`*^9, 3.857904133691585*^9}, 3.857904980231278*^9, {3.8579050644653263`*^9, 3.8579050950901237`*^9}, 3.857905207088758*^9, {3.857905322898846*^9, 3.857905326227235*^9}, 3.857905398685835*^9, {3.85790546690773*^9, 3.8579054912531605`*^9}, 3.8579055438864384`*^9, {3.8579057725756235`*^9, 3.85790577271636*^9}, { 3.8579058511122923`*^9, 3.8579059453120403`*^9}, {3.8579067612682753`*^9, 3.8579067764903345`*^9}, {3.857910764350592*^9, 3.857910856722734*^9}, { 3.857911089022894*^9, 3.857911090928539*^9}, {3.8579111731477613`*^9, 3.857911174975456*^9}, {3.8579114288035126`*^9, 3.857911446877577*^9}, { 3.857911726143736*^9, 3.8579117376116457`*^9}, {3.8579117689702835`*^9, 3.8579117813267546`*^9}, {3.8579130841939487`*^9, 3.8579130852562046`*^9}, 3.8579137684512067`*^9, {3.857917023619361*^9, 3.857917023760105*^9}, { 3.8579174799008636`*^9, 3.857917486131113*^9}, {3.8579175589856043`*^9, 3.85791756048508*^9}, 3.8579176363155613`*^9, {3.858086084630415*^9, 3.8580860908966293`*^9}, {3.858093224672691*^9, 3.8580932420000796`*^9}},ExpressionUUID->"8033d4c8-1e64-412d-8f08-\ 758cf361dd7e"], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"etabla", "[", "t_", "]"}], ":=", " ", RowBox[{"Table", "[", " ", RowBox[{ RowBox[{ RowBox[{ FractionBox["1", "2"], SuperscriptBox["\[Lambda]", "2"], SuperscriptBox[ RowBox[{ RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "'"}], "[", "t", "]"}], "2"]}], " ", "+", " ", RowBox[{ FractionBox["1", "2"], "\[Mu]", RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{ RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "'"}], "[", "t", "]"}], "2"], "+", " ", RowBox[{"2", " ", "\[Lambda]", " ", RowBox[{ RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "'"}], "[", "t", "]"}], " ", RowBox[{ RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "'"}], "[", "t", "]"}], " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "[", "t", "]"}], " ", "-", " ", RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "[", "t", "]"}]}], "]"}]}]}], ")"}]}], " ", "-", " ", RowBox[{"\[Lambda]", " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "[", "t", "]"}], "]"}]}], " ", "-", " ", RowBox[{"\[Mu]", " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "[", "t", "]"}], "]"}]}]}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", RowBox[{"Length", "@", "tablasol"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"e0tabla", " ", "=", " ", RowBox[{"etabla", "[", "0", "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.8579068830607147`*^9, 3.8579068874529467`*^9}, { 3.8579069244620056`*^9, 3.8579069895893354`*^9}, {3.8579070985834217`*^9, 3.8579071283939643`*^9}},ExpressionUUID->"d4b157ec-bd79-40cb-a7cb-\ b1e5afeb106a"], Cell[BoxData[ RowBox[{"Table", "[", RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{ FractionBox[ RowBox[{ RowBox[{"etabla", "[", "t", "]"}], "[", RowBox[{"[", "ipr", "]"}], "]"}], RowBox[{"e0tabla", "[", RowBox[{"[", "ipr", "]"}], "]"}]], "-", "1"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "0", ",", " ", "tmax"}], "}"}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"ipr", ",", " ", "1", ",", " ", RowBox[{"Length", "@", "tablasol"}]}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.8579069962619014`*^9, 3.857907043497657*^9}, { 3.857907137142177*^9, 3.857907154548627*^9}, {3.857910555566225*^9, 3.8579105774421377`*^9}, {3.857910612294837*^9, 3.857910625854333*^9}, { 3.8579137788704753`*^9, 3.8579137790735188`*^9}, {3.8579170108212175`*^9, 3.8579170163067045`*^9}},ExpressionUUID->"f001b523-e069-46e8-81fe-\ 4dd62354d2eb"], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"len", " ", "=", " ", RowBox[{"Length", "@", "tablasol"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{"DynamicModule", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{ "\[Phi]1tabla", ",", " ", "\[Phi]2tabla", ",", " ", "r1tabla", ",", " ", "r2tabla", ",", " ", "trazatabla", ",", " ", "tpr"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"Manipulate", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Phi]1tabla", " ", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "[", "tpr", "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", "len"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"\[Phi]2tabla", " ", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"tablasol", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "[", "tpr", "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", "len"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"r0", " ", "=", " ", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"r1tabla", " ", "=", " ", RowBox[{"\[Lambda]", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", RowBox[{"\[Phi]1tabla", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", RowBox[{"-", RowBox[{"Cos", "[", RowBox[{"\[Phi]1tabla", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", "len"}], "}"}]}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"r2tabla", " ", "=", " ", RowBox[{"r1tabla", " ", "+", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", RowBox[{"\[Phi]2tabla", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", RowBox[{"-", RowBox[{"Cos", "[", RowBox[{"\[Phi]2tabla", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", "len"}], "}"}]}], "]"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"trazatabla", " ", "=", " ", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"r2soltabla", "[", "tt", "]"}], ",", " ", RowBox[{"{", RowBox[{"tt", ",", " ", RowBox[{"tpr", "-", "4"}], ",", " ", "tpr"}], "}"}], ",", " ", RowBox[{"ColorFunction", "\[Rule]", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"x", ",", " ", "y", ",", " ", "u"}], "}"}], ",", RowBox[{"Directive", "[", RowBox[{"Opacity", "[", SqrtBox["u"], "]"}], "]"}]}], "]"}]}]}], "]"}]}], ";"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"gtemp", " ", "=", " ", RowBox[{"Show", "[", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"trazatabla", ","}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Graphics", "[", "\[IndentingNewLine]", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"AbsoluteThickness", "[", "3", "]"}], ",", "\[IndentingNewLine]", "Orange", ",", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"Line", "[", RowBox[{"{", RowBox[{"r0", ",", " ", RowBox[{"r1tabla", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", "len"}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", "Green", ",", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"r1tabla", "[", RowBox[{"[", "i", "]"}], "]"}], ",", " ", RowBox[{"r2tabla", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", "len"}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", "Red", ",", " ", RowBox[{"Disk", "[", RowBox[{"r0", ",", " ", "0.03"}], "]"}], ",", "\[IndentingNewLine]", "Blue", ",", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"Disk", "[", RowBox[{ RowBox[{"r2tabla", "[", RowBox[{"[", "i", "]"}], "]"}], ",", " ", RowBox[{"0.05", " ", FractionBox["\[Mu]", RowBox[{"1", " ", "-", " ", "\[Mu]"}]]}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", "len"}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", "White", ",", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"Disk", "[", RowBox[{ RowBox[{"r1tabla", "[", RowBox[{"[", "i", "]"}], "]"}], ",", " ", "0.05"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", "len"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", "}"}], "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", RowBox[{"PlotRange", "\[Rule]", " ", RowBox[{"1.2", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "1"}], "-", "\[Lambda]"}], ",", " ", RowBox[{"1", "+", "\[Lambda]"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "1"}], "-", "\[Lambda]"}], ",", " ", RowBox[{"1", "+", "\[Lambda]"}]}], "}"}]}], "}"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"Axes", "\[Rule]", " ", "False"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Background", "\[Rule]", " ", "Black"}], ",", " ", "\[IndentingNewLine]", RowBox[{"ImageSize", "\[Rule]", " ", "500"}]}], "\[IndentingNewLine]", "]"}]}]}], ",", " ", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ "tpr", ",", " ", "0", ",", " ", "tmax", ",", " ", "0.01", ",", " ", RowBox[{"Appearance", "\[Rule]", " ", "\"\\""}], ",", " ", RowBox[{"AnimationRate", "\[Rule]", " ", "2"}]}], "}"}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}]}], "Input", CellChangeTimes->{{3.857904270521346*^9, 3.8579042773010054`*^9}, { 3.857904311293007*^9, 3.857904497131488*^9}, {3.857904569366128*^9, 3.8579046586677847`*^9}, {3.8579046985975647`*^9, 3.85790477977835*^9}, { 3.8579051859843473`*^9, 3.857905222682843*^9}, {3.8579059913825183`*^9, 3.8579059967427635`*^9}, {3.857906030104517*^9, 3.857906351481996*^9}, { 3.8579064482281675`*^9, 3.8579065099950333`*^9}, {3.857906655750167*^9, 3.857906689283537*^9}, {3.8579106552536716`*^9, 3.8579107419946404`*^9}, { 3.8579108967289467`*^9, 3.857910922401399*^9}, {3.857910985101269*^9, 3.857911021499003*^9}, {3.8579115633313*^9, 3.857911572488371*^9}, { 3.857911680834355*^9, 3.85791168220903*^9}, {3.857911851701544*^9, 3.857911852092079*^9}, {3.857913096909708*^9, 3.8579130973627234`*^9}, { 3.8579131382749767`*^9, 3.8579131610819263`*^9}, {3.857917415840505*^9, 3.857917439891062*^9}, {3.857917533735338*^9, 3.8579175663752975`*^9}, { 3.857917614395178*^9, 3.8579176173050256`*^9}, {3.857936786891897*^9, 3.857936788409273*^9}, 3.858086075342352*^9, {3.858092764236987*^9, 3.8580927683310204`*^9}},ExpressionUUID->"2f871cd0-fc64-4406-8790-\ ab623e1be351"] }, Open ]] }, WindowSize->{1904, 989}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"11.2 for Microsoft Windows (64-bit) (September 10, 2017)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[558, 20, 200, 5, 28, "Input",ExpressionUUID->"9b358bc8-61b5-45f2-98ae-a2f0f3b6e844"], Cell[761, 27, 824, 22, 48, "Input",ExpressionUUID->"ef9c96d3-dc19-445e-a48a-bb9c27e27e0a"], Cell[CellGroupData[{ Cell[1610, 53, 417, 12, 75, "Section",ExpressionUUID->"3e90e1f8-09cb-4dac-a915-db3d64c0c680"], Cell[2030, 67, 11247, 254, 898, "Input",ExpressionUUID->"8c33cc18-b350-403a-9e82-7edf0bd84570"], Cell[13280, 323, 1101, 25, 59, "Input",ExpressionUUID->"20627f30-a074-4ad5-a85b-de7b1268e934"], Cell[14384, 350, 1308, 31, 93, "Input",ExpressionUUID->"a3bba777-776f-444c-837c-8e5c07fa006c"], Cell[15695, 383, 9210, 191, 1097, "Input",ExpressionUUID->"662f790d-04c3-4b02-8091-44675ec4add0"], Cell[CellGroupData[{ Cell[24930, 578, 500, 9, 53, "Subsection",ExpressionUUID->"f6ad20cb-733c-40e5-a05f-214a4b7b8fc4"], Cell[25433, 589, 8113, 165, 889, "Input",ExpressionUUID->"2411b3af-b027-482f-b721-0a5c44b3d0ad"], Cell[33549, 756, 886, 20, 48, "Input",ExpressionUUID->"146fa7ae-c10c-42ff-81ac-fd08b6480725"] }, Open ]], Cell[CellGroupData[{ Cell[34472, 781, 168, 3, 53, "Subsection",ExpressionUUID->"f747a3e6-cedc-4410-af73-e22310b3bec4"], Cell[34643, 786, 913, 21, 86, "Input",ExpressionUUID->"7f711c5b-f110-477d-ace2-ed407a0389b2"], Cell[35559, 809, 1596, 39, 143, "Input",ExpressionUUID->"83146592-6f2d-4f24-98a9-a347e6c2a18c"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[37204, 854, 350, 7, 67, "Section",ExpressionUUID->"976a21aa-2f5b-4d36-b242-67ac31b5e568"], Cell[37557, 863, 8052, 177, 587, "Input",ExpressionUUID->"8033d4c8-1e64-412d-8f08-758cf361dd7e"], Cell[45612, 1042, 2640, 77, 103, "Input",ExpressionUUID->"d4b157ec-bd79-40cb-a7cb-b1e5afeb106a"], Cell[48255, 1121, 938, 23, 55, "Input",ExpressionUUID->"f001b523-e069-46e8-81fe-4dd62354d2eb"], Cell[49196, 1146, 8568, 193, 793, "Input",ExpressionUUID->"2f871cd0-fc64-4406-8790-ab623e1be351"] }, Open ]] } ] *)