RBC_Mathematica_Imperative.nb

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

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

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     28373,        650]
NotebookOptionsPosition[     26847,        601]
NotebookOutlinePosition[     27235,        618]
CellTagsIndexPosition[     27192,        615]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"Clear", "[", "\"\<Global`*\>\"", "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"Timing", "[", "\n", 
  RowBox[{"(*", " ", 
   RowBox[{"1.", " ", "Calibration"}], "*)"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"\[Alpha]", "=", " ", "0.333333333333"}], ";", 
   "\[IndentingNewLine]", 
   RowBox[{"\[Beta]", "=", " ", "0.95"}], ";", "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"Productivity", " ", "values"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{"vProductivity", " ", "=", 
    RowBox[{"{", 
     RowBox[{
     "0.9792", ",", "0.9896", ",", " ", "1.0000", ",", "1.0106", ",", 
      "1.0212"}], "}"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
   
   RowBox[{"(*", " ", 
    RowBox[{"Transition", " ", "matrix"}], " ", "*)"}], "\[IndentingNewLine]",
    "\[IndentingNewLine]", 
   RowBox[{"mTransition", " ", "=", " ", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "0.9727", ",", "0.0273", ",", "0.0000", ",", "0.0000", ",", "0.0000"}],
        "}"}], ",", "\n", 
      RowBox[{"{", 
       RowBox[{
       "0.0041", ",", "0.9806", ",", "0.0153", ",", "0.0000", ",", "0.0000"}],
        "}"}], ",", "\n", 
      RowBox[{"{", 
       RowBox[{
       "0.0000", ",", "0.0082", ",", "0.9837", ",", "0.0082", ",", "0.0000"}],
        "}"}], ",", "\n", 
      RowBox[{"{", 
       RowBox[{
       "0.0000", ",", "0.0000", ",", "0.0153", ",", "0.9806", ",", "0.0041"}],
        "}"}], ",", "\n", 
      RowBox[{"{", 
       RowBox[{
       "0.0000", ",", "0.0000", ",", "0.0000", ",", "0.0273", ",", "0.9727"}],
        "}"}]}], "}"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"2.", " ", "Steady", " ", "State"}], "*)"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    SubscriptBox["k", "ss"], "=", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"\[Alpha]", "*", "\[Beta]"}], ")"}], 
     FractionBox["1", 
      RowBox[{"1", "-", "\[Alpha]"}]]]}], ";", "\n", "  ", 
   RowBox[{
    SubscriptBox["y", "ss"], "=", 
    SuperscriptBox[
     SubscriptBox["k", "ss"], "\[Alpha]"]}], ";", "\n", "  ", 
   RowBox[{
    SubscriptBox["c", "ss"], "=", 
    RowBox[{
     SubscriptBox["y", "ss"], "-", 
     SubscriptBox["k", "ss"]}]}], ";", "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "We", " ", "generate", " ", "the", " ", "grid", " ", "of", " ", 
     "capital"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{"vGridCapital", "=", 
    RowBox[{"Range", "[", 
     RowBox[{"0.089099143696263", ",", 
      RowBox[{"1.5", "*", 
       SubscriptBox["k", "ss"]}], ",", "0.00001"}], "]"}]}], ";", 
   "\[IndentingNewLine]", 
   RowBox[{"nGridCapital", "=", 
    RowBox[{"Length", "[", "vGridCapital", "]"}]}], ";", "\n", "   ", 
   RowBox[{"nGridProductivity", "=", 
    RowBox[{"Length", "[", "vProductivity", "]"}]}], ";", 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "3.", " ", "Required", " ", "matrices", " ", "and", " ", "vectors"}], 
    "*)"}], "\[IndentingNewLine]", "\n", "  ", 
   RowBox[{"mOutput", "=", 
    RowBox[{"ConstantArray", "[", 
     RowBox[{"0", ",", 
      RowBox[{"{", 
       RowBox[{"nGridCapital", ",", "nGridProductivity"}], "}"}]}], "]"}]}], 
   ";", "\n", "  ", 
   RowBox[{"mValueFunction", "=", 
    RowBox[{"ConstantArray", "[", 
     RowBox[{"0", ",", 
      RowBox[{"{", 
       RowBox[{"nGridCapital", ",", "nGridProductivity"}], "}"}]}], "]"}]}], 
   ";", "\n", "  ", 
   RowBox[{"mValueFunctionNew", "=", 
    RowBox[{"ConstantArray", "[", 
     RowBox[{"0", ",", 
      RowBox[{"{", 
       RowBox[{"nGridCapital", ",", "nGridProductivity"}], "}"}]}], "]"}]}], 
   ";", "\n", "  ", 
   RowBox[{"mPolicyFunction", "=", 
    RowBox[{"ConstantArray", "[", 
     RowBox[{"0", ",", 
      RowBox[{"{", 
       RowBox[{"nGridCapital", ",", "nGridProductivity"}], "}"}]}], "]"}]}], 
   ";", 
   RowBox[{"expectedValueFunction", "=", 
    RowBox[{"ConstantArray", "[", 
     RowBox[{"0", ",", 
      RowBox[{"{", 
       RowBox[{"nGridCapital", ",", "nGridProductivity"}], "}"}]}], "]"}]}], 
   ";", "\n", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"4.", " ", "We", " ", "pre"}], "-", 
     RowBox[{
     "build", " ", "output", " ", "for", " ", "each", " ", "point", " ", "in",
       " ", "the", " ", "grid"}]}], "*)"}], "\[IndentingNewLine]", " ", 
   RowBox[{"mOutput", " ", "=", " ", 
    RowBox[{
     RowBox[{"Transpose", "[", 
      RowBox[{"{", 
       SuperscriptBox["vGridCapital", "\[Alpha]"], "}"}], "]"}], ".", 
     RowBox[{"{", "vProductivity", "}"}]}]}], ";", "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"5.", " ", "Main", " ", "iteration"}], "*)"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"maxDifference", "=", "10.0"}], ";", "\n", "  ", 
   RowBox[{"tolerance", "=", "0.0000001"}], ";", "\n", "  ", 
   RowBox[{"iteration", "=", "0"}], ";", "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"While", " ", "[", 
    RowBox[{
     RowBox[{"maxDifference", ">", "tolerance"}], ",", "\[IndentingNewLine]", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"expectedValueFunction", "=", 
       RowBox[{"mValueFunction", ".", 
        RowBox[{"Transpose", "[", "mTransition", "]"}]}]}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"nProductivity", "=", "1"}], ",", 
        RowBox[{"nProductivity", "<=", "nGridProductivity"}], ",", 
        RowBox[{"nProductivity", "++"}], ",", "\[IndentingNewLine]", 
        "\[IndentingNewLine]", 
        RowBox[{"(*", 
         RowBox[{
         "We", " ", "start", " ", "from", " ", "previous", " ", "choice", " ", 
          RowBox[{"(", 
           RowBox[{
           "monotonicity", " ", "of", " ", "policy", " ", "function"}], 
           ")"}]}], "*)"}], "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"gridCapitalNextPeriod", "=", "1"}], ";", 
         "\[IndentingNewLine]", "\[IndentingNewLine]", 
         RowBox[{"For", " ", "[", 
          RowBox[{
           RowBox[{"nCapital", "=", "1"}], ",", 
           RowBox[{"nCapital", "<=", "nGridCapital"}], ",", 
           RowBox[{"nCapital", "++"}], ",", "\[IndentingNewLine]", 
           "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"valueHighSoFar", "=", 
             RowBox[{"-", "1000.0"}]}], ";", "\[IndentingNewLine]", 
            RowBox[{"capitalChoice", "=", 
             RowBox[{"vGridCapital", "[", 
              RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
            "\[IndentingNewLine]", 
            RowBox[{"For", " ", "[", 
             RowBox[{
              RowBox[{"nCapitalNextPeriod", "=", "gridCapitalNextPeriod"}], 
              ",", " ", 
              RowBox[{"nCapitalNextPeriod", "<=", "nGridCapital"}], ",", 
              RowBox[{"nCapitalNextPeriod", "++"}], ",", 
              "\[IndentingNewLine]", "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{"consumption", "=", " ", 
                RowBox[{
                 RowBox[{"mOutput", "[", 
                  RowBox[{"[", 
                   RowBox[{"nCapital", ",", "nProductivity"}], "]"}], "]"}], 
                 "-", 
                 RowBox[{"vGridCapital", "[", 
                  RowBox[{"[", "nCapitalNextPeriod", "]"}], "]"}]}]}], ";", 
               "\n", "\t       ", 
               RowBox[{"valueProvisional", "=", 
                RowBox[{
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"1", "-", "\[Beta]"}], ")"}], "*", 
                  RowBox[{"Log", "[", "consumption", "]"}]}], "+", 
                 RowBox[{"\[Beta]", "*", 
                  RowBox[{"expectedValueFunction", "[", 
                   RowBox[{"[", 
                    RowBox[{"nCapitalNextPeriod", ",", "nProductivity"}], 
                    "]"}], "]"}]}]}]}], ";", "\n", "\[IndentingNewLine]", 
               RowBox[{"If", "[", 
                RowBox[{
                 RowBox[{"valueProvisional", ">", "valueHighSoFar"}], ",", 
                 "\n", "                   ", 
                 RowBox[{
                  RowBox[{"valueHighSoFar", "=", "valueProvisional"}], ";", 
                  "\n", "                   ", 
                  RowBox[{"capitalChoice", "=", 
                   RowBox[{"vGridCapital", "[", 
                    RowBox[{"[", "nCapitalNextPeriod", "]"}], "]"}]}], ";", 
                  "\n", "                   ", 
                  RowBox[{
                  "gridCapitalNextPeriod", "=", "nCapitalNextPeriod"}], ";"}],
                  ",", "\n", "\t\t ", 
                 RowBox[{"Break", "[", "]"}]}], " ", 
                RowBox[{"(*", 
                 RowBox[{
                 "We", " ", "break", " ", "when", " ", "we", " ", "have", " ",
                   "achieved", " ", "the", " ", "max"}], "*)"}], 
                "\[IndentingNewLine]", "]"}], ";"}]}], "\[IndentingNewLine]", 
             "\n", "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
            RowBox[{
             RowBox[{"mValueFunctionNew", "[", 
              RowBox[{"[", 
               RowBox[{"nCapital", ",", "nProductivity"}], "]"}], "]"}], "=", 
             "valueHighSoFar"}], ";", "\[IndentingNewLine]", 
            RowBox[{
             RowBox[{"mPolicyFunction", "[", 
              RowBox[{"[", 
               RowBox[{"nCapital", ",", "nProductivity"}], "]"}], "]"}], " ", 
             "=", " ", "capitalChoice"}], ";"}]}], "\[IndentingNewLine]", 
          "\[IndentingNewLine]", "]"}], ";"}]}], "\[IndentingNewLine]", 
       "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{"maxDifference", "=", 
       RowBox[{"Max", "[", 
        RowBox[{"Abs", "[", 
         RowBox[{"mValueFunctionNew", "-", "mValueFunction"}], "]"}], "]"}]}],
       ";", "\[IndentingNewLine]", 
      RowBox[{"mValueFunction", "=", "mValueFunctionNew"}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{"iteration", " ", "=", " ", 
       RowBox[{"iteration", "+", "1"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"Mod", "[", 
           RowBox[{"iteration", ",", "10"}], "]"}], "\[Equal]", "0"}], " ", "||",
          "  ", 
         RowBox[{"iteration", " ", "\[Equal]", " ", "1"}]}], ",", 
        RowBox[{"Print", "[", 
         RowBox[{"StringForm", "[", 
          RowBox[{
          "\"\<Iteration = ``, Sup Diff = ``\>\"", ",", " ", "iteration", ",",
            "maxDifference"}], "]"}], "]"}]}], "]"}]}]}], 
    "\[IndentingNewLine]", "\[IndentingNewLine]", "]"}], ";", 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Print", "[", 
    RowBox[{"StringForm", "[", 
     RowBox[{
     "\"\<Iteration = ``, Sup Diff = ``\>\"", ",", " ", "iteration", ",", 
      "maxDifference"}], "]"}], "]"}], ";", "\n", "     ", 
   RowBox[{"Print", "[", 
    RowBox[{"StringForm", "[", 
     RowBox[{"\"\<My check = ``\>\"", ",", 
      RowBox[{"mPolicyFunction", "[", 
       RowBox[{"[", 
        RowBox[{"1000", ",", "3"}], "]"}], "]"}]}], "]"}], "]"}], ";"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]", "]"}]}], "Input",
 CellChangeTimes->{{3.5846306682079163`*^9, 3.584630735618572*^9}, {
   3.584630776252343*^9, 3.584630879969584*^9}, {3.584630923608036*^9, 
   3.584631095703539*^9}, {3.584631266024363*^9, 3.5846313062995777`*^9}, 
   3.584631350801435*^9, {3.5846315403193893`*^9, 3.584631638513904*^9}, {
   3.5846316835027122`*^9, 3.584631696438438*^9}, {3.5846317848009777`*^9, 
   3.584631826256612*^9}, {3.584631890111478*^9, 3.5846319919441977`*^9}, {
   3.584632266434929*^9, 3.584632267489297*^9}, {3.584632830309516*^9, 
   3.584632844450457*^9}, {3.584632886610779*^9, 3.5846328869951878`*^9}, {
   3.584633479970128*^9, 3.584633487384995*^9}, 3.584633632727571*^9, {
   3.584633725008789*^9, 3.584633764102736*^9}, 3.584634072641005*^9, {
   3.584634197282297*^9, 3.584634319219901*^9}, {3.5846343560783463`*^9, 
   3.5846343986028137`*^9}, {3.5846344989571123`*^9, 3.584634503307355*^9}, {
   3.584634567122217*^9, 3.5846345899181137`*^9}, {3.5846346395341463`*^9, 
   3.584634684603146*^9}, {3.584634725724728*^9, 3.584634737018153*^9}, 
   3.584634922641015*^9, {3.5846349793243*^9, 3.584635055686584*^9}, {
   3.584635179164765*^9, 3.584635179485058*^9}, {3.584635209627819*^9, 
   3.584635253564048*^9}, {3.584635300618244*^9, 3.5846353159655952`*^9}, 
   3.5846353648356333`*^9, {3.584635469043132*^9, 3.5846355248219748`*^9}, {
   3.5846355664729433`*^9, 3.5846356157066183`*^9}, {3.58463576651217*^9, 
   3.584635880330806*^9}, {3.584635932822137*^9, 3.5846359941312532`*^9}, {
   3.584636050518354*^9, 3.5846360731446953`*^9}, {3.584636111569872*^9, 
   3.584636153068194*^9}, 3.5846361919155807`*^9, {3.584636228871002*^9, 
   3.584636246953108*^9}, {3.5846363005485067`*^9, 3.584636301530733*^9}, {
   3.584636368732749*^9, 3.584636400257228*^9}, {3.5846364310065947`*^9, 
   3.584636631138157*^9}, {3.584636721127604*^9, 3.584636860417894*^9}, {
   3.5846369047013083`*^9, 3.5846369079224*^9}, {3.584637019219225*^9, 
   3.584637039779928*^9}, {3.584637177221635*^9, 3.5846372155455503`*^9}, {
   3.5846373578185263`*^9, 3.584637360437106*^9}, {3.584637414188767*^9, 
   3.584637507911108*^9}, {3.584637542210458*^9, 3.58463755782082*^9}, {
   3.584637622703624*^9, 3.5846376229493303`*^9}, {3.58463770262567*^9, 
   3.584637709538363*^9}, {3.584637747240013*^9, 3.5846377475438213`*^9}, {
   3.584637806081025*^9, 3.584637809711335*^9}, {3.58463788987212*^9, 
   3.584638009353263*^9}, {3.584638102213913*^9, 3.5846381483551283`*^9}, 
   3.5846383034701643`*^9, {3.584638349388652*^9, 3.584638388403916*^9}, {
   3.584638671081318*^9, 3.584638777613783*^9}, {3.584638934780575*^9, 
   3.5846389391716347`*^9}, {3.584638970633809*^9, 3.584638971460761*^9}, {
   3.584653746235396*^9, 3.584653746977951*^9}, {3.584653903426301*^9, 
   3.58465393237663*^9}, {3.584654108582027*^9, 3.584654216449576*^9}, {
   3.5846543087056427`*^9, 3.584654346428083*^9}, {3.584654435083221*^9, 
   3.584654441722088*^9}, 3.584654534974717*^9, {3.584654567938138*^9, 
   3.584654589452366*^9}, {3.58465466464844*^9, 3.584654689688072*^9}, {
   3.5846548685196543`*^9, 3.584654902357936*^9}, {3.5846549410236893`*^9, 
   3.5846549417554483`*^9}, {3.5846549743315983`*^9, 3.584655050897027*^9}, {
   3.584655170796872*^9, 3.584655206942562*^9}, {3.584655253031302*^9, 
   3.584655270042478*^9}, {3.584655361093749*^9, 3.584655402903014*^9}, {
   3.584655452691639*^9, 3.5846554536941853`*^9}, {3.584655631582487*^9, 
   3.584655632286879*^9}, 3.5846557525665693`*^9, {3.584655904448348*^9, 
   3.58465594506713*^9}, {3.584656069234003*^9, 3.584656077367999*^9}, {
   3.584656109882598*^9, 3.584656165457362*^9}, {3.584656212050213*^9, 
   3.5846562606149178`*^9}, 3.584656458532299*^9, {3.5846565245865307`*^9, 
   3.584656678524419*^9}, {3.5846567097900953`*^9, 3.5846567558360653`*^9}, {
   3.5846569022553377`*^9, 3.584656905517365*^9}, 3.584656943684619*^9, {
   3.584657000686372*^9, 3.584657004094679*^9}, {3.5846570519495363`*^9, 
   3.584657053663232*^9}, {3.584657147987688*^9, 3.584657152052088*^9}, {
   3.584657239090324*^9, 3.584657240824882*^9}, {3.584657359313053*^9, 
   3.584657389462154*^9}, {3.584657454420285*^9, 3.584657466734055*^9}, 
   3.58465776915736*^9, {3.5846579528881197`*^9, 3.584657954553339*^9}, 
   3.584658380612658*^9, {3.584658412891932*^9, 3.584658415562222*^9}, {
   3.584658535968289*^9, 3.584658536964333*^9}, {3.584658734246684*^9, 
   3.5846587474888887`*^9}, {3.584658794955592*^9, 3.584658899563531*^9}, {
   3.584658946479395*^9, 3.5846589468863173`*^9}, {3.5846589896888027`*^9, 
   3.5846590396214457`*^9}, {3.584659077542511*^9, 3.584659081813241*^9}, {
   3.584689657975124*^9, 3.584689674786891*^9}, {3.584689750248437*^9, 
   3.584689751573558*^9}, {3.584689872903369*^9, 3.5846900051924267`*^9}, {
   3.58469007372206*^9, 3.584690080294009*^9}, {3.58469012523381*^9, 
   3.584690194631998*^9}, {3.5846902615698223`*^9, 3.5846902661447153`*^9}, 
   3.584690385174383*^9, {3.584692297709449*^9, 3.584692355559654*^9}, {
   3.584692389682088*^9, 3.5846923954044847`*^9}, {3.58469245957306*^9, 
   3.584692461220242*^9}, 3.584692568293254*^9, {3.584692624018064*^9, 
   3.584692626266182*^9}, {3.584692667634684*^9, 3.5846926712411537`*^9}, {
   3.584692737515375*^9, 3.58469285827071*^9}, {3.584692904360652*^9, 
   3.584692972784519*^9}, 3.584693036193777*^9, {3.5846930746338787`*^9, 
   3.584693099733122*^9}, {3.584693296565711*^9, 3.584693327754486*^9}, {
   3.584693664394355*^9, 3.584693673723754*^9}, {3.584693972230011*^9, 
   3.5846940003639812`*^9}, {3.5848673225051622`*^9, 
   3.5848673760401897`*^9}, {3.584960321475877*^9, 3.584960515046008*^9}, {
   3.58496054638894*^9, 3.584960617204257*^9}, {3.584960675159724*^9, 
   3.58496067634107*^9}, {3.5849613271301727`*^9, 3.58496132903345*^9}, {
   3.584961362717745*^9, 3.584961396946703*^9}, {3.5849614360117493`*^9, 
   3.584961438170683*^9}, {3.584961568622175*^9, 3.584961583567399*^9}, 
   3.584964818314152*^9, {3.585524509126898*^9, 3.585524510150824*^9}, 
   3.5856063311571693`*^9, {3.5856068155354757`*^9, 3.5856068176147623`*^9}, {
   3.585607842673071*^9, 3.585607844813765*^9}, {3.585607913584051*^9, 
   3.585607948202743*^9}, {3.585608002619741*^9, 3.5856080076892633`*^9}, {
   3.5856081122347813`*^9, 3.585608116860621*^9}, {3.585608278630849*^9, 
   3.5856082788220167`*^9}, {3.585608317609104*^9, 3.58560831798881*^9}, {
   3.585608416491969*^9, 3.585608441041428*^9}, {3.585609866009879*^9, 
   3.585609886741611*^9}, {3.585609962103615*^9, 3.585609969700995*^9}, {
   3.5856106230263*^9, 3.585610635796145*^9}, {3.585610874159697*^9, 
   3.585610875236269*^9}, {3.5856110902940693`*^9, 3.5856110927223*^9}, {
   3.585752963106114*^9, 3.585752969528968*^9}, {3.585753074600917*^9, 
   3.585753081652508*^9}, {3.585753164861755*^9, 3.585753181029497*^9}, {
   3.5857532114509153`*^9, 3.585753212402307*^9}, {3.58575338118172*^9, 
   3.585753411710149*^9}, {3.585753444877486*^9, 3.58575346059914*^9}, {
   3.585753492066987*^9, 3.585753500636058*^9}, {3.58575357923068*^9, 
   3.58575358968198*^9}, {3.5857573250649357`*^9, 3.58575732527171*^9}, {
   3.5857644958828783`*^9, 3.585764499712311*^9}, 3.585767566337879*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(1\\), Sup Diff = \
\\!\\(0.05274159340728605`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 1, 0.05274159340728605],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828194092165*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(10\\), Sup Diff = \
\\!\\(0.03134694926582504`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 10, 0.03134694926582504],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828211560184*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(20\\), Sup Diff = \
\\!\\(0.01870345989334099`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 20, 0.01870345989334099],
  Editable->False]], "Print",
 CellChangeTimes->{3.6088282315360622`*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(30\\), Sup Diff = \
\\!\\(0.0111655120339611`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 30, 0.0111655120339611],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828251269455*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(40\\), Sup Diff = \
\\!\\(0.006668541708126696`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 40, 0.006668541708126696],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828271054023*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(50\\), Sup Diff = \
\\!\\(0.0039842927487135915`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 50, 0.0039842927487135915`],
  Editable->False]], "Print",
 CellChangeTimes->{3.6088282907674*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(60\\), Sup Diff = \
\\!\\(0.0023813118039306413`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 60, 0.0023813118039306413`],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828310569254*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(70\\), Sup Diff = \
\\!\\(0.0014236586450968591`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 70, 0.0014236586450968591`],
  Editable->False]], "Print",
 CellChangeTimes->{3.6088283303096247`*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(80\\), Sup Diff = \
\\!\\(0.0008513397747197393`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 80, 0.0008513397747197393],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828350146817*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(90\\), Sup Diff = \
\\!\\(0.0005092051752285665`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 90, 0.0005092051752285665],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828369860079*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(100\\), Sup Diff = \
\\!\\(0.00030462324421454134`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 100, 0.00030462324421454134`],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828389616458*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(110\\), Sup Diff = \
\\!\\(0.000182264856322778`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 110, 0.000182264856322778],
  Editable->False]], "Print",
 CellChangeTimes->{3.6088284094282084`*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(120\\), Sup Diff = \
\\!\\(0.00010906950872624499`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 120, 0.00010906950872624499`],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828429057581*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(130\\), Sup Diff = \
\\!\\(0.00006527643736298216`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 130, 0.00006527643736298216],
  Editable->False]], "Print",
 CellChangeTimes->{3.6088284486767197`*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(140\\), Sup Diff = \
\\!\\(0.000039071082119868095`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 140, 0.000039071082119868095`],
  Editable->False]], "Print",
 CellChangeTimes->{3.6088284684654837`*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(150\\), Sup Diff = \
\\!\\(0.00002338807712010116`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 150, 0.00002338807712010116],
  Editable->False]], "Print",
 CellChangeTimes->{3.6088284881837482`*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(160\\), Sup Diff = \
\\!\\(0.000014008644637297785`\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 160, 0.000014008644637297785`],
  Editable->False]], "Print",
 CellChangeTimes->{3.6088285079257193`*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(170\\), Sup Diff = \
\\!\\(8.391317202982584`*^-6\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 170, 8.391317202982584*^-6],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828527624156*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(180\\), Sup Diff = \
\\!\\(5.026474538150083`*^-6\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 180, 5.026474538150083*^-6],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828547262946*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(190\\), Sup Diff = \
\\!\\(3.010899863653549`*^-6\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 190, 3.010899863653549*^-6],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828567001528*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(200\\), Sup Diff = \
\\!\\(1.8035522481030242`*^-6\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 200, 1.8035522481030242`*^-6],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828586839779*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(210\\), Sup Diff = \
\\!\\(1.0803409160597965`*^-6\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 210, 1.0803409160597965`*^-6],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828606608573*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(220\\), Sup Diff = \
\\!\\(6.471316944534067`*^-7\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 220, 6.471316944534067*^-7],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828626274186*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(230\\), Sup Diff = \
\\!\\(3.876361938104367`*^-7\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 230, 3.876361938104367*^-7],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828645998354*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(240\\), Sup Diff = \
\\!\\(2.3219657929729465`*^-7\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 240, 2.3219657929729465`*^-7],
  Editable->False]], "Print",
 CellChangeTimes->{3.608828665589617*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(250\\), Sup Diff = \
\\!\\(1.3908720952748865`*^-7\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 250, 1.3908720952748865`*^-7],
  Editable->False]], "Print",
 CellChangeTimes->{3.6088286852239428`*^9}],

Cell[BoxData[
 InterpretationBox["\<\"Iteration = \\!\\(257\\), Sup Diff = \
\\!\\(9.716035653806188`*^-8\\)\"\>",
  StringForm["Iteration = ``, Sup Diff = ``", 257, 9.716035653806188*^-8],
  Editable->False]], "Print",
 CellChangeTimes->{3.6088286989960537`*^9}],

Cell[BoxData[
 InterpretationBox["\<\"My check = \\!\\(0.146549143696263`\\)\"\>",
  StringForm["My check = ``", 0.146549143696263],
  Editable->False]], "Print",
 CellChangeTimes->{3.60882869899737*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "588.56787799999995058897184208035469055176`8.790396469398598", ",", 
   "Null"}], "}"}]], "Output",
 CellChangeTimes->{3.608828698999681*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{3.6088280230153637`*^9}]
},
WindowSize->{1088, 765},
WindowMargins->{{156, Automatic}, {Automatic, 3}},
Magnification:>1.5 Inherited,
FrontEndVersion->"9.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (November 20, \
2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 18617, 368, 2516, "Input"],
Cell[CellGroupData[{
Cell[19221, 394, 253, 5, 30, "Print"],
Cell[19477, 401, 255, 5, 30, "Print"],
Cell[19735, 408, 257, 5, 30, "Print"],
Cell[19995, 415, 253, 5, 30, "Print"],
Cell[20251, 422, 257, 5, 30, "Print"],
Cell[20511, 429, 258, 5, 30, "Print"],
Cell[20772, 436, 260, 5, 30, "Print"],
Cell[21035, 443, 262, 5, 30, "Print"],
Cell[21300, 450, 259, 5, 30, "Print"],
Cell[21562, 457, 259, 5, 30, "Print"],
Cell[21824, 464, 264, 5, 30, "Print"],
Cell[22091, 471, 261, 5, 30, "Print"],
Cell[22355, 478, 264, 5, 30, "Print"],
Cell[22622, 485, 265, 5, 30, "Print"],
Cell[22890, 492, 268, 5, 30, "Print"],
Cell[23161, 499, 265, 5, 30, "Print"],
Cell[23429, 506, 268, 5, 30, "Print"],
Cell[23700, 513, 261, 5, 30, "Print"],
Cell[23964, 520, 261, 5, 30, "Print"],
Cell[24228, 527, 261, 5, 30, "Print"],
Cell[24492, 534, 264, 5, 30, "Print"],
Cell[24759, 541, 264, 5, 30, "Print"],
Cell[25026, 548, 261, 5, 30, "Print"],
Cell[25290, 555, 261, 5, 30, "Print"],
Cell[25554, 562, 264, 5, 30, "Print"],
Cell[25821, 569, 266, 5, 30, "Print"],
Cell[26090, 576, 263, 5, 30, "Print"],
Cell[26356, 583, 203, 4, 30, "Print"]
}, Open  ]],
Cell[26574, 590, 184, 5, 77, "Output"]
}, Open  ]],
Cell[26773, 598, 70, 1, 43, "Input"]
}
]
*)

(* End of internal cache information *)