(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 10.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     48268,       1335]
NotebookOptionsPosition[     45979,       1261]
NotebookOutlinePosition[     46321,       1276]
CellTagsIndexPosition[     46278,       1273]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"Quit", "[", "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.718117210165*^9, 3.718117211616*^9}}],

Cell[TextData[{
 "Se trata del sistema de 4 masas iguales unidas por resortes en un \
c\[IAcute]rculo. El lagrangiano (adimensionalizado) es:\n\n",
 Cell[BoxData[
  FormBox[
   RowBox[{"L", " ", "=", " ", 
    RowBox[{
     RowBox[{
      FractionBox["1", "2"], "\[CapitalSigma]", " ", 
      SuperscriptBox[
       OverscriptBox[
        SubscriptBox["x", "i"], "."], "2"]}], "-", 
     RowBox[{
      FractionBox["1", "2"], "[", 
      RowBox[{
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{
          SubscriptBox["x", "1"], "-", 
          SubscriptBox["x", "2"]}], ")"}], "2"], "+", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{
          SubscriptBox["x", "2"], "-", 
          SubscriptBox["x", "3"]}], ")"}], "2"], "+", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{
          SubscriptBox["x", "3"], "-", 
          SubscriptBox["x", "4"]}], ")"}], "2"], "+", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{
          SubscriptBox["x", "4"], "-", 
          SubscriptBox["x", "1"]}], ")"}], "2"]}], "]"}]}]}], 
   TraditionalForm]],
  FontSize->18]
}], "Text",
 CellChangeTimes->{{3.718117021442*^9, 3.7181171537790003`*^9}, {
  3.718127911415*^9, 3.7181279591926003`*^9}}],

Cell[CellGroupData[{

Cell["La matriz energ\[IAcute]a cin\[EAcute]tica:", "Subsection",
 CellChangeTimes->{{3.7181121620472*^9, 3.7181121721736*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n", " ", "=", " ", "4"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", " ", "=", " ", 
   RowBox[{"IdentityMatrix", "[", "n", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"T", " ", "//", " ", "MatrixForm"}]}], "Input",
 CellChangeTimes->{{3.7181121763831997`*^9, 3.7181122030394*^9}, {
  3.718112499836*^9, 3.7181125032986*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["La matriz energ\[IAcute]a potencial:", "Subsection",
 CellChangeTimes->{{3.7181122074960003`*^9, 3.7181122189*^9}}],

Cell["Construida a mano:", "Text",
 CellChangeTimes->{{3.7181122327984*^9, 3.7181122395657997`*^9}, {
  3.718116955623*^9, 3.718116955727*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"V", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"2", ",", " ", 
       RowBox[{"-", "1"}], ",", " ", "0", ",", " ", 
       RowBox[{"-", "1"}]}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "1"}], ",", " ", "2", ",", " ", 
       RowBox[{"-", "1"}], ",", " ", "0"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"0", ",", " ", 
       RowBox[{"-", "1"}], ",", " ", "2", ",", " ", 
       RowBox[{"-", "1"}]}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "1"}], ",", " ", "0", ",", " ", 
       RowBox[{"-", "1"}], ",", " ", "2"}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"V", " ", "//", " ", "MatrixForm"}]}], "Input",
 CellChangeTimes->{{3.7181122461028*^9, 3.71811231591*^9}}],

Cell["\<\
Construida usando comandos f\[AAcute]ciles de generalizar a problemas con m\
\[AAcute]s masas:\
\>", "Text",
 CellChangeTimes->{{3.7181123247084*^9, 3.7181123513782*^9}, {
  3.718116958316*^9, 3.718116958527*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n", " ", "=", " ", "4"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"V", " ", "=", " ", 
   RowBox[{
    RowBox[{"2", " ", 
     RowBox[{"IdentityMatrix", "[", "n", "]"}]}], " ", "+", " ", 
    RowBox[{"SparseArray", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{
         RowBox[{"Band", "[", 
          RowBox[{"{", 
           RowBox[{"1", ",", "2"}], "}"}], "]"}], "\[Rule]", 
         RowBox[{"-", "1"}]}], ",", 
        RowBox[{
         RowBox[{"Band", "[", 
          RowBox[{"{", 
           RowBox[{"2", ",", "1"}], "}"}], "]"}], "\[Rule]", 
         RowBox[{"-", "1"}]}], ",", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"1", ",", "n"}], "}"}], "\[Rule]", 
         RowBox[{"-", "1"}]}], ",", " ", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"n", ",", " ", "1"}], "}"}], "\[Rule]", " ", 
         RowBox[{"-", "1"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n", ",", "n"}], "}"}]}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"V", " ", "//", " ", "MatrixForm"}]}], "Input",
 CellChangeTimes->{{3.7181123547790003`*^9, 3.718112540141*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "El sistema de autovalores y autovectores: ",
 Cell[BoxData[
  FormBox[
   RowBox[{"\[Lambda]", " ", "=", " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"\[Omega]", "/", 
       SubscriptBox["\[Omega]", "o"]}], ")"}], "2"]}], TraditionalForm]]]
}], "Subsection",
 CellFrame->{{0, 0}, {0, 2}},
 CellChangeTimes->{{3.718112568249*^9, 3.7181125755628*^9}, {
  3.7181130785144*^9, 3.7181130788584003`*^9}}],

Cell[CellGroupData[{

Cell["\<\
Primer m\[EAcute]todo: buscando las soluciones de Det[V - \[Lambda] T] = 0 y \
resolviendo para cada autovalor el sistema lineal.\
\>", "Subsubsection",
 CellChangeTimes->{{3.718112896715*^9, 3.7181129099432*^9}, {
   3.7181129545917997`*^9, 3.718112967619*^9}, {3.7181137300154*^9, 
   3.7181137388482*^9}, 3.718117252135*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n", " ", "=", " ", "4"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", " ", "=", " ", 
   RowBox[{"IdentityMatrix", "[", "n", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"V", " ", "=", " ", 
    RowBox[{
     RowBox[{"2", " ", 
      RowBox[{"IdentityMatrix", "[", "n", "]"}]}], " ", "+", " ", 
     RowBox[{"SparseArray", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"Band", "[", 
           RowBox[{"{", 
            RowBox[{"1", ",", "2"}], "}"}], "]"}], "\[Rule]", 
          RowBox[{"-", "1"}]}], ",", 
         RowBox[{
          RowBox[{"Band", "[", 
           RowBox[{"{", 
            RowBox[{"2", ",", "1"}], "}"}], "]"}], "\[Rule]", 
          RowBox[{"-", "1"}]}], ",", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"1", ",", "n"}], "}"}], "\[Rule]", 
          RowBox[{"-", "1"}]}], ",", " ", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"n", ",", " ", "1"}], "}"}], "\[Rule]", " ", 
          RowBox[{"-", "1"}]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"n", ",", "n"}], "}"}]}], "]"}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"solpre", " ", "=", " ", 
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"Det", "[", 
      RowBox[{"V", " ", "-", " ", 
       RowBox[{"\[Lambda]", " ", "T"}]}], "]"}], " ", "\[Equal]", " ", "0"}], 
    ",", " ", "\[Lambda]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"sol\[Lambda]", " ", "=", " ", 
  RowBox[{"solpre", "[", 
   RowBox[{"[", 
    RowBox[{"All", ",", " ", "1", ",", "2"}], "]"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.7181125799172*^9, 3.7181126851879997`*^9}, {
  3.7181138773022003`*^9, 3.7181138775362*^9}, {3.7181230136046*^9, 
  3.7181230240254*^9}}],

Cell["\<\
Para los autovalores con multiplicidad 1, fijamos la primera componente igual \
a 1\
\>", "Text",
 CellChangeTimes->{{3.7181134192496*^9, 3.7181134418008003`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"\[Lambda]", " ", "=", " ", "0"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"x", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{"1", ",", " ", "b", ",", " ", "c", ",", " ", "d"}], "}"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{"solx", " ", "=", " ", 
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"V", " ", "-", " ", 
          RowBox[{
           RowBox[{"sol\[Lambda]", "[", 
            RowBox[{"[", "1", "]"}], "]"}], " ", "T"}]}], ")"}], " ", ".", 
        " ", "x"}], " ", "\[Equal]", " ", "0"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"b", ",", " ", "c", ",", " ", "d"}], "}"}]}], "]"}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"x1", " ", "=", " ", 
    RowBox[{"x", " ", "/.", "  ", 
     RowBox[{"solx", "[", 
      RowBox[{"[", "1", "]"}], "]"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.7181132010453997`*^9, 3.7181132020931997`*^9}, {
   3.7181133010494003`*^9, 3.7181134075664*^9}, {3.7181134443702*^9, 
   3.7181135195332003`*^9}, 3.7181136581452*^9, {3.7181138808416*^9, 
   3.7181138810912*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"\[Lambda]", " ", "=", " ", "4"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"x", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{"1", ",", " ", "b", ",", " ", "c", ",", " ", "d"}], "}"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{"solx", " ", "=", " ", 
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"V", " ", "-", " ", 
          RowBox[{
           RowBox[{"sol\[Lambda]", "[", 
            RowBox[{"[", "4", "]"}], "]"}], " ", "T"}]}], ")"}], " ", ".", 
        " ", "x"}], " ", "\[Equal]", " ", "0"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"b", ",", " ", "c", ",", " ", "d"}], "}"}]}], "]"}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"x2", " ", "=", " ", 
    RowBox[{"x", " ", "/.", "  ", 
     RowBox[{"solx", "[", 
      RowBox[{"[", "1", "]"}], "]"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.7181135312018003`*^9, 3.7181135354282*^9}, {
  3.7181136384124002`*^9, 3.7181136558576*^9}, {3.718113882913*^9, 
  3.718113883158*^9}}],

Cell["\<\
Para el autovalor con multiplicidad 2 debemos fijar dos componentes para \
conseguir dos soluciones L.I. 

Lo m\[AAcute]s c\[OAcute]modo es elegir una componente igual a 1 y otra a 0.\
\>", "Text",
 CellChangeTimes->{{3.7181135490378*^9, 3.7181136110852003`*^9}, {
  3.718117279943*^9, 3.71811728084*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"\[Lambda]", " ", "=", " ", 
    RowBox[{"2", " ", "-", " ", 
     RowBox[{"Primer", " ", "autovector"}]}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"x", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{"1", ",", " ", "0", ",", " ", "c", ",", " ", "d"}], "}"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{"solx", " ", "=", " ", 
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"V", " ", "-", " ", 
          RowBox[{
           RowBox[{"sol\[Lambda]", "[", 
            RowBox[{"[", "2", "]"}], "]"}], " ", "T"}]}], ")"}], " ", ".", 
        " ", "x"}], " ", "\[Equal]", " ", "0"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"c", ",", " ", "d"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"x3", " ", "=", " ", 
    RowBox[{"x", " ", "/.", "  ", 
     RowBox[{"solx", "[", 
      RowBox[{"[", "1", "]"}], "]"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.7181136184004*^9, 3.7181136750746*^9}, {
  3.718113885788*^9, 3.7181138859908*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"\[Lambda]", " ", "=", " ", 
    RowBox[{"2", " ", "-", " ", 
     RowBox[{"Segundo", " ", "autovector"}]}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"x", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{"0", ",", " ", "1", ",", " ", "c", ",", " ", "d"}], "}"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{"solx", " ", "=", " ", 
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"V", " ", "-", " ", 
          RowBox[{
           RowBox[{"sol\[Lambda]", "[", 
            RowBox[{"[", "2", "]"}], "]"}], " ", "T"}]}], ")"}], " ", ".", 
        " ", "x"}], " ", "\[Equal]", " ", "0"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"c", ",", " ", "d"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"x4", "=", " ", 
    RowBox[{"x", " ", "/.", "  ", 
     RowBox[{"solx", "[", 
      RowBox[{"[", "1", "]"}], "]"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.7181136778768*^9, 3.7181136962939997`*^9}, {
  3.7181138879656*^9, 3.7181138881848*^9}}],

Cell["Finalmente, estos son los autovectores:", "Text",
 CellChangeTimes->{{3.7181137785466003`*^9, 3.7181137851064*^9}}],

Cell[BoxData[{"x1", "\[IndentingNewLine]", "x2", "\[IndentingNewLine]", "x3", \
"\[IndentingNewLine]", "x4"}], "Input",
 CellChangeTimes->{{3.718113787127*^9, 3.7181138192771997`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "El sistema de autovalores y autovectores, segundo m\[EAcute]todo: \
reescribiendo la ecuaci\[OAcute]n caracter\[IAcute]stica como \n\n",
 Cell[BoxData[
  FormBox[
   RowBox[{
    RowBox[{"V", ".", 
     SuperscriptBox["T", 
      RowBox[{"-", "1"}]]}], "=", " ", 
    RowBox[{"\[Lambda]", " ", "x"}]}], TraditionalForm]],
  FormatType->"TraditionalForm"],
 "\n\ny usando el comando Eigensystem.  Se lee todo directamente y es f\
\[AAcute]cil de generalizar para cualquier n\[UAcute]mero de masas."
}], "Subsubsection",
 CellFrame->{{0, 0}, {0, 0.5}},
 CellChangeTimes->{{3.7181129370104*^9, 3.7181129421376*^9}, {
  3.7181131065786*^9, 3.7181131610444*^9}, {3.7181137154018*^9, 
  3.7181137567014*^9}, {3.7181138282046003`*^9, 3.7181138511314*^9}, {
  3.718117320142*^9, 3.7181173395439997`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"T", " ", "=", " ", 
   RowBox[{"IdentityMatrix", "[", "n", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"V", " ", "=", " ", 
   RowBox[{
    RowBox[{"2", " ", 
     RowBox[{"IdentityMatrix", "[", "n", "]"}]}], " ", "+", " ", 
    RowBox[{"SparseArray", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{
         RowBox[{"Band", "[", 
          RowBox[{"{", 
           RowBox[{"1", ",", "2"}], "}"}], "]"}], "\[Rule]", 
         RowBox[{"-", "1"}]}], ",", 
        RowBox[{
         RowBox[{"Band", "[", 
          RowBox[{"{", 
           RowBox[{"2", ",", "1"}], "}"}], "]"}], "\[Rule]", 
         RowBox[{"-", "1"}]}], ",", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"1", ",", "n"}], "}"}], "\[Rule]", 
         RowBox[{"-", "1"}]}], ",", " ", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"n", ",", " ", "1"}], "}"}], "\[Rule]", " ", 
         RowBox[{"-", "1"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n", ",", "n"}], "}"}]}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"solESpre", " ", "=", " ", 
  RowBox[{"Eigensystem", "[", 
   RowBox[{"V", " ", ".", " ", 
    RowBox[{"Inverse", "[", "T", "]"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.7181127881205997`*^9, 3.7181127989578*^9}, {
  3.7181130362911997`*^9, 3.7181130510442*^9}, {3.7181138592892*^9, 
  3.7181138610601997`*^9}, {3.7181164316018*^9, 3.7181164318378*^9}, {
  3.7181230313146*^9, 3.718123031848*^9}}],

Cell["\<\
Como el orden en que aparecen los autovalores no suele ser siempre el mismo \
(el autovalor 0 no siempre es el \[UAcute]ltimo valor que da Eigensystem), \
aqu\[IAcute] abajo reordeno autovalores y autovectores de modo que el primero \
de todos corresponda siempre al modo de rotaci\[OAcute]n uniforme, es decir, \
al que tiene \[Lambda] = 0,\
\>", "Text",
 CellChangeTimes->{{3.7181165032166*^9, 3.7181165604018*^9}, {
  3.7181173616870003`*^9, 3.71811737334*^9}}],

Cell[BoxData[{
 RowBox[{"solLambda", " ", "=", 
  RowBox[{"Sort", "[", 
   RowBox[{"solESpre", "[", 
    RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"solEquis", " ", "=", " ", 
  RowBox[{"solESpre", "[", 
   RowBox[{"[", 
    RowBox[{"2", ",", 
     RowBox[{"Ordering", "[", 
      RowBox[{"solESpre", "[", 
       RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "]"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.718113928183*^9, 3.7181139802981997`*^9}, {
  3.7181164453808002`*^9, 3.7181164947734003`*^9}, {3.718117384203*^9, 
  3.718117391517*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "La soluci\[OAcute]n general: \n\n",
 StyleBox["Notar que depende de 2n constantes. Eligiendo ci[[ j, . ]] \
\[Proportional] ",
  FontSize->18],
 Cell[BoxData[
  FormBox[
   SubscriptBox["\[Delta]", "jl"], TraditionalForm]],
  FormatType->"TraditionalForm",
  FontSize->18],
 StyleBox[" se aisla el modo l.",
  FontSize->18],
 "\n\n",
 StyleBox["Notar que sumamos  el vector con las posiciones de equilibrio.\n\n\
Notar que la soluci\[OAcute]n del modo con frecuencia cero ",
  FontSize->18],
 StyleBox["no",
  FontSize->18,
  FontWeight->"Bold",
  FontSlant->"Italic"],
 StyleBox[" es una trigonom\[EAcute]trica sino una lineal.",
  FontSize->18]
}], "Subsection",
 CellChangeTimes->{{3.7181139947524*^9, 3.718113998528*^9}, {
   3.718114124784*^9, 3.7181142346142*^9}, {3.7181146685782003`*^9, 
   3.718114684014*^9}, {3.7181156243406*^9, 3.7181156641493998`*^9}, 
   3.718117415259*^9, {3.718117457639*^9, 3.718117500966*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"\[Omega]s", " ", "=", " ", 
    RowBox[{"Sqrt", "[", "solLambda", "]"}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"las", " ", "frecuencias"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Xequilibrio", " ", "=", " ", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"i", "-", "1"}], ")"}], " ", 
       FractionBox[
        RowBox[{"2", "\[Pi]"}], "n"]}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", " ", "n"}], "}"}]}], "]"}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"las", " ", "posiciones", " ", "de", " ", "equilibrio"}], "*)"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"X", "[", 
    RowBox[{"t_", ",", " ", "ci_"}], "]"}], ":=", "\[IndentingNewLine]", " ", 
   RowBox[{"Xequilibrio", " ", "+", " ", "\[IndentingNewLine]", " ", 
    RowBox[{
     RowBox[{"ci", "[", 
      RowBox[{"[", "1", "]"}], "]"}], 
     RowBox[{"solEquis", "[", 
      RowBox[{"[", "1", "]"}], "]"}], " ", "t"}], " ", "+", 
    "\[IndentingNewLine]", "  ", 
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"solEquis", "[", 
        RowBox[{"[", "i", "]"}], "]"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"ci", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], " ", 
          RowBox[{"Cos", "[", 
           RowBox[{
            RowBox[{"\[Omega]s", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "t"}], "]"}]}], "+", " ", 
         RowBox[{
          RowBox[{"ci", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "2"}], "]"}], "]"}], " ", 
          RowBox[{"Sin", "[", 
           RowBox[{
            RowBox[{"\[Omega]s", "[", 
             RowBox[{"[", "i", "]"}], "]"}], " ", "t"}], "]"}]}]}], ")"}]}], 
      ",", " ", 
      RowBox[{"{", 
       RowBox[{"i", ",", " ", "2", ",", " ", "n"}], "}"}]}], "]"}]}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.718114004856*^9, 3.718114121512*^9}, 
   3.7181141516232*^9, {3.7181146954639997`*^9, 3.7181147459104*^9}, {
   3.7181154228004*^9, 3.7181155034892*^9}, {3.7181156508244*^9, 
   3.7181156677664003`*^9}, {3.7181157022246*^9, 3.7181157479684*^9}, {
   3.7181165744516*^9, 3.7181166175067997`*^9}, {3.71811751937*^9, 
   3.7181175475360003`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Gr\[AAcute]ficos:", "Subsection",
 CellChangeTimes->{{3.7181142516140003`*^9, 3.7181142583162003`*^9}}],

Cell["Primero preparamos el escenario:", "Text",
 CellChangeTimes->{{3.7181142660692*^9, 3.718114274104*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Graphics", "[", "\[IndentingNewLine]", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Hue", "[", 
         FractionBox["i", "n"], "]"}], ",", " ", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"Cos", "[", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"i", "-", "1"}], ")"}], 
              FractionBox[
               RowBox[{"2.", "\[Pi]"}], "n"]}], "]"}], ",", " ", 
            RowBox[{"Sin", "[", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"i", "-", "1"}], ")"}], 
              FractionBox[
               RowBox[{"2.", "\[Pi]"}], "n"]}], "]"}]}], "}"}], ",", "0.1"}], 
         "]"}]}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"i", ",", " ", "1", ",", " ", "n"}], "}"}]}], "]"}], "]"}], 
   ",", "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Graphics", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", 
      RowBox[{"Circle", "[", "]"}]}], "}"}], "]"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"Axes", "\[Rule]", " ", "True"}], ",", "\[IndentingNewLine]", 
   RowBox[{"Ticks", "\[Rule]", " ", "False"}], ",", "\[IndentingNewLine]", 
   RowBox[{"PlotRange", "\[Rule]", " ", 
    RowBox[{"1.2", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}]}]}], 
  "\[IndentingNewLine]", "]"}]], "Input",
 CellChangeTimes->{{3.718114275922*^9, 3.718114553262*^9}, 3.7181149370982*^9}
  ],

Cell["\<\
Ahora reemplazamos las posiciones de los discos por las funciones de t y \
hacemos la animaci\[OAcute]n:

Fijo todas las constantes que acompa\[NTilde]an a los Sin[\[Omega] t] iguales \
a cero, de modo que los controles se reducen a 4 constantes, cada una \
corresponde a la amplitud de un modo. 

Hay un control extra, Amp, que multiplica la amplitud de todos los modos por \
el mismo factor. 

La animaci\[OAcute]n empieza con el modo '4' activo y los otros apagados.\
\>", "Text",
 CellChangeTimes->{{3.7181145846938*^9, 3.7181146075118*^9}, {
  3.7181166784714003`*^9, 3.718116776296*^9}, {3.718117585111*^9, 
  3.718117587377*^9}}],

Cell[BoxData[
 RowBox[{"Manipulate", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Show", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Graphics", "[", "\[IndentingNewLine]", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Hue", "[", 
           FractionBox["i", "n"], "]"}], ",", " ", 
          RowBox[{"Disk", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"punto", " ", "=", " ", 
                RowBox[{
                 RowBox[{"X", "[", 
                  RowBox[{"t", ",", " ", 
                   RowBox[{"Amp", 
                    RowBox[{"{", 
                    RowBox[{"c1", ",", " ", 
                    RowBox[{"{", 
                    RowBox[{"c2", ",", "0"}], "}"}], ",", " ", 
                    RowBox[{"{", 
                    RowBox[{"c3", ",", "0"}], "}"}], ",", 
                    RowBox[{"{", 
                    RowBox[{"c4", ",", "0"}], "}"}]}], "}"}]}]}], "]"}], "[", 
                 RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", " ", 
              RowBox[{"Sin", "[", "punto", "]"}]}], "}"}], ",", "0.1"}], 
           "]"}]}], "}"}], ",", " ", 
        RowBox[{"{", 
         RowBox[{"i", ",", " ", "1", ",", " ", "n"}], "}"}]}], "]"}], "]"}], 
     ",", "\[IndentingNewLine]", "\[IndentingNewLine]", 
     RowBox[{"Graphics", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", 
        RowBox[{"Circle", "[", "]"}]}], "}"}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"Axes", "\[Rule]", " ", "True"}], ",", "\[IndentingNewLine]", 
     RowBox[{"Ticks", "\[Rule]", " ", "False"}], ",", "\[IndentingNewLine]", 
     RowBox[{"PlotRange", "\[Rule]", " ", 
      RowBox[{"1.2", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}]}]}], 
    "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"t", ",", " ", "0.", ",", " ", "\[Infinity]", ",", " ", 
     FractionBox["\[Pi]", "200."], ",", " ", 
     RowBox[{"Appearance", "\[Rule]", " ", "\"\<Open\>\""}], ",", " ", 
     RowBox[{"AnimationRate", "\[Rule]", " ", "2"}]}], "}"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"Amp", ",", " ", "0.25"}], "}"}], ",", " ", "0.", ",", " ", "4",
      ",", " ", "0.01", ",", " ", 
     RowBox[{"Appearance", "\[Rule]", " ", "\"\<Open\>\""}]}], "}"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"c1", ",", " ", "0"}], "}"}], ",", " ", "0", ",", " ", "1", ",",
      " ", "0.01", ",", " ", 
     RowBox[{"Appearance", "\[Rule]", " ", "\"\<Open\>\""}]}], "}"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"c2", ",", " ", "0"}], "}"}], ",", " ", "0", ",", " ", "1", ",",
      " ", "0.01", ",", " ", 
     RowBox[{"Appearance", "\[Rule]", " ", "\"\<Open\>\""}]}], "}"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"c3", ",", " ", "0"}], "}"}], ",", " ", "0", ",", " ", "1", ",",
      " ", "0.01", ",", " ", 
     RowBox[{"Appearance", "\[Rule]", " ", "\"\<Open\>\""}]}], "}"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"c4", ",", " ", "1"}], "}"}], ",", " ", "0", ",", " ", "1", ",",
      " ", "0.01", ",", " ", 
     RowBox[{"Appearance", "\[Rule]", " ", "\"\<Open\>\""}]}], "}"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"TrackedSymbols", "\[RuleDelayed]", " ", 
    RowBox[{"{", 
     RowBox[{
     "Amp", ",", " ", "t", ",", " ", "c1", ",", " ", "c2", ",", " ", "c3", 
      ",", " ", "c3"}], "}"}]}], ",", "\[IndentingNewLine]", 
   RowBox[{"ControlPlacement", "\[Rule]", "Left"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7181146128368*^9, 3.7181146144032*^9}, {
   3.7181148384960003`*^9, 3.7181149341210003`*^9}, {3.7181149957005997`*^9, 
   3.71811503438*^9}, {3.7181150687574*^9, 3.7181152329264*^9}, {
   3.7181152956988*^9, 3.718115315277*^9}, 3.718115357767*^9, 
   3.7181155184274*^9, 3.7181155552335997`*^9, {3.7181158096622*^9, 
   3.7181158385408*^9}, {3.7181159011534*^9, 3.7181159237094*^9}, {
   3.7181159588798*^9, 3.7181160361006002`*^9}, {3.7181166317402*^9, 
   3.7181166341982*^9}, {3.7181166676792*^9, 3.7181166725042*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Para terminar, construyo una tabla con valores predeterminados de las \
condiciones iniciales, para hacer una animaci\[OAcute]n, y en todo caso \
exportar el archivo gif.\
\>", "Subsection",
 CellChangeTimes->{{3.7181178744098*^9, 3.7181179122167997`*^9}}],

Cell[BoxData[
 RowBox[{"Monitor", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"tpuntos", " ", "=", " ", 
     RowBox[{"Table", "[", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"temp", " ", "=", " ", 
        RowBox[{"Show", "[", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"Graphics", "[", "\[IndentingNewLine]", 
           RowBox[{"Table", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"Hue", "[", 
                FractionBox["i", "n"], "]"}], ",", " ", 
               RowBox[{"Disk", "[", 
                RowBox[{
                 RowBox[{"{", 
                  RowBox[{
                   RowBox[{"Cos", "[", 
                    RowBox[{"punto", " ", "=", " ", 
                    RowBox[{
                    RowBox[{"X", "[", 
                    RowBox[{"t", ",", " ", 
                    RowBox[{"0.25", 
                    RowBox[{"{", 
                    RowBox[{"1", ",", " ", 
                    RowBox[{"{", 
                    RowBox[{".4", ",", "0"}], "}"}], ",", " ", 
                    RowBox[{"{", 
                    RowBox[{".9", ",", "0"}], "}"}], ",", 
                    RowBox[{"{", 
                    RowBox[{"1", ",", "0"}], "}"}]}], "}"}]}]}], "]"}], "[", 
                    RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", " ", 
                   RowBox[{"Sin", "[", "punto", "]"}]}], "}"}], ",", "0.1"}], 
                "]"}]}], "}"}], ",", " ", 
             RowBox[{"{", 
              RowBox[{"i", ",", " ", "1", ",", " ", "n"}], "}"}]}], "]"}], 
           "]"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", 
          RowBox[{"Graphics", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", 
             RowBox[{"Circle", "[", "]"}]}], "}"}], "]"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Axes", "\[Rule]", " ", "True"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Ticks", "\[Rule]", " ", "False"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"PlotRange", "\[Rule]", " ", 
           RowBox[{"1.2", " ", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", " ", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}]}]}], 
         "\[IndentingNewLine]", "]"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"t", ",", " ", "0.", ",", 
         RowBox[{"20", " ", "\[Pi]"}], ",", " ", 
         FractionBox["\[Pi]", "50."]}], "}"}]}], "]"}]}], ";"}], ",", " ", 
   RowBox[{"{", 
    RowBox[{"t", ",", " ", "temp"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7181176535298*^9, 3.718117756509*^9}, {
  3.7181178205828*^9, 3.7181178227148*^9}}],

Cell[BoxData[
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"tpuntos", "[", 
    RowBox[{"[", "i", "]"}], "]"}], ",", " ", 
   RowBox[{"{", 
    RowBox[{"i", ",", " ", "1", ",", " ", 
     RowBox[{"Length", "[", "tpuntos", "]"}], ",", " ", "1", ",", " ", 
     RowBox[{"Appearance", "\[Rule]", " ", "\"\<Open\>\""}], ",", " ", 
     RowBox[{"AnimationRate", "\[Rule]", " ", "25"}]}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.7181177693710003`*^9, 3.7181178051338*^9}}],

Cell["\<\
Aqu\[IAcute] exporto los frames y hago el gif aparte. Tambi\[EAcute]n se \
puede exportar directamente como archivo gif toda la tabla de \
gr\[AAcute]ficos.\
\>", "Text",
 CellChangeTimes->{{3.718122873719*^9, 3.7181228987826*^9}}],

Cell[BoxData[
 RowBox[{"dir", " ", "=", " ", 
  RowBox[{
  "CreateDirectory", "[", "\"\<..\\\\Desktop\\\\masas\\\\\>\"", 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.7181205002876*^9, 3.718120568009*^9}}],

Cell[BoxData[
 RowBox[{"Monitor", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"i", " ", "=", " ", "1"}], ",", " ", 
      RowBox[{"i", " ", "\[LessEqual]", " ", 
       RowBox[{"Length", "[", "tpuntos", "]"}]}], ",", " ", 
      RowBox[{"++", "i"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"itemp", " ", "=", " ", 
        RowBox[{"tpuntos", "[", 
         RowBox[{"[", "i", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"Export", "[", 
        RowBox[{
         RowBox[{"StringJoin", "[", 
          RowBox[{"dir", ",", " ", "\"\<frame_\>\"", ",", " ", 
           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\>\""}], 
          "]"}], ",", " ", "itemp", ",", " ", "\"\<PNG\>\""}], "]"}], ";"}]}],
      "\[IndentingNewLine]", "]"}], ";"}], ",", "i"}], "]"}]], "Input",
 CellChangeTimes->{{3.7181205827862*^9, 3.71812069869*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["En 3D", "Subsection",
 CellFrame->{{0, 0}, {0, 2}},
 CellChangeTimes->{{3.7181198662757998`*^9, 3.7181198675218*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"antialias", "[", 
    RowBox[{"g_", ",", 
     RowBox[{"n_:", " ", "1"}]}], "]"}], ":=", 
   RowBox[{"ImageResize", "[", 
    RowBox[{
     RowBox[{"Rasterize", "[", 
      RowBox[{"g", ",", "\"\<Image\>\"", ",", 
       RowBox[{"ImageResolution", "\[Rule]", 
        RowBox[{"n", " ", "72"}]}], ",", " ", 
       RowBox[{"Background", "\[Rule]", "None"}]}], "]"}], ",", 
     RowBox[{"Scaled", "[", 
      RowBox[{"1", "/", "n"}], "]"}]}], "]"}]}], ";"}]], "Input"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"aro", " ", "=", " ", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"RGBColor", "[", 
       RowBox[{"0.16", ",", "0.5", ",", "0.83"}], "]"}], ",", 
      RowBox[{"Specularity", "[", 
       RowBox[{"White", ",", "100"}], "]"}], ",", " ", 
      RowBox[{"Tube", "[", 
       RowBox[{
        RowBox[{"Table", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"Cos", "[", "x", "]"}], ",", " ", 
            RowBox[{"Sin", "[", "x", "]"}], ",", " ", "0"}], "}"}], ",", " ", 
          RowBox[{"{", 
           RowBox[{"x", ",", " ", "0", ",", " ", 
            RowBox[{"2.", "\[Pi]"}], ",", " ", 
            RowBox[{"2", 
             FractionBox["\[Pi]", "200."]}]}], "}"}]}], "]"}], ",", " ", 
        "0.02"}], "]"}]}], "}"}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"color", " ", "=", " ", 
    RowBox[{"{", 
     RowBox[{
     "Orange", ",", " ", "Red", ",", " ", "Yellow", ",", " ", "Green"}], 
     "}"}]}], ";"}], "\[IndentingNewLine]"}], "\n", 
 RowBox[{
  RowBox[{"marco", " ", "=", " ", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"RGBColor", "[", 
         RowBox[{"0.16", ",", "0.5", ",", "0.83"}], "]"}], ",", " ", 
        RowBox[{"Line", "[", 
         RowBox[{"0.75", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", " ", 
            RowBox[{"{", 
             RowBox[{"1", ",", "0", ",", "0"}], "}"}]}], "}"}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Line", "[", 
         RowBox[{"0.75", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", 
              RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", " ", 
            RowBox[{"{", 
             RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], "}"}]}], "]"}]}], 
       "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"RGBColor", "[", 
         RowBox[{"0.16", ",", "0.5", ",", "0.83"}], "]"}], ",", " ", 
        RowBox[{"Line", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"-", "1.5"}], ",", "0", ",", "0"}], "}"}], ",", " ", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"-", "1.25"}], ",", "0", ",", "0"}], "}"}]}], "}"}], 
         "]"}], ",", "\[IndentingNewLine]", 
        RowBox[{"Line", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"0", ",", 
             RowBox[{"-", "1.5"}], ",", "0"}], "}"}], ",", " ", 
           RowBox[{"{", 
            RowBox[{"0", ",", 
             RowBox[{"-", "1.25"}], ",", "0"}], "}"}]}], "}"}], "]"}]}], 
       "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"RGBColor", "[", 
         RowBox[{"0.16", ",", "0.5", ",", "0.83"}], "]"}], ",", " ", 
        RowBox[{"Line", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"1.5", ",", "0", ",", "0"}], "}"}], ",", " ", 
           RowBox[{"{", 
            RowBox[{"1.25", ",", "0", ",", "0"}], "}"}]}], "}"}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Line", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"0", ",", "1.5", ",", "0"}], "}"}], ",", " ", 
           RowBox[{"{", 
            RowBox[{"0", ",", "1.25", ",", "0"}], "}"}]}], "}"}], "]"}]}], 
       "}"}]}], "}"}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.7181199894308*^9, 3.7181204485156*^9}, 
   3.7181208674309998`*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"tpr", " ", "=", " ", "0.2"}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"graf", " ", "=", " ", 
  RowBox[{"Show", "[", "\[IndentingNewLine]", 
   RowBox[{
   "marco", ",", "\[IndentingNewLine]", "aro", ",", "\[IndentingNewLine]", 
    RowBox[{"Graphics3D", "[", "\[IndentingNewLine]", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", " ", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"color", "[", 
           RowBox[{"[", "i", "]"}], "]"}], ",", 
          RowBox[{"Specularity", "[", 
           RowBox[{"White", ",", "100"}], "]"}], ",", " ", 
          RowBox[{"Sphere", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"punto", " ", "=", " ", 
                RowBox[{
                 RowBox[{"X", "[", 
                  RowBox[{"tpr", ",", " ", 
                   RowBox[{"0.25", 
                    RowBox[{"{", 
                    RowBox[{"1", ",", " ", 
                    RowBox[{"{", 
                    RowBox[{".4", ",", "0"}], "}"}], ",", " ", 
                    RowBox[{"{", 
                    RowBox[{".9", ",", "0"}], "}"}], ",", 
                    RowBox[{"{", 
                    RowBox[{"1", ",", "0"}], "}"}]}], "}"}]}]}], "]"}], "[", 
                 RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", " ", 
              RowBox[{"Sin", "[", "punto", "]"}], ",", "0"}], "}"}], ",", 
            "0.1"}], "]"}]}], "}"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"i", ",", " ", "1", ",", " ", "n"}], "}"}]}], "]"}], "]"}], 
    ",", "\[IndentingNewLine]", 
    RowBox[{"Axes", "\[Rule]", " ", "False"}], ",", "\[IndentingNewLine]", 
    RowBox[{"Ticks", "\[Rule]", " ", "False"}], ",", "\[IndentingNewLine]", 
    RowBox[{"Boxed", "\[Rule]", " ", "False"}], ",", "\[IndentingNewLine]", 
    RowBox[{"Background", "\[Rule]", " ", "Black"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{"PlotRange", "\[Rule]", " ", 
     RowBox[{"2.5", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", " ", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", " ", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}]}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{"ViewVector", "\[Rule]", " ", 
     RowBox[{"{", 
      RowBox[{"5", ",", "2", ",", " ", "3"}], "}"}]}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{"ViewVertical", "\[Rule]", " ", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], ",", "\[IndentingNewLine]", 
    RowBox[{"ViewAngle", "\[Rule]", 
     RowBox[{"30", "\[Degree]"}]}], ",", "\[IndentingNewLine]", 
    RowBox[{"ImageSize", "\[Rule]", " ", "800"}]}], "\[IndentingNewLine]", 
   "]"}]}]}], "Input",
 CellChangeTimes->{{3.7181199894308*^9, 3.7181204485156*^9}, {
  3.7181208674309998`*^9, 3.7181208734224*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"$HistoryLength", " ", "=", " ", "0"}], ";", " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "para", " ", "que", " ", "no", " ", "guarde", " ", "resultados", " ", "y", 
    " ", "sature", " ", "la", " ", "memoria", " ", "al", " ", "hacer", " ", 
    "los", " ", "gr\[AAcute]ficos", " ", "3", "D"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.7181209862511997`*^9, 3.7181209903424*^9}, {
  3.7181212258906*^9, 3.7181212477846003`*^9}}],

Cell[BoxData[
 RowBox[{"dir", " ", "=", " ", 
  RowBox[{
  "SetDirectory", "[", "\"\<..\\\\Desktop\\\\masas\\\\\>\"", "]"}]}]], "Input",
 CellChangeTimes->{{3.7181219232606*^9, 3.7181219282284*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"nframes", " ", "=", " ", "500"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Monitor", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"i", " ", "=", " ", "80"}], ",", " ", 
      RowBox[{"i", " ", "\[LessEqual]", " ", "nframes"}], ",", " ", 
      RowBox[{"++", "i"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"tpr", " ", "=", " ", 
        RowBox[{
         FractionBox[
          RowBox[{"10.", "\[Pi]"}], "nframes"], "i"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"graf", " ", "=", " ", 
        RowBox[{"Show", "[", "\[IndentingNewLine]", 
         RowBox[{
         "marco", ",", "\[IndentingNewLine]", "aro", ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Graphics3D", "[", "\[IndentingNewLine]", 
           RowBox[{"Table", "[", 
            RowBox[{
             RowBox[{"{", " ", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"color", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], ",", 
                RowBox[{"Specularity", "[", 
                 RowBox[{"White", ",", "100"}], "]"}], ",", " ", 
                RowBox[{"Sphere", "[", 
                 RowBox[{
                  RowBox[{"{", 
                   RowBox[{
                    RowBox[{"Cos", "[", 
                    RowBox[{"punto", " ", "=", " ", 
                    RowBox[{
                    RowBox[{"X", "[", 
                    RowBox[{"tpr", ",", " ", 
                    RowBox[{"0.25", 
                    RowBox[{"{", 
                    RowBox[{"1", ",", " ", 
                    RowBox[{"{", 
                    RowBox[{".4", ",", "0"}], "}"}], ",", " ", 
                    RowBox[{"{", 
                    RowBox[{".9", ",", "0"}], "}"}], ",", 
                    RowBox[{"{", 
                    RowBox[{"1", ",", "0"}], "}"}]}], "}"}]}]}], "]"}], "[", 
                    RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", " ", 
                    RowBox[{"Sin", "[", "punto", "]"}], ",", "0"}], "}"}], 
                  ",", "0.1"}], "]"}]}], "}"}], "}"}], ",", " ", 
             RowBox[{"{", 
              RowBox[{"i", ",", " ", "1", ",", " ", "n"}], "}"}]}], "]"}], 
           "]"}], ",", "\[IndentingNewLine]", 
          RowBox[{"Axes", "\[Rule]", " ", "False"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Ticks", "\[Rule]", " ", "False"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Boxed", "\[Rule]", " ", "False"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Background", "\[Rule]", " ", "Black"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"PlotRange", "\[Rule]", " ", 
           RowBox[{"2.5", " ", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", " ", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", " ", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"ViewVector", "\[Rule]", " ", 
           RowBox[{"{", 
            RowBox[{"5", ",", "2", ",", " ", "3"}], "}"}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"ViewVertical", "\[Rule]", " ", 
           RowBox[{"{", 
            RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"ViewAngle", "\[Rule]", 
           RowBox[{"30", "\[Degree]"}]}], ",", "\[IndentingNewLine]", 
          RowBox[{"ImageSize", "\[Rule]", " ", "800"}]}], 
         "\[IndentingNewLine]", "]"}]}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{"itemp", " ", "=", " ", 
        RowBox[{"antialias", "[", 
         RowBox[{"graf", ",", " ", "2"}], "]"}]}], ";", "\[IndentingNewLine]",
        "\[IndentingNewLine]", 
       RowBox[{"Export", "[", 
        RowBox[{
         RowBox[{"StringJoin", "[", 
          RowBox[{"dir", ",", " ", "\"\<\\\\frame_3D_\>\"", ",", " ", 
           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\>\""}], 
          "]"}], ",", " ", "itemp", ",", " ", "\"\<PNG\>\""}], "]"}], ";"}]}],
      "\[IndentingNewLine]", "]"}], ";"}], ",", 
   RowBox[{"{", 
    RowBox[{"i", ",", " ", "itemp"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.7181210042460003`*^9, 3.7181211385088*^9}, {
  3.7181218916918*^9, 3.7181218917854*^9}, {3.7181219532714*^9, 
  3.718121953443*^9}}]
}, Open  ]]
},
WindowSize->{1440, 830},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
FrontEndVersion->"10.0 for Microsoft Windows (32-bit) (July 1, 2014)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 130, 3, 31, "Input"],
Cell[691, 25, 1249, 39, 82, "Text"],
Cell[CellGroupData[{
Cell[1965, 68, 127, 1, 43, "Subsection"],
Cell[2095, 71, 379, 8, 72, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2511, 84, 121, 1, 43, "Subsection"],
Cell[2635, 87, 143, 2, 30, "Text"],
Cell[2781, 91, 832, 23, 52, "Input"],
Cell[3616, 116, 223, 5, 30, "Text"],
Cell[3842, 123, 1173, 34, 72, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5052, 162, 431, 12, 53, "Subsection"],
Cell[CellGroupData[{
Cell[5508, 178, 337, 6, 34, "Subsubsection"],
Cell[5848, 186, 1839, 51, 132, "Input"],
Cell[7690, 239, 173, 4, 30, "Text"],
Cell[7866, 245, 1185, 32, 92, "Input"],
Cell[9054, 279, 1107, 31, 92, "Input"],
Cell[10164, 312, 315, 7, 68, "Text"],
Cell[10482, 321, 1111, 31, 92, "Input"],
Cell[11596, 354, 1107, 31, 92, "Input"],
Cell[12706, 387, 121, 1, 30, "Text"],
Cell[12830, 390, 184, 2, 92, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13051, 397, 814, 18, 133, "Subsubsection"],
Cell[13868, 417, 1505, 41, 72, "Input"],
Cell[15376, 460, 474, 8, 49, "Text"],
Cell[15853, 470, 573, 14, 52, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[16475, 490, 946, 26, 183, "Subsection"],
Cell[17424, 518, 2379, 65, 168, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19840, 588, 109, 1, 43, "Subsection"],
Cell[19952, 591, 109, 1, 30, "Text"],
Cell[20064, 594, 1796, 50, 231, "Input"],
Cell[21863, 646, 646, 15, 144, "Text"],
Cell[22512, 663, 4600, 109, 440, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27149, 777, 266, 5, 67, "Subsection"],
Cell[27418, 784, 2891, 69, 292, "Input"],
Cell[30312, 855, 483, 11, 31, "Input"],
Cell[30798, 868, 241, 5, 30, "Text"],
Cell[31042, 875, 204, 5, 31, "Input"],
Cell[31249, 882, 1443, 33, 92, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32729, 920, 123, 2, 53, "Subsection"],
Cell[32855, 924, 513, 14, 31, "Input"],
Cell[33371, 940, 3756, 105, 206, "Input"],
Cell[37130, 1047, 3049, 74, 352, "Input"],
Cell[40182, 1123, 462, 9, 31, "Input"],
Cell[40647, 1134, 198, 4, 31, "Input"],
Cell[40848, 1140, 5115, 118, 503, "Input"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)