(************** Content-type: application/mathematica ************** Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 22943, 710]*) (*NotebookOutlinePosition[ 23649, 734]*) (* CellTagsIndexPosition[ 23605, 730]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["Basic RBC", "Title"], Cell["Hansen Model of Indivisible Labour", "Subtitle"], Cell[BoxData[ \(Off[General::spell]; \ Off[General::spell1];\)], "Input"], Cell[CellGroupData[{ Cell["Set Working Directory, Load Package and Parameters", "Section"], Cell[BoxData[ \(\(SetDirectory["\"];\)\)], "Input"], Cell[BoxData[ \(<< SymbOpt`\)], "Input"], Cell[BoxData[ \(\(dim = {2, 1, 1, 1, 3};\)\)], "Input"], Cell[BoxData[ \(\[Beta] = 0.99; \ \[Delta] = 0.025; \[Alpha] = 0.36; \[Rho]1 = 0.9; Amean = Exp[1]; \ h = 0.333;\)], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["Relevant Functions", "Section"], Cell[BoxData[ \(\(u = Log[ct] - bcap*ht;\)\)], "Input"], Cell[BoxData[ \(\(Y = Log[At]*kt\^\[Alpha]*ht\^\(1 - \[Alpha]\);\)\)], "Input"], Cell[BoxData[ \(\(Y1 = Log[At1]*kt1\^\[Alpha]*ht1\^\(1 - \[Alpha]\);\)\)], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["First Order Nonlinear Dynamic System", "Section"], Cell[BoxData[ \(\(f1 = \[PartialD]\_ct u - \[Lambda]t;\)\)], "Input"], Cell[BoxData[ \(\(f2 = \[PartialD]\_ht u + \[Lambda]t*\((1 - \[Alpha])\)* Log[At]*\((kt\/ht)\)\^\[Alpha]\ ;\)\)], "Input"], Cell[BoxData[ \(\(f3 = \[Lambda]t\/\(\[Beta]*\[Lambda]t1\) - 1 + \[Delta] - \[Alpha]* Log[At1]*\((ht1\/kt1)\)\^\(1 - \[Alpha]\);\)\)], "Input"], Cell[BoxData[ \(\(f4 = ct + kt1 - \((1 - \[Delta])\) kt - Log[At]*kt\^\[Alpha]*ht\^\(1 - \[Alpha]\);\)\)], "Input"], Cell[BoxData[ \(\(f5 = yt - Y;\)\)], "Input"], Cell[BoxData[ \(\(f6 = ivt + ct - Y;\)\)], "Input"], Cell[BoxData[ \(\(f7 = LPt - yt\/ht;\)\)], "Input"], Cell[BoxData[ \(\(f8 = Log[At1] - \[Rho]1\ *Log[At] - \((1 - \[Rho]1)\)* Log[Amean];\)\)], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["Defining Functions and Lists", "Section"], Cell[BoxData[ \(\(f = {f1, f2, f3, f4, f5, f6, f7, f8};\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(MatrixForm[f]\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {\(1\/ct - \[Lambda]t\)}, {\(\(-bcap\) + 0.64`\ \((kt\/ht)\)\^0.36`\ \[Lambda]t\ Log[At]\)}, {\(\(-0.975`\) + \(1.0101010101010102`\ \[Lambda]t\)\/\[Lambda]t1 \ - 0.36`\ \((ht1\/kt1)\)\^0.64`\ Log[At1]\)}, {\(ct - 0.975`\ kt + kt1 - ht\^0.64`\ kt\^0.36`\ Log[At]\)}, {\(yt - ht\^0.64`\ kt\^0.36`\ Log[At]\)}, {\(ct + ivt - ht\^0.64`\ kt\^0.36`\ Log[At]\)}, {\(LPt - yt\/ht\)}, {\(\(-0.09999999999999998`\) - 0.9`\ Log[At] + Log[At1]\)} }], "\[NoBreak]", ")"}], MatrixForm[ { Plus[ Power[ ct, -1], Times[ -1, \[Lambda]t]], Plus[ Times[ -1, bcap], Times[ .64000000000000001, Power[ Times[ Power[ ht, -1], kt], .35999999999999999], \[Lambda]t, Log[ At]]], Plus[ -.97499999999999998, Times[ 1.0101010101010102, \[Lambda]t, Power[ \[Lambda]t1, -1]], Times[ -.35999999999999999, Power[ Times[ ht1, Power[ kt1, -1]], .64000000000000001], Log[ At1]]], Plus[ ct, Times[ -.97499999999999998, kt], kt1, Times[ -1, Power[ ht, .64000000000000001], Power[ kt, .35999999999999999], Log[ At]]], Plus[ yt, Times[ -1, Power[ ht, .64000000000000001], Power[ kt, .35999999999999999], Log[ At]]], Plus[ ct, ivt, Times[ -1, Power[ ht, .64000000000000001], Power[ kt, .35999999999999999], Log[ At]]], Plus[ LPt, Times[ -1, Power[ ht, -1], yt]], Plus[ -.099999999999999978, Times[ -.90000000000000002, Log[ At]], Log[ At1]]}]]], "Output"] }, Open ]], Cell[BoxData[ \(\(vt = {ct, ht, kt, \[Lambda]t, At, yt, ivt, LPt};\)\)], "Input"], Cell[BoxData[ \(\(vt1 = {ct1, ht1, kt1, \[Lambda]t1, At1, yt1, ivt1, LPt1};\)\)], "Input"], Cell[BoxData[ \(\(vht = {cht, hht, kht, \[Lambda]ht, Aht, yht, ivht, LPht};\)\)], "Input"], Cell[BoxData[ \(\(vht1 = {cht1, hht1, kht1, \[Lambda]ht1, Aht1, yht1, ivht1, LPht1};\)\)], "Input"], Cell[BoxData[ \(\(v = {c, h, k, \[Lambda], A, y, iv, LP};\)\)], "Input"], Cell[BoxData[ \(\(v\&^ = {c\&^, h\&^, k\&^, \[Lambda]\&^, A\&^, y\&^, iv\&^, LP\&^};\)\)], "Input"], Cell[BoxData[ \(\(ceros = Table[0, {Length[v]}];\)\)], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["Computation of the Steady State", "Section"], Cell[BoxData[ \(\(nvar = Length[v];\)\)], "Input"], Cell[BoxData[ \(\(xss0 = Table[vt[\([i]\)] \[Rule] v[\([i]\)], {i, nvar}];\)\)], "Input"], Cell[BoxData[ \(\(xss1 = Table[vt1[\([i]\)] \[Rule] v[\([i]\)], {i, nvar}];\)\)], "Input"], Cell[BoxData[ \(\(xss\ = \ Flatten[{xss0, xss1}];\)\)], "Input"], Cell[BoxData[ \(\(fss = f /. xss;\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(MatrixForm[fss]\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {\(1\/c - \[Lambda]\)}, {\(\(-bcap\) + 0.9508238129969964`\ k\^0.36`\ \[Lambda]\ Log[A]\)}, {\(\(\(0.03510101010101019`\)\(\[InvisibleSpace]\)\) - 0.17810118547199988`\ \((1\/k)\)\^0.64`\ Log[A]\)}, {\(c + 0.025000000000000022`\ k - 0.4947255151999997`\ k\^0.36`\ Log[A]\)}, {\(y - 0.4947255151999997`\ k\^0.36`\ Log[A]\)}, {\(c + iv - 0.4947255151999997`\ k\^0.36`\ Log[A]\)}, {\(LP - 3.003003003003003`\ y\)}, {\(\(-0.09999999999999998`\) + 0.09999999999999998`\ Log[A]\)} }], "\[NoBreak]", ")"}], MatrixForm[ { Plus[ Power[ c, -1], Times[ -1, \[Lambda]]], Plus[ Times[ -1, bcap], Times[ .95082381299699636, Power[ k, .35999999999999999], \[Lambda], Log[ A]]], Plus[ .035101010101010188, Times[ -.17810118547199988, Power[ Power[ k, -1], .64000000000000001], Log[ A]]], Plus[ c, Times[ .025000000000000022, k], Times[ -.49472551519999969, Power[ k, .35999999999999999], Log[ A]]], Plus[ y, Times[ -.49472551519999969, Power[ k, .35999999999999999], Log[ A]]], Plus[ c, iv, Times[ -.49472551519999969, Power[ k, .35999999999999999], Log[ A]]], Plus[ LP, Times[ -3.0030030030030028, y]], Plus[ -.099999999999999978, Times[ .099999999999999978, Log[ A]]]}]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(sstate = FindRoot[fss \[Equal] ceros, {c, 0.4}, {bcap, 0.05}, {k, 0.2}, {\[Lambda], 0.5}, {A, 1.5}, {y, 2}, {iv, 0.25}, {LP, 3}, MaxIterations -> 5000]\)], "Input"], Cell[BoxData[ \({c \[Rule] 0.9171910485591523`, bcap \[Rule] 2.584627971546402`, k \[Rule] 12.650421426900639`, \[Lambda] \[Rule] 1.090285389909367`, A \[Rule] 2.718281828459045`, y \[Rule] 1.2334515842316687`, iv \[Rule] 0.3162605356725163`, LP \[Rule] 3.704058811506512`}\)], "Output"] }, Open ]], Cell[BoxData[ \(\(v1 = Table[sstate[\([i, 2]\)], {i, nvar}];\)\)], "Input"], Cell[BoxData[ \(\(v1[\([2]\)] = h;\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(MatrixForm[v1]\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0.9171910485591523`"}, {"0.333`"}, {"12.650421426900639`"}, {"1.090285389909367`"}, {"2.718281828459045`"}, {"1.2334515842316687`"}, {"0.3162605356725163`"}, {"3.704058811506512`"} }], "\[NoBreak]", ")"}], MatrixForm[ {.91719104855915234, .33300000000000002, 12.650421426900639, 1.090285389909367, 2.7182818284590451, 1.2334515842316687, .31626053567251627, 3.7040588115065121}]]], "Output"] }, Open ]], Cell[BoxData[ \(\(bcap = sstate[\([2, 2]\)];\)\)], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["Transforming Variables: I", "Section"], Cell[CellGroupData[{ Cell[BoxData[ \(transf = Table[If[v1[\([i]\)] =!= 0, v\&^[\([i]\)] \[Rule] v1[\([i]\)]\/Abs[v1[\([i]\)]]*v1[\([i]\)], v\&^[\([i]\)] \[Rule] v1[\([i]\)]], {i, nvar}]\)], "Input"], Cell[BoxData[ \({\(c\&^\) \[Rule] 0.9171910485591523`, \(0.333`\&^\) \[Rule] 0.333`, \(k\&^\) \[Rule] 12.650421426900639`, \(\[Lambda]\&^\) \[Rule] 1.090285389909367`, \(A\&^\) \[Rule] 2.718281828459045`, \(y\&^\) \[Rule] 1.2334515842316687`, \(iv\&^\) \[Rule] 0.3162605356725163`, \(LP\&^\) \[Rule] 3.7040588115065116`}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(v\&^ = Table[transf[\([i, 2]\)], {i, nvar}]\)], "Input"], Cell[BoxData[ \({0.9171910485591523`, 0.333`, 12.650421426900639`, 1.090285389909367`, 2.718281828459045`, 1.2334515842316687`, 0.3162605356725163`, 3.7040588115065116`}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(MatrixForm[%]\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0.9171910485591523`"}, {"0.333`"}, {"12.650421426900639`"}, {"1.090285389909367`"}, {"2.718281828459045`"}, {"1.2334515842316687`"}, {"0.3162605356725163`"}, {"3.7040588115065116`"} }], "\[NoBreak]", ")"}], MatrixForm[ {.91719104855915234, .33300000000000002, 12.650421426900639, 1.090285389909367, 2.7182818284590451, 1.2334515842316687, .31626053567251627, 3.7040588115065116}]]], "Output"] }, Open ]], Cell[BoxData[ \(\(Table[ If[v1[\([i]\)] =!= v\&^[\([i]\)], StylePrint[\ \ WARNING : \ cambi\[OAcute]\ de\ signo\ v[\([i]\)]\ , "\"]], {i, nvar}];\)\)], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["Transforming Variables: II", "Section"], Cell[CellGroupData[{ Cell[BoxData[ \(xhat0 = Table[If[v1[\([i]\)] =!= 0, If[v1[\([i]\)] < 0, vt[\([i]\)] \[Rule] \((\(-v\&^[\([i]\)]\)*Exp[vht[\([i]\)]])\), vt[\([i]\)] \[Rule] v\&^[\([i]\)]*Exp[vht[\([i]\)]]], vt[\([i]\)] \[Rule] v\&^[\([i]\)]], {i, nvar}]\)], "Input"], Cell[BoxData[ \({ct \[Rule] 0.9171910485591523`\ \[ExponentialE]\^cht, ht \[Rule] 0.333`\ \[ExponentialE]\^hht, kt \[Rule] 12.650421426900639`\ \[ExponentialE]\^kht, \[Lambda]t \[Rule] 1.090285389909367`\ \[ExponentialE]\^\[Lambda]ht, At \[Rule] 2.718281828459045`\ \[ExponentialE]\^Aht, yt \[Rule] 1.2334515842316687`\ \[ExponentialE]\^yht, ivt \[Rule] 0.3162605356725163`\ \[ExponentialE]\^ivht, LPt \[Rule] 3.7040588115065116`\ \[ExponentialE]\^LPht}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(xhat1 = Table[If[v1[\([i]\)] =!= 0, If[v1[\([i]\)] < 0, vt1[\([i]\)] \[Rule] \((\(-v\&^[\([i]\)]\)*Exp[vht1[\([i]\)]])\), vt1[\([i]\)] \[Rule] v\&^[\([i]\)]*Exp[vht1[\([i]\)]]], vt1[\([i]\)] \[Rule] v\&^[\([i]\)]], {i, nvar}]\)], "Input"], Cell[BoxData[ \({ct1 \[Rule] 0.9171910485591523`\ \[ExponentialE]\^cht1, ht1 \[Rule] 0.333`\ \[ExponentialE]\^hht1, kt1 \[Rule] 12.650421426900639`\ \[ExponentialE]\^kht1, \[Lambda]t1 \[Rule] 1.090285389909367`\ \[ExponentialE]\^\[Lambda]ht1, At1 \[Rule] 2.718281828459045`\ \[ExponentialE]\^Aht1, yt1 \[Rule] 1.2334515842316687`\ \[ExponentialE]\^yht1, ivt1 \[Rule] 0.3162605356725163`\ \[ExponentialE]\^ivht1, LPt1 \[Rule] 3.7040588115065116`\ \[ExponentialE]\^LPht1}\)], "Output"] }, Open ]], Cell[BoxData[ \(\(xhat = Flatten[{xhat0, xhat1}];\)\)], "Input"], Cell[BoxData[ \(\(fh = f /. xhat;\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(MatrixForm[%]\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {\(1.0902853899096978`\ \[ExponentialE]\^\(-cht\) - 1.090285389909367`\ \[ExponentialE]\^\[Lambda]ht\)}, {\(\(-2.584627971546402`\) + 2.584627971514956`\ \[ExponentialE]\^\[Lambda]ht\ \((\ \[ExponentialE]\^\(\(-hht\) + kht\))\)\^0.36`\ Log[ 2.718281828459045`\ \[ExponentialE]\^Aht]\)}, {\(\(-0.975`\) + 1.0101010101010102`\ \[ExponentialE]\^\(\[Lambda]ht - \ \[Lambda]ht1\) - 0.035101010103325926`\ \((\[ExponentialE]\^\(hht1 - kht1\))\)\ \^0.64`\ Log[2.718281828459045`\ \[ExponentialE]\^Aht1]\)}, {\(0.9171910485591523`\ \[ExponentialE]\^cht - 12.334160891228123`\ \[ExponentialE]\^kht + 12.650421426900639`\ \[ExponentialE]\^kht1 - 1.2334515842138059`\ \((\[ExponentialE]\^hht)\)\^0.64`\ \((\ \[ExponentialE]\^kht)\)\^0.36`\ Log[ 2.718281828459045`\ \[ExponentialE]\^Aht]\)}, {\(1.2334515842316687`\ \[ExponentialE]\^yht - 1.2334515842138059`\ \((\[ExponentialE]\^hht)\)\^0.64`\ \((\ \[ExponentialE]\^kht)\)\^0.36`\ Log[ 2.718281828459045`\ \[ExponentialE]\^Aht]\)}, {\(0.9171910485591523`\ \[ExponentialE]\^cht + 0.3162605356725163`\ \[ExponentialE]\^ivht - 1.2334515842138059`\ \((\[ExponentialE]\^hht)\)\^0.64`\ \((\ \[ExponentialE]\^kht)\)\^0.36`\ Log[ 2.718281828459045`\ \[ExponentialE]\^Aht]\)}, {\(3.7040588115065116`\ \[ExponentialE]\^LPht - 3.7040588115065125`\ \[ExponentialE]\^\(\(-hht\) + yht\)\)}, {\(\(-0.09999999999999998`\) - 0.9`\ Log[2.718281828459045`\ \[ExponentialE]\^Aht] + Log[2.718281828459045`\ \[ExponentialE]\^Aht1]\)} }], "\[NoBreak]", ")"}], MatrixForm[ { Plus[ Times[ 1.0902853899096978, Power[ E, Times[ -1, cht]]], Times[ -1.090285389909367, Power[ E, \[Lambda]ht]]], Plus[ -2.5846279715464018, Times[ 2.5846279715149558, Power[ E, \[Lambda]ht], Power[ Power[ E, Plus[ Times[ -1, hht], kht]], .35999999999999999], Log[ Times[ 2.7182818284590451, Power[ E, Aht]]]]], Plus[ -.97499999999999998, Times[ 1.0101010101010102, Power[ E, Plus[ \[Lambda]ht, Times[ -1, \[Lambda]ht1]]]], Times[ -.035101010103325926, Power[ Power[ E, Plus[ hht1, Times[ -1, kht1]]], .64000000000000001], Log[ Times[ 2.7182818284590451, Power[ E, Aht1]]]]], Plus[ Times[ .91719104855915234, Power[ E, cht]], Times[ -12.334160891228123, Power[ E, kht]], Times[ 12.650421426900639, Power[ E, kht1]], Times[ -1.2334515842138059, Power[ Power[ E, hht], .64000000000000001], Power[ Power[ E, kht], .35999999999999999], Log[ Times[ 2.7182818284590451, Power[ E, Aht]]]]], Plus[ Times[ 1.2334515842316687, Power[ E, yht]], Times[ -1.2334515842138059, Power[ Power[ E, hht], .64000000000000001], Power[ Power[ E, kht], .35999999999999999], Log[ Times[ 2.7182818284590451, Power[ E, Aht]]]]], Plus[ Times[ .91719104855915234, Power[ E, cht]], Times[ .31626053567251627, Power[ E, ivht]], Times[ -1.2334515842138059, Power[ Power[ E, hht], .64000000000000001], Power[ Power[ E, kht], .35999999999999999], Log[ Times[ 2.7182818284590451, Power[ E, Aht]]]]], Plus[ Times[ 3.7040588115065116, Power[ E, LPht]], Times[ -3.7040588115065125, Power[ E, Plus[ Times[ -1, hht], yht]]]], Plus[ -.099999999999999978, Times[ -.90000000000000002, Log[ Times[ 2.7182818284590451, Power[ E, Aht]]]], Log[ Times[ 2.7182818284590451, Power[ E, Aht1]]]]}]]], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Log Linear Approximation", "Section"], Cell[BoxData[ \(\(xhatss0 = Table[vht[\([i]\)] \[Rule] 0, {i, nvar}];\)\)], "Input"], Cell[BoxData[ \(\(xhatss1 = Table[vht1[\([i]\)] \[Rule] 0, {i, nvar}];\)\)], "Input"], Cell[BoxData[ \(\(xhatss = Flatten[{xhatss0, xhatss1}];\)\)], "Input"], Cell[BoxData[ \(\(fhss = fh /. xhatss;\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(MatrixForm[fhss]\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"3.3084646133829665`*^-13"}, {\(-3.1445956949482934`*^-11\)}, {\(-2.315737879232671`*^-12\)}, {"1.786215619858922`*^-11"}, {"1.7862822332403994`*^-11"}, {"1.7862822332403994`*^-11"}, {\(-8.881784197001252`*^-16\)}, {"0.`"} }], "\[NoBreak]", ")"}], MatrixForm[ {.33084646133829665*^-12, -.31445956949482934*^-10, \ -.2315737879232671*^-11, .17862156198589219*^-10, .17862822332403994*^-10, \ .17862822332403994*^-10, -.88817841970012523*^-15, 0.0}]]], "Output"] }, Open ]], Cell[BoxData[ \(\(aux1 = Transpose[Grad[fh, vht1] /. xhatss];\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(MatrixForm[aux1]\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0", "0", "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", "0", "0"}, {"0", \(-0.022464646466128592`\), "0.022464646466128592`", \(-1.0101010101010102`\), \ \(-0.035101010103325926`\), "0", "0", "0"}, {"0", "0", "12.650421426900639`", "0", "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", "0", "0"}, {"0", "0", "0", "0", "1.`", "0", "0", "0"} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[BoxData[ \(\(aux = Transpose[Grad[fh, vht] /. xhatss];\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(MatrixForm[aux]\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {\(-1.0902853899096978`\), "0", "0", \(-1.090285389909367`\), "0", "0", "0", "0"}, {"0", \(-0.930466069745384`\), "0.930466069745384`", "2.584627971514956`", "2.584627971514956`", "0", "0", "0"}, {"0", "0", "0", "1.0101010101010102`", "0", "0", "0", "0"}, { "0.9171910485591523`", \(-0.7894090138968357`\), \ \(-12.778203461545093`\), "0", \(-1.2334515842138059`\), "0", "0", "0"}, {"0", \(-0.7894090138968357`\), \(-0.4440425703169701`\), "0", \(-1.2334515842138059`\), "1.2334515842316687`", "0", "0"}, { "0.9171910485591523`", \(-0.7894090138968357`\), \ \(-0.4440425703169701`\), "0", \(-1.2334515842138059`\), "0", "0.3162605356725163`", "0"}, {"0", "3.7040588115065125`", "0", "0", "0", \(-3.7040588115065125`\), "0", "3.7040588115065116`"}, {"0", "0", "0", "0", \(-0.9`\), "0", "0", "0"} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Export Matrices", "Section"], Cell[BoxData[ \(\(Export["\", aux, "\"];\)\)], "Input"], Cell[BoxData[ \(\(Export["\", aux1, "\"];\)\)], "Input"], Cell[BoxData[ \(\(Export["\", v1, "\"];\)\)], "Input"], Cell[BoxData[ \(\(Export["\", dim, "\"];\)\)], "Input"] }, Open ]] }, Open ]] }, FrontEndVersion->"4.1 for Microsoft Windows", ScreenRectangle->{{0, 1024}, {0, 670}}, WindowSize->{1015, 626}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic} ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1727, 52, 26, 0, 115, "Title"], Cell[1756, 54, 54, 0, 64, "Subtitle"], Cell[1813, 56, 77, 1, 30, "Input"], Cell[CellGroupData[{ Cell[1915, 61, 69, 0, 59, "Section"], Cell[1987, 63, 114, 2, 30, "Input"], Cell[2104, 67, 44, 1, 30, "Input"], Cell[2151, 70, 59, 1, 30, "Input"], Cell[2213, 73, 134, 2, 30, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[2384, 80, 37, 0, 59, "Section"], Cell[2424, 82, 59, 1, 30, "Input"], Cell[2486, 85, 83, 1, 30, "Input"], Cell[2572, 88, 87, 1, 30, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[2696, 94, 55, 0, 59, "Section"], Cell[2754, 96, 73, 1, 30, "Input"], Cell[2830, 99, 138, 2, 42, "Input"], Cell[2971, 103, 170, 3, 45, "Input"], Cell[3144, 108, 139, 3, 30, "Input"], Cell[3286, 113, 49, 1, 30, "Input"], Cell[3338, 116, 55, 1, 30, "Input"], Cell[3396, 119, 55, 1, 42, "Input"], Cell[3454, 122, 124, 3, 30, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[3615, 130, 47, 0, 59, "Section"], Cell[3665, 132, 74, 1, 30, "Input"], Cell[CellGroupData[{ Cell[3764, 137, 46, 1, 30, "Input"], Cell[3813, 140, 2011, 55, 197, "Output"] }, Open ]], Cell[5839, 198, 85, 1, 30, "Input"], Cell[5927, 201, 105, 2, 30, "Input"], Cell[6035, 205, 105, 2, 30, "Input"], Cell[6143, 209, 114, 2, 30, "Input"], Cell[6260, 213, 76, 1, 30, "Input"], Cell[6339, 216, 114, 2, 31, "Input"], Cell[6456, 220, 67, 1, 30, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[6560, 226, 50, 0, 59, "Section"], Cell[6613, 228, 54, 1, 30, "Input"], Cell[6670, 231, 93, 1, 30, "Input"], Cell[6766, 234, 103, 2, 30, "Input"], Cell[6872, 238, 69, 1, 30, "Input"], Cell[6944, 241, 52, 1, 30, "Input"], Cell[CellGroupData[{ Cell[7021, 246, 48, 1, 30, "Input"], Cell[7072, 249, 1763, 46, 181, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8872, 300, 211, 4, 30, "Input"], Cell[9086, 306, 319, 5, 29, "Output"] }, Open ]], Cell[9420, 314, 79, 1, 30, "Input"], Cell[9502, 317, 53, 1, 30, "Input"], Cell[CellGroupData[{ Cell[9580, 322, 47, 1, 30, "Input"], Cell[9630, 325, 614, 15, 151, "Output"] }, Open ]], Cell[10259, 343, 63, 1, 30, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[10359, 349, 44, 0, 59, "Section"], Cell[CellGroupData[{ Cell[10428, 353, 210, 4, 44, "Input"], Cell[10641, 359, 401, 8, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11079, 372, 76, 1, 30, "Input"], Cell[11158, 375, 201, 3, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11396, 383, 46, 1, 30, "Input"], Cell[11445, 386, 615, 15, 151, "Output"] }, Open ]], Cell[12075, 404, 220, 5, 30, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[12332, 414, 45, 0, 59, "Section"], Cell[CellGroupData[{ Cell[12402, 418, 308, 6, 50, "Input"], Cell[12713, 426, 531, 9, 48, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13281, 440, 313, 6, 50, "Input"], Cell[13597, 448, 547, 9, 48, "Output"] }, Open ]], Cell[14159, 460, 68, 1, 30, "Input"], Cell[14230, 463, 52, 1, 30, "Input"], Cell[CellGroupData[{ Cell[14307, 468, 46, 1, 30, "Input"], Cell[14356, 471, 4742, 118, 179, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[19147, 595, 43, 0, 59, "Section"], Cell[19193, 597, 88, 1, 30, "Input"], Cell[19284, 600, 89, 1, 30, "Input"], Cell[19376, 603, 74, 1, 30, "Input"], Cell[19453, 606, 57, 1, 30, "Input"], Cell[CellGroupData[{ Cell[19535, 611, 49, 1, 30, "Input"], Cell[19587, 614, 655, 14, 171, "Output"] }, Open ]], Cell[20257, 631, 80, 1, 30, "Input"], Cell[CellGroupData[{ Cell[20362, 636, 49, 1, 30, "Input"], Cell[20414, 639, 735, 15, 151, "Output"] }, Open ]], Cell[21164, 657, 78, 1, 30, "Input"], Cell[CellGroupData[{ Cell[21267, 662, 48, 1, 30, "Input"], Cell[21318, 665, 1189, 23, 151, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[22556, 694, 34, 0, 59, "Section"], Cell[22593, 696, 78, 1, 30, "Input"], Cell[22674, 699, 79, 1, 30, "Input"], Cell[22756, 702, 79, 1, 30, "Input"], Cell[22838, 705, 77, 1, 30, "Input"] }, Open ]] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)