(* Content-type: application/mathematica *)

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

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     32001,        909]
NotebookOptionsPosition[     30121,        841]
NotebookOutlinePosition[     30464,        856]
CellTagsIndexPosition[     30421,        853]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[TextData[{
 StyleBox["Some Examples of Rossler MPSM\n",
  FontWeight->"Bold"],
 "Providing and good initial guess will ensure a very fast convergence to \
true set of cycle points on the Poincar\[EAcute] section. Note that the \
guesses I provide always start on the Poincar\[EAcute] section; my choice was \
x=0.  Take one of the cycle points on the last iteration of the MPSM routine \
and try it as an initial condition for your rossler flow.  You can use the \
cycle visualizer for this."
}], "Text",
 CellChangeTimes->{{3.4323054295360003`*^9, 3.432305483946*^9}, {
  3.4343818084556*^9, 3.4343818437335997`*^9}, {3.4343819769316*^9, 
  3.4343820061085997`*^9}, {3.4343821187046003`*^9, 3.4343822417736*^9}, {
  3.4343841860466003`*^9, 3.4343842874695997`*^9}, {3.4343845198386*^9, 
  3.4343845594336*^9}, {3.4343848273396*^9, 3.4343848310006*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", 
    RowBox[{"2", " ", "cycle"}], "*)"}],
   FontSize->14], "\n", 
  RowBox[{
   RowBox[{
    RowBox[{"initguess", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "3.4", ",", "3.5"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "6.8", ",", "0.02"}], "}"}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"MPSA", "[", 
     RowBox[{"initguess", ",", "2", ",", "3", ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", "0.2", ",", "0.2", ",", 
      "5.7"}], "]"}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.432305719748*^9, 3.432305763299*^9}, 
   3.4323061064519997`*^9, 3.43230627234*^9, 3.432307065898*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "6.822592174552639`*^-17", ",", "3.886329380096718`", ",", 
     "3.8662822525702643`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "6.927627990584362`", ",", "0.06751059961069122`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{{3.4323058447980003`*^9, 3.4323058853570004`*^9}, 
   3.4343837145116*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "6.446700320166513`*^-17", ",", "3.915731495212483`", ",", 
     "3.6910347089178486`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "6.914988037132254`", ",", "0.07592062950036563`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{{3.4323058447980003`*^9, 3.4323058853570004`*^9}, 
   3.4343837145826*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "6.441216742595173`*^-17", ",", "3.9158040540049206`", ",", 
     "3.692833409247779`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "6.9149828484700695`", ",", "0.07571693959676527`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{{3.4323058447980003`*^9, 3.4323058853570004`*^9}, 
   3.4343837146586*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", 
    RowBox[{"3", " ", "cycle"}], "*)"}],
   FontSize->14], "\n", 
  RowBox[{
   RowBox[{
    RowBox[{"initguess", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "4.5", ",", "0.04"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "7.0", ",", "0.1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "3.4", ",", "7.5"}], "}"}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"MPSA", "[", 
     RowBox[{"initguess", ",", "3", ",", "4", ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", "0.2", ",", "0.2", ",", 
      "5.7"}], "]"}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.4343843246205997`*^9, 3.4343843279416*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "3.18746648180364`*^-17", ",", "4.550926387454075`", ",", 
     "0.016868861927562213`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.9265488869537105`*^-16", ",", "7.578450031676873`", ",", 
     "0.14672007174334348`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.139512244140624`*^-16", ",", "2.2123148413892793`", ",", 
     "8.379250515315302`"}], "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, 3.4343838554606*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "6.465082170554957`*^-17", ",", "4.508009312598946`", ",", 
     "0.05210534185831685`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.5266525899125628`*^-16", ",", "7.548765768036511`", ",", 
     "0.13943466000402008`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.082732441986027`*^-16", ",", "2.2807989570056666`", ",", 
     "7.397690796705879`"}], "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, 3.4343838555636*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "6.505510233003501`*^-17", ",", "4.511764765596628`", ",", 
     "0.043486078077556556`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.5088261618385766`*^-16", ",", "7.549970827460294`", ",", 
     "0.14042764116925616`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.085456894849353`*^-16", ",", "2.2782735475013602`", ",", 
     "7.416495410901792`"}], "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, 3.4343838556745996`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "6.505796637344468`*^-17", ",", "4.511776520802628`", ",", 
     "0.043420717553711026`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.5088131248275317`*^-16", ",", "7.549968117232429`", ",", 
     "0.14043237844402642`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.0854479536165534`*^-16", ",", "2.278280788629349`", ",", 
     "7.416482942544927`"}], "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, 3.4343838557855997`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", 
    RowBox[{"3", " ", "cycle"}], "*)"}],
   FontSize->14], "\n", 
  RowBox[{
   RowBox[{
    RowBox[{"initguess", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "2", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "8", ",", "0.1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "6", ",", "7.5"}], "}"}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"MPSA", "[", 
     RowBox[{"initguess", ",", "3", ",", "6", ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", "0.2", ",", "0.2", ",", 
      "5.7"}], "]"}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.4343843246205997`*^9, 3.4343844102746*^9}, {
  3.4343844421216*^9, 3.4343844763866*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.0123955812917586`*^-15"}], ",", "2.8815099267097377`", 
     ",", "9.623464245545987`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "9.3779530429926`*^-16", ",", "5.491515245300818`", ",", 
     "0.11402958841647215`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "7.607842741193942`", ",", "4.896725442730467`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{{3.4343843393466*^9, 3.4343844109636*^9}, {
  3.4343844429256*^9, 3.4343844774486*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.0482802198538626`*^-15"}], ",", "2.985338446046383`", 
     ",", "5.36957918083095`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"9.491024137429142`*^-16", ",", "5.46543532593493`", ",", 
     RowBox[{"-", "0.014902138439324836`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "7.286904493678966`", ",", "0.4150195516986761`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{{3.4343843393466*^9, 3.4343844109636*^9}, {
  3.4343844429256*^9, 3.4343844775615997`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.0561370620886242`*^-15"}], ",", "2.9336221128014373`", 
     ",", "5.665416912823261`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "9.487989390037388`*^-16", ",", "5.463798160430649`", ",", 
     "0.04225298095718515`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "7.294113633787718`", ",", "0.5555513072299314`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{{3.4343843393466*^9, 3.4343844109636*^9}, {
  3.4343844429256*^9, 3.4343844776756*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.0560587706186812`*^-15"}], ",", "2.9328780744730256`", 
     ",", "5.670803605141563`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "9.487942902713644`*^-16", ",", "5.463634949544131`", ",", 
     "0.048540717646444374`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "7.294429431906955`", ",", "0.5561692608713531`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{{3.4343843393466*^9, 3.4343844109636*^9}, {
  3.4343844429256*^9, 3.4343844777896*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.0560587641768304`*^-15"}], ",", "2.9328775591425873`", 
     ",", "5.670805940754447`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "9.487943483764872`*^-16", ",", "5.463634783185518`", ",", 
     "0.04855142379651915`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "7.294429652091356`", ",", "0.5561696289013798`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{{3.4343843393466*^9, 3.4343844109636*^9}, {
  3.4343844429256*^9, 3.4343844779096003`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.0560587641768338`*^-15"}], ",", "2.9328775591421663`", 
     ",", "5.670805940755713`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "9.487943483766016`*^-16", ",", "5.463634783185219`", ",", 
     "0.04855142380538011`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "7.294429652091529`", ",", "0.5561696289020088`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{{3.4343843393466*^9, 3.4343844109636*^9}, {
  3.4343844429256*^9, 3.4343844780326*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", 
    RowBox[{"4", " ", "cycle"}], "*)"}],
   FontSize->14], "\n", 
  RowBox[{
   RowBox[{
    RowBox[{"initguess", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "6.5", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "5.7", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "7", ",", "0.7"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "3.5", ",", "4.7"}], "}"}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"MPSA", "[", 
     RowBox[{"initguess", ",", "4", ",", "4", ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", "0.2", ",", "0.2", ",", 
      "5.7"}], "]"}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.432305719748*^9, 3.432305772915*^9}, 
   3.4323058158859997`*^9, 3.432307410618*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "1.450457692740836`*^-16", ",", "6.280467163821883`", ",", 
     "0.05938314035370501`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "5.614362465012736`", ",", "1.7590921250294123`"}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "6.052665264884778`*^-17", ",", "7.09011223206952`", ",", 
     "0.6859784274292406`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2.54553005099823`*^-17"}], ",", "3.4681990807017913`", ",",
      "4.559706287906965`"}], "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, {3.432307257024*^9, 
   3.4323072716429996`*^9}, {3.43230733588*^9, 3.432307369115*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "1.569608544350317`*^-16", ",", "6.27655454880066`", ",", 
     "0.05792470452763149`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "5.6096724057033445`", ",", "1.6699080756740476`"}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "5.703880315521535`*^-17", ",", "7.089973524318443`", ",", 
     "0.6750027878564403`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2.446402700996127`*^-17"}], ",", "3.4667567315160177`", 
     ",", "4.506188528767023`"}], "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, {3.432307257024*^9, 
   3.4323072716429996`*^9}, {3.43230733588*^9, 3.432307369252*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "1.569399190908553`*^-16", ",", "6.276567398337326`", ",", 
     "0.057941114889608375`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "6.048754030374849`*^-21", ",", "5.6096518331212755`", ",", 
     "1.6704332105794308`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "5.703859392057977`*^-17", ",", "7.089973935023751`", ",", 
     "0.6751139475970422`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2.4463896486148028`*^-17"}], ",", "3.466758711290957`", 
     ",", "4.506217513142985`"}], "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, {3.432307257024*^9, 
   3.4323072716429996`*^9}, {3.43230733588*^9, 3.432307369394*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "1.5693991834013115`*^-16", ",", "6.276567398702727`", ",", 
     "0.05794111499502585`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "6.048754030374849`*^-21", ",", "5.609651832516891`", ",", 
     "1.6704332224425338`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "5.703859396882852`*^-17", ",", "7.089973934963729`", ",", 
     "0.6751139530552761`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2.446389651589922`*^-17"}], ",", "3.4667587115290392`", 
     ",", "4.506217512846343`"}], "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, {3.432307257024*^9, 
   3.4323072716429996`*^9}, {3.43230733588*^9, 3.4323073695360003`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", 
    RowBox[{"5", " ", "cycle"}], "*)"}],
   FontSize->14], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"initguess", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "2.28", ",", "0.058"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "5.59", ",", "1.68"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "7.11", ",", "0.66"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "3.41", ",", "4.6"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "6.2", ",", "0.05"}], "}"}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"MPSA", "[", 
     RowBox[{"initguess", ",", "5", ",", "4", ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", "0.2", ",", "0.2", ",", 
      "5.7"}], "]"}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.432301623572*^9, 3.432301628086*^9}, {3.432301680934*^9,
    3.432301757935*^9}, {3.4323021853710003`*^9, 3.432302207751*^9}, {
   3.4323024903*^9, 3.432302506224*^9}, {3.432302796047*^9, 
   3.432302802599*^9}, 3.4323028332539997`*^9, 3.432303090416*^9, 
   3.4323031362720003`*^9, 3.43230429535*^9, {3.4323043786359997`*^9, 
   3.432304394729*^9}, {3.4323045624519997`*^9, 3.4323045650889997`*^9}, {
   3.432304597917*^9, 3.4323046243719997`*^9}, 3.432304767656*^9, {
   3.432305264398*^9, 3.43230529416*^9}, {3.432305512564*^9, 
   3.432305516486*^9}, {3.4323063430030003`*^9, 3.432306357553*^9}, {
   3.432306436006*^9, 3.4323064655*^9}, {3.4323064986809998`*^9, 
   3.432306575546*^9}, {3.4323067171610003`*^9, 3.432306753742*^9}, {
   3.432306786245*^9, 3.4323068247060003`*^9}, {3.4323068626029997`*^9, 
   3.432306875035*^9}, {3.432306962549*^9, 3.432307042943*^9}, {
   3.432307250801*^9, 3.4323073684820004`*^9}, {3.4323074532*^9, 
   3.4323074942390003`*^9}, {3.432307551959*^9, 3.432307609487*^9}, 
   3.432307657631*^9, {3.4323077244230003`*^9, 3.432307725046*^9}, 
   3.4343845076356*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.1830088243398063`*^-17"}], ",", "3.1453885427497648`", 
     ",", "0.5823971461543026`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.0331403392826883`*^-16", ",", "5.7529392664639625`", ",", 
     "0.2541573146688598`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.404462966510609`*^-16"}], ",", "6.8517560446330705`", 
     ",", "0.8129663265918152`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "7.6905640802231`*^-18", ",", "4.094489432488347`", ",", 
     "3.2424309896777865`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "7.223322832176509`", ",", "0.07586934901006212`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, {3.432307257024*^9, 
   3.4323072716429996`*^9}, {3.43230733588*^9, 3.432307369115*^9}, 
   3.432307495066*^9, {3.4323075525959997`*^9, 3.4323076134779997`*^9}, 
   3.432307658462*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "7.026647445814769`*^-17", ",", "3.179730073188907`", ",", 
     "5.317636706772244`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.0331403392826883`*^-16", ",", "5.743122262935466`", ",", 
     "0.3083222975679813`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.4606617496721124`*^-16"}], ",", "6.836600135851605`", 
     ",", "0.8227468057758998`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.978151778444336`*^-17", ",", "4.142781226318884`", ",", 
     "3.3361433888080985`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "7.197142896372644`", ",", "0.09159386332143944`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, {3.432307257024*^9, 
   3.4323072716429996`*^9}, {3.43230733588*^9, 3.432307369115*^9}, 
   3.432307495066*^9, {3.4323075525959997`*^9, 3.4323076134779997`*^9}, 
   3.432307658635*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "7.019095515527195`*^-17", ",", "3.1164211934102184`", ",", 
     "5.246994615703803`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.0331403392826883`*^-16", ",", "5.744767555603709`", ",", 
     "0.05251120879054627`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.470300356501473`*^-16"}], ",", "6.829139587156327`", ",",
      "0.8269323648868594`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "2.081895265768906`*^-17", ",", "4.1627571014903175`", ",", 
     "3.303823848709906`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "7.218956254749241`", ",", "0.0938313958301701`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, {3.432307257024*^9, 
   3.4323072716429996`*^9}, {3.43230733588*^9, 3.432307369115*^9}, 
   3.432307495066*^9, {3.4323075525959997`*^9, 3.4323076134779997`*^9}, 
   3.432307658809*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "7.039583345750858`*^-17", ",", "3.1167353549075307`", ",", 
     "5.244258048286008`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "1.0331403392826883`*^-16", ",", "5.7446363677291625`", ",", 
     "0.05097250272476415`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.470377606706866`*^-16"}], ",", "6.829125370390308`", ",",
      "0.8269469399406458`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "2.082096328978922`*^-17", ",", "4.162798758664551`", ",", 
     "3.303903340124721`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "7.218822289289799`", ",", "0.09386689928436276`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, {3.432307257024*^9, 
   3.4323072716429996`*^9}, {3.43230733588*^9, 3.432307369115*^9}, 
   3.432307495066*^9, {3.4323075525959997`*^9, 3.4323076134779997`*^9}, 
   3.432307658981*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", 
    RowBox[{"5", " ", "cycle"}], "*)"}],
   FontSize->14], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"initguess", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "6.28", ",", "0.058"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "5.59", ",", "1.68"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "7.11", ",", "0.66"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "3.41", ",", "4.6"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "6.2", ",", "0.05"}], "}"}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"MPSA", "[", 
     RowBox[{"initguess", ",", "5", ",", "4", ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", "0.2", ",", "0.2", ",", 
      "5.7"}], "]"}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.432301623572*^9, 3.432301628086*^9}, {3.432301680934*^9,
    3.432301757935*^9}, {3.4323021853710003`*^9, 3.432302207751*^9}, {
   3.4323024903*^9, 3.432302506224*^9}, {3.432302796047*^9, 
   3.432302802599*^9}, 3.4323028332539997`*^9, 3.432303090416*^9, 
   3.4323031362720003`*^9, 3.43230429535*^9, {3.4323043786359997`*^9, 
   3.432304394729*^9}, {3.4323045624519997`*^9, 3.4323045650889997`*^9}, {
   3.432304597917*^9, 3.4323046243719997`*^9}, 3.432304767656*^9, {
   3.432305264398*^9, 3.43230529416*^9}, {3.432305512564*^9, 
   3.432305516486*^9}, {3.4323063430030003`*^9, 3.432306357553*^9}, {
   3.432306436006*^9, 3.4323064655*^9}, {3.4323064986809998`*^9, 
   3.432306575546*^9}, {3.4323067171610003`*^9, 3.432306753742*^9}, {
   3.432306786245*^9, 3.4323068247060003`*^9}, {3.4323068626029997`*^9, 
   3.432306875035*^9}, {3.432306962549*^9, 3.432307042943*^9}, {
   3.432307250801*^9, 3.4323073684820004`*^9}, {3.4323074532*^9, 
   3.4323074942390003`*^9}, {3.432307551959*^9, 3.432307609487*^9}, {
   3.4343845110746*^9, 3.4343845119716*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "1.8820083974819038`*^-17", ",", "6.308923119211599`", ",", 
     "1.7071412438013303`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "6.872365566252963`*^-18", ",", "5.567639817822242`", ",", 
     "1.70673700931661`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "3.9545500975227124`*^-18", ",", "7.1581518161045565`", ",", 
     "0.6355953887744481`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.7510421722377997`*^-17"}], ",", "3.2870374578583497`", 
     ",", "4.866832908406179`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "6.005963458034405`", ",", "0.05302163090042547`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, {3.432307257024*^9, 
   3.4323072716429996`*^9}, {3.43230733588*^9, 3.432307369115*^9}, 
   3.432307495066*^9, {3.4323075525959997`*^9, 3.4323076134779997`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "5.821257251538306`*^-17"}], ",", "6.303022516923222`", ",",
      "1.143622459407824`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "5.640634626589349`*^-18", ",", "5.565198356893565`", ",", 
     "1.7086408532887467`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "5.386927435049947`*^-18", ",", "7.161451914248247`", ",", 
     "0.6337984888149486`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.906646789089413`*^-17"}], ",", "3.278716851365018`", ",",
      "4.890827406543269`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "5.992294440430584`", ",", "0.05371781508323967`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, {3.432307257024*^9, 
   3.4323072716429996`*^9}, {3.43230733588*^9, 3.432307369115*^9}, 
   3.432307495066*^9, {3.4323075525959997`*^9, 3.432307613651*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "5.798438179772571`*^-17"}], ",", "6.302599058000384`", ",",
      "1.1550952496191929`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "5.641799365026893`*^-18", ",", "5.565239534498521`", ",", 
     "1.7078756221721576`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "5.3907981561984074`*^-18", ",", "7.161377354620107`", ",", 
     "0.6338399154036314`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.903902673828299`*^-17"}], ",", "3.27891425966039`", ",", 
     "4.890453106372313`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "5.992594977711655`", ",", "0.05372965088962727`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, {3.432307257024*^9, 
   3.4323072716429996`*^9}, {3.43230733588*^9, 3.432307369115*^9}, 
   3.432307495066*^9, {3.4323075525959997`*^9, 3.432307613828*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "5.798444864026705`*^-17"}], ",", "6.302598854074691`", ",",
      "1.1551009587648406`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "5.6417811046903774`*^-18", ",", "5.565239557712278`", ",", 
     "1.7078761272650915`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "5.3907783855336514`*^-18", ",", "7.161377317529087`", ",", 
     "0.633839936673508`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.9039038453168214`*^-17"}], ",", "3.2789143575926056`", 
     ",", "4.890452911551636`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "5.9925951247234925`", ",", "0.0537296521775598`"}], 
    "}"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.432306466368*^9, {3.432306500807*^9, 3.432306526132*^9}, 
   3.4323065773640003`*^9, 3.4323067552609997`*^9, {3.4323067892460003`*^9, 
   3.4323068250480003`*^9}, {3.4323068645*^9, 3.432306875916*^9}, {
   3.432306964238*^9, 3.432307043328*^9}, {3.432307257024*^9, 
   3.4323072716429996`*^9}, {3.43230733588*^9, 3.432307369115*^9}, 
   3.432307495066*^9, {3.4323075525959997`*^9, 3.432307614002*^9}}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{627, 597},
WindowMargins->{{266, Automatic}, {-21, Automatic}},
FrontEndVersion->"6.0 for Microsoft Windows (32-bit) (June 19, 2007)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[568, 21, 858, 14, 101, "Text"],
Cell[CellGroupData[{
Cell[1451, 39, 767, 23, 72, "Input"],
Cell[CellGroupData[{
Cell[2243, 66, 395, 11, 26, "Print"],
Cell[2641, 79, 395, 11, 26, "Print"],
Cell[3039, 92, 396, 11, 26, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[3484, 109, 789, 24, 72, "Input"],
Cell[CellGroupData[{
Cell[4298, 137, 756, 19, 48, "Print"],
Cell[5057, 158, 756, 19, 48, "Print"],
Cell[5816, 179, 761, 19, 48, "Print"],
Cell[6580, 200, 761, 19, 48, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[7390, 225, 825, 25, 72, "Input"],
Cell[CellGroupData[{
Cell[8240, 254, 559, 15, 48, "Print"],
Cell[8802, 271, 574, 14, 48, "Print"],
Cell[9379, 287, 562, 15, 48, "Print"],
Cell[9944, 304, 563, 15, 48, "Print"],
Cell[10510, 321, 566, 15, 48, "Print"],
Cell[11079, 338, 562, 15, 48, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[11690, 359, 904, 27, 72, "Input"],
Cell[CellGroupData[{
Cell[12619, 390, 947, 23, 48, "Print"],
Cell[13569, 415, 949, 23, 48, "Print"],
Cell[14521, 440, 977, 24, 48, "Print"],
Cell[15501, 466, 981, 24, 48, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[16531, 496, 2035, 44, 92, "Input"],
Cell[CellGroupData[{
Cell[18591, 544, 1201, 29, 70, "Print"],
Cell[19795, 575, 1183, 29, 70, "Print"],
Cell[20981, 606, 1183, 29, 70, "Print"],
Cell[22167, 637, 1184, 29, 70, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[23400, 672, 1989, 43, 92, "Input"],
Cell[CellGroupData[{
Cell[25414, 719, 1161, 28, 70, "Print"],
Cell[26578, 749, 1168, 28, 70, "Print"],
Cell[27749, 779, 1169, 28, 70, "Print"],
Cell[28921, 809, 1172, 28, 70, "Print"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)