(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Wolfram 14.3' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 154, 7] NotebookDataLength[ 24335, 580] NotebookOptionsPosition[ 22033, 534] NotebookOutlinePosition[ 22429, 550] CellTagsIndexPosition[ 22386, 547] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{"GetBeta", ",", " ", "Cost"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"GetBeta", "[", RowBox[{"{", RowBox[{"cc_", ",", "epsilon_", ",", "k1_", ",", "k2_"}], "}"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"beta", ",", "a", ",", "b", ",", "c", ",", "d"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"d", "=", "cc"}], ";", "\[IndentingNewLine]", RowBox[{"c", "=", RowBox[{ RowBox[{"k2", "^", "2"}], "+", RowBox[{"epsilon", "^", "2"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"b", "=", RowBox[{"k1", "*", "k2"}]}], ";", "\[IndentingNewLine]", RowBox[{"a", "=", RowBox[{ RowBox[{"k1", "^", "2"}], "/", "3"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Return", "[", RowBox[{"{", RowBox[{"d", ",", "c", ",", "b", ",", "a"}], "}"}], "]"}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Cost", "[", RowBox[{"pars_", ",", "lam2_", ",", "x_", ",", "y_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"beta", ",", "X", ",", "nllf", ",", "s", ",", "i"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"beta", " ", "=", " ", RowBox[{"GetBeta", "[", "pars", "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"s", "=", RowBox[{ RowBox[{"2", "*", "y"}], "-", "1"}]}], ";", "\[IndentingNewLine]", RowBox[{"X", "=", RowBox[{"Transpose", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"x", "^", "i"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "3"}], "}"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"nllf", " ", "=", " ", RowBox[{"Total", "[", RowBox[{"Log", "[", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "s"}], "*", RowBox[{"(", RowBox[{"X", ".", "beta"}], ")"}]}], "]"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Return", "[", RowBox[{"nllf", " ", "+", " ", RowBox[{"lam2", " ", RowBox[{"Total", "[", RowBox[{"pars", "^", "2"}], "]"}]}]}], "]"}]}]}], "\[IndentingNewLine]", "]"}]}]}], "Input", CellChangeTimes->{{3.980215196159321*^9, 3.98021520112125*^9}, { 3.980215465357512*^9, 3.9802155187192497`*^9}, {3.980215561493145*^9, 3.98021556877605*^9}, {3.980215603281454*^9, 3.980215766491699*^9}, { 3.980215809381616*^9, 3.980215824734212*^9}, {3.9802159883378687`*^9, 3.980216092393737*^9}, {3.9802179711796503`*^9, 3.980217994933426*^9}, { 3.980218037526173*^9, 3.9802180508589983`*^9}, {3.980218114391039*^9, 3.9802182123715487`*^9}, {3.980218668097128*^9, 3.9802186800955877`*^9}, { 3.980219900715466*^9, 3.980219900917904*^9}, {3.9802199440804777`*^9, 3.980219950481944*^9}, {3.980220901657138*^9, 3.9802209090949917`*^9}},ExpressionUUID->"db312a24-878f-40c8-962f-\ 6aedbdf27fe4"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"TRIM", ",", " ", "log"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x", "=", RowBox[{"{", RowBox[{ "0.7444588568509963", ",", "0.4294306635860587", ",", "0.2424791244329582", ",", "0.1722166015253219", ",", "0.5155957608380367", ",", RowBox[{"-", "0.0044000248522457"}], ",", "0.2619898405306257", ",", "0.3707014579741952", ",", "0.3458779082346233", ",", RowBox[{"-", "0.0230997734892581"}], ",", "1.4235506263993014", ",", "0.1974151389300383", ",", "0.3479428293749495", ",", "0.038497594143739", ",", RowBox[{"-", "0.0494408800705622"}], ",", "0.1562561841152657", ",", "0.2440810426279506", ",", "0.738605475021583", ",", "0.3388466870335315", ",", "0.2730058878681247", ",", "0.2349466740981531", ",", "0.0152345391191637", ",", RowBox[{"-", "0.1306983645886834"}], ",", "0.3863126943346241", ",", "0.32696785342031", ",", "0.3960496407131636", ",", "0.116543550839155", ",", "0.3745678322136136", ",", "0.3902449371047537", ",", RowBox[{"-", "0.0224728544974076"}], ",", "0.3120049574170508", ",", "0.140909818171296", ",", "0.4065884615650837", ",", "0.1347990203891564", ",", "0.1747544412028135", ",", "0.8025435019583259", ",", "0.3424702625117166", ",", RowBox[{"-", "0.2576749555237702"}], ",", "0.2098741187455145", ",", "0.1827383448113199", ",", "0.4344451063845957", ",", "0.3628263770555691", ",", "0.3418422883292066", ",", "0.2690071065545839", ",", "0.3785734614926337", ",", "0.2619248908950089", ",", "0.232085349459071", ",", "0.1599028123673964", ",", "0.5800908847567773", ",", "0.208084169807583", ",", RowBox[{"-", "0.0239525163747711"}], ",", "0.2010508033108026", ",", "0.4280507666730073", ",", "0.4184257480672208", ",", "0.444582693795233", ",", "0.2107665896841349", ",", "0.3642384221371574"}], "}"}]}], "\[IndentingNewLine]", RowBox[{"y", "=", RowBox[{"{", RowBox[{ "1", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]}], "\[IndentingNewLine]", RowBox[{"p", "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "1.1319747348643707"}], " ", RowBox[{"10", "^", "1"}]}], ",", RowBox[{ RowBox[{"-", "6.5316161792581431"}], "*", " ", RowBox[{"10", "^", RowBox[{"-", "5"}]}]}], ",", "1.6944026199144295", ",", "3.9357124963564112"}], "}"}]}], "\[IndentingNewLine]", RowBox[{"lam2", "=", RowBox[{"10", "^", RowBox[{"-", "3"}]}]}]}]}]], "Input", CellChangeTimes->{{3.980219550729746*^9, 3.980219552713852*^9}, { 3.980219691095812*^9, 3.9802197156948*^9}, {3.98022297307031*^9, 3.980222979562756*^9}},ExpressionUUID->"f22ff9ef-193c-434b-91da-\ e9e8b3175259"], Cell[BoxData[ RowBox[{"{", RowBox[{ "0.7444588568509963`", ",", "0.4294306635860587`", ",", "0.2424791244329582`", ",", "0.1722166015253219`", ",", "0.5155957608380367`", ",", RowBox[{"-", "0.0044000248522457`"}], ",", "0.2619898405306257`", ",", "0.3707014579741952`", ",", "0.3458779082346233`", ",", RowBox[{"-", "0.0230997734892581`"}], ",", "1.4235506263993014`", ",", "0.1974151389300383`", ",", "0.3479428293749495`", ",", "0.038497594143739`", ",", RowBox[{"-", "0.0494408800705622`"}], ",", "0.1562561841152657`", ",", "0.2440810426279506`", ",", "0.738605475021583`", ",", "0.3388466870335315`", ",", "0.2730058878681247`", ",", "0.2349466740981531`", ",", "0.0152345391191637`", ",", RowBox[{"-", "0.1306983645886834`"}], ",", "0.3863126943346241`", ",", "0.32696785342031`", ",", "0.3960496407131636`", ",", "0.116543550839155`", ",", "0.3745678322136136`", ",", "0.3902449371047537`", ",", RowBox[{"-", "0.0224728544974076`"}], ",", "0.3120049574170508`", ",", "0.140909818171296`", ",", "0.4065884615650837`", ",", "0.1347990203891564`", ",", "0.1747544412028135`", ",", "0.8025435019583259`", ",", "0.3424702625117166`", ",", RowBox[{"-", "0.2576749555237702`"}], ",", "0.2098741187455145`", ",", "0.1827383448113199`", ",", "0.4344451063845957`", ",", "0.3628263770555691`", ",", "0.3418422883292066`", ",", "0.2690071065545839`", ",", "0.3785734614926337`", ",", "0.2619248908950089`", ",", "0.232085349459071`", ",", "0.1599028123673964`", ",", "0.5800908847567773`", ",", "0.208084169807583`", ",", RowBox[{"-", "0.0239525163747711`"}], ",", "0.2010508033108026`", ",", "0.4280507666730073`", ",", "0.4184257480672208`", ",", "0.444582693795233`", ",", "0.2107665896841349`", ",", "0.3642384221371574`"}], "}"}]], "Output", CellChangeTimes->{3.9802197180368643`*^9, 3.9802229692825193`*^9}, CellLabel->"Out[59]=",ExpressionUUID->"a7f74a5e-642c-4b21-9dc6-dc4532fcd2dd"], Cell[BoxData[ RowBox[{"{", RowBox[{ "1", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]], "Output", CellChangeTimes->{3.9802197180368643`*^9, 3.980222969285266*^9}, CellLabel->"Out[60]=",ExpressionUUID->"6fe63fdd-9529-4ff3-95a7-2b315610cc0a"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "11.319747348643707`"}], ",", RowBox[{"-", "0.00006531616179258143`"}], ",", "1.6944026199144295`", ",", "3.935712496356411`"}], "}"}]], "Output", CellChangeTimes->{3.9802197180368643`*^9, 3.98022296928763*^9}, CellLabel->"Out[61]=",ExpressionUUID->"9174505c-bb9d-4cfc-a204-c6c464fab952"], Cell[BoxData[ FractionBox["1", "1000"]], "Output", CellChangeTimes->{3.9802197180368643`*^9, 3.980222969288666*^9}, CellLabel->"Out[62]=",ExpressionUUID->"dfb9c680-45f0-495e-bc68-7c102a014c4a"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Cost", "[", RowBox[{"p", ",", "lam2", ",", "x", ",", "y"}], "]"}]], "Input", CellChangeTimes->{{3.9802197403954973`*^9, 3.980219744640079*^9}}, CellLabel->"In[17]:=",ExpressionUUID->"73502887-6057-4d13-881a-8932ae266804"], Cell[BoxData["2.9061537536536557`"], "Output", CellChangeTimes->{{3.980219744932208*^9, 3.9802197515000153`*^9}}, CellLabel->"Out[17]=",ExpressionUUID->"603f6210-ea83-425d-9d61-d60dffe5c701"] }, Open ]], Cell[BoxData[ RowBox[{"Needs", "[", "\"\\"", "]"}]], "Input", CellChangeTimes->{{3.9802198372526817`*^9, 3.980219843824862*^9}}, CellLabel->"In[18]:=",ExpressionUUID->"417b4445-9108-40f3-a80b-fda472fd9d54"], Cell[CellGroupData[{ Cell[BoxData["p"], "Input", CellChangeTimes->{3.9802203983793707`*^9}, CellLabel->"In[38]:=",ExpressionUUID->"3ca07280-6eb4-41d4-8ab2-dfe9bbf7b8be"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "11.319747348643707`"}], ",", RowBox[{"-", "0.00006531616179258143`"}], ",", "1.6944026199144295`", ",", "3.935712496356411`"}], "}"}]], "Output", CellChangeTimes->{3.980220398702986*^9}, CellLabel->"Out[38]=",ExpressionUUID->"029d324b-6436-44a3-8803-14c8381b5896"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", "gradient", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"DCost", "[", RowBox[{"i_", ",", "p_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"u", ",", "ii"}], "}"}], ",", RowBox[{"ND", "[", RowBox[{ RowBox[{"Cost", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "==", "ii"}], ",", "u", ",", RowBox[{"p", "[", RowBox[{"[", "ii", "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"ii", ",", RowBox[{"Length", "[", "p", "]"}]}], "}"}]}], "]"}], ",", "lam2", ",", "x", ",", "y"}], "]"}], ",", "u", ",", RowBox[{"p", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"Terms", "->", "10"}]}], "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"DCost", "[", RowBox[{"i", ",", "p"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "p", "]"}]}], "}"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.98021979703826*^9, 3.980219827377369*^9}, { 3.980219863431608*^9, 3.9802198819049683`*^9}, {3.980219995402974*^9, 3.980219997759408*^9}, {3.980220196481331*^9, 3.980220213434679*^9}, { 3.980220348281417*^9, 3.980220363235091*^9}, {3.980220425600909*^9, 3.9802205459722*^9}, {3.980220589175432*^9, 3.980220594692799*^9}, { 3.9802207059692497`*^9, 3.9802207262248*^9}, {3.9802210428429956`*^9, 3.980221072838937*^9}, {3.980221280858418*^9, 3.980221289567711*^9}, { 3.980221422741997*^9, 3.98022142313131*^9}},ExpressionUUID->"b0ded771-dccd-45f6-a46c-\ 9b2ff466801f"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "2.02166414330068`*^-6"}], ",", RowBox[{"-", "3.137769860503177`*^-8"}], ",", "4.938462760785317`*^-7", ",", RowBox[{"-", "2.546636364252448`*^-6"}]}], "}"}]], "Output", CellChangeTimes->{{3.980219870760887*^9, 3.9802199094790907`*^9}, 3.980219962691937*^9, {3.980219998174837*^9, 3.9802200191940937`*^9}, { 3.980220197643147*^9, 3.980220214115034*^9}, 3.980220363622348*^9, { 3.98022042632763*^9, 3.9802204383576603`*^9}, 3.980220547299224*^9, { 3.980220712553927*^9, 3.9802207272890863`*^9}, {3.980221066313427*^9, 3.980221073167602*^9}, 3.980221290310155*^9}, CellLabel->"Out[49]=",ExpressionUUID->"2d1c80cf-dfd3-4a34-8f1d-5ca8d4f13801"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", "Hessian", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"HCost", "[", RowBox[{"i_", ",", "j_", ",", "p"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"u", ",", "v", ",", "k", ",", "F", ",", "G"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"F", "[", "v_", "]"}], ":=", RowBox[{"ND", "[", RowBox[{ RowBox[{"Cost", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "==", "k"}], ",", "u", ",", RowBox[{"If", "[", RowBox[{ RowBox[{"j", "==", "k"}], ",", "v", ",", RowBox[{"p", "[", RowBox[{"[", "k", "]"}], "]"}]}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", RowBox[{"Length", "[", "p", "]"}]}], "}"}]}], "]"}], ",", "lam2", ",", "x", ",", "y"}], "]"}], ",", "u", ",", RowBox[{"p", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"Terms", "->", "10"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"G", ":=", RowBox[{"ND", "[", RowBox[{ RowBox[{"Cost", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "==", "k"}], ",", "u", ",", RowBox[{"p", "[", RowBox[{"[", "k", "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", RowBox[{"Length", "[", "p", "]"}]}], "}"}]}], "]"}], ",", "lam2", ",", "x", ",", "y"}], "]"}], ",", RowBox[{"{", RowBox[{"u", ",", "2"}], "}"}], ",", RowBox[{"p", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"Terms", "->", "10"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "!=", "j"}], ",", "\[IndentingNewLine]", RowBox[{"ND", "[", RowBox[{ RowBox[{"F", "[", "v", "]"}], ",", "v", ",", RowBox[{"p", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", "\[IndentingNewLine]", "G"}], "]"}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"H0", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"HCost", "[", RowBox[{"i", ",", "j", ",", "p"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "p", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", "p", "]"}]}], "}"}]}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.980220845242152*^9, 3.9802210388102617`*^9}, { 3.980221083420125*^9, 3.980221210743389*^9}, {3.9802212959394293`*^9, 3.980221347003302*^9}, {3.9802214368899813`*^9, 3.980221498114513*^9}, { 3.980221552212894*^9, 3.980221718470955*^9}, {3.980221751831498*^9, 3.980221754488792*^9}, {3.980222174396463*^9, 3.9802222588377943`*^9}, { 3.9802223438606653`*^9, 3.9802223486758127`*^9}},ExpressionUUID->"db853e79-d0e4-4a86-8783-\ 98f6bc5933c4"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0.7673411749127825`", ",", RowBox[{"-", "0.00004960801306647878`"}], ",", "0.889749135171745`", ",", "3.3212700921235605`"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.0000496078733072483`"}], ",", "0.0004804061483556116`", ",", RowBox[{"-", "0.0000624478957252519`"}], ",", RowBox[{"-", "0.00022472631280071956`"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0.8897491344529453`", ",", RowBox[{"-", "0.00006244763646673682`"}], ",", "1.2193302465332434`", ",", "4.2116021632604665`"}], "}"}], ",", RowBox[{"{", RowBox[{"3.3212689431121136`", ",", RowBox[{"-", "0.0002247255362663118`"}], ",", "4.211598490108457`", ",", "15.106151994325286`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.980221740709745*^9, 3.9802217766213427`*^9, 3.9802222890883713`*^9}, CellLabel->"Out[56]=",ExpressionUUID->"9ac024a3-d751-44de-bc53-6ad3de884f93"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"MatrixForm", "[", "H0", "]"}]], "Input", CellChangeTimes->{{3.980221756962586*^9, 3.980221761796362*^9}, { 3.9802218237372227`*^9, 3.98022182432813*^9}}, CellLabel->"In[57]:=",ExpressionUUID->"c25363e3-6dd7-423a-81cd-c2c937573e5b"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0.7673411749127825`", RowBox[{"-", "0.00004960801306647878`"}], "0.889749135171745`", "3.3212700921235605`"}, { RowBox[{"-", "0.0000496078733072483`"}], "0.0004804061483556116`", RowBox[{"-", "0.0000624478957252519`"}], RowBox[{"-", "0.00022472631280071956`"}]}, {"0.8897491344529453`", RowBox[{"-", "0.00006244763646673682`"}], "1.2193302465332434`", "4.2116021632604665`"}, {"3.3212689431121136`", RowBox[{"-", "0.0002247255362663118`"}], "4.211598490108457`", "15.106151994325286`"} }, GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.980221824708803*^9, 3.9802222891345053`*^9}, CellLabel-> "Out[57]//MatrixForm=",ExpressionUUID->"60a8e5e9-338f-4169-8f24-\ 577a50a6f32e"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"Inverse", "[", "H0", "]"}], "]"}]], "Input", CellChangeTimes->{{3.980222302130584*^9, 3.980222304864068*^9}}, CellLabel->"In[58]:=",ExpressionUUID->"3c1e775c-af2e-4c18-8894-85a24d299c82"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"124.2107378939491`", "0.008819362051593838`", "99.68227366780131`", RowBox[{"-", "55.10069583812897`"}]}, {"0.008763138449467222`", "2081.5865173231978`", RowBox[{"-", "0.002461036440763724`"}], "0.029726128780833826`"}, {"99.68361653826865`", RowBox[{"-", "0.002420312442278083`"}], "102.15456747369905`", RowBox[{"-", "50.39738859190163`"}]}, { RowBox[{"-", "55.10103654542312`"}], "0.029702306550748403`", RowBox[{"-", "50.39706092462403`"}], "26.23155138113158`"} }, GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.98022230519285*^9}, CellLabel-> "Out[58]//MatrixForm=",ExpressionUUID->"635785cc-da30-4ada-b0a8-\ 311ff9635453"] }, Open ]] }, WindowSize->{1366.5, 683.25}, WindowMargins->{{7.5, Automatic}, {Automatic, 24}}, FrontEndVersion->"14.3 for Linux x86 (64-bit) (July 8, 2025)", StyleDefinitions->"Default.nb", ExpressionUUID->"9ea52935-3f71-45ad-8b9c-8ebfe3a2f1d5" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[554, 20, 3290, 81, 413, "Input",ExpressionUUID->"db312a24-878f-40c8-962f-6aedbdf27fe4"], Cell[CellGroupData[{ Cell[3869, 105, 3458, 67, 258, "Input",ExpressionUUID->"f22ff9ef-193c-434b-91da-e9e8b3175259"], Cell[7330, 174, 2011, 34, 97, "Output",ExpressionUUID->"a7f74a5e-642c-4b21-9dc6-dc4532fcd2dd"], Cell[9344, 210, 799, 12, 34, "Output",ExpressionUUID->"6fe63fdd-9529-4ff3-95a7-2b315610cc0a"], Cell[10146, 224, 354, 7, 34, "Output",ExpressionUUID->"9174505c-bb9d-4cfc-a204-c6c464fab952"], Cell[10503, 233, 196, 3, 60, "Output",ExpressionUUID->"dfb9c680-45f0-495e-bc68-7c102a014c4a"] }, Open ]], Cell[CellGroupData[{ Cell[10736, 241, 250, 4, 30, "Input",ExpressionUUID->"73502887-6057-4d13-881a-8932ae266804"], Cell[10989, 247, 193, 2, 33, "Output",ExpressionUUID->"603f6210-ea83-425d-9d61-d60dffe5c701"] }, Open ]], Cell[11197, 252, 230, 3, 30, "Input",ExpressionUUID->"417b4445-9108-40f3-a80b-fda472fd9d54"], Cell[CellGroupData[{ Cell[11452, 259, 150, 2, 29, "Input",ExpressionUUID->"3ca07280-6eb4-41d4-8ab2-dfe9bbf7b8be"], Cell[11605, 263, 331, 7, 34, "Output",ExpressionUUID->"029d324b-6436-44a3-8803-14c8381b5896"] }, Open ]], Cell[CellGroupData[{ Cell[11973, 275, 1869, 46, 73, "Input",ExpressionUUID->"b0ded771-dccd-45f6-a46c-9b2ff466801f"], Cell[13845, 323, 723, 12, 45, "Output",ExpressionUUID->"2d1c80cf-dfd3-4a34-8f1d-5ca8d4f13801"] }, Open ]], Cell[CellGroupData[{ Cell[14605, 340, 3534, 88, 247, "Input",ExpressionUUID->"db853e79-d0e4-4a86-8783-98f6bc5933c4"], Cell[18142, 430, 995, 22, 57, "Output",ExpressionUUID->"9ac024a3-d751-44de-bc53-6ad3de884f93"] }, Open ]], Cell[CellGroupData[{ Cell[19174, 457, 259, 4, 30, "Input",ExpressionUUID->"c25363e3-6dd7-423a-81cd-c2c937573e5b"], Cell[19436, 463, 1181, 30, 111, "Output",ExpressionUUID->"60a8e5e9-338f-4169-8f24-577a50a6f32e"] }, Open ]], Cell[CellGroupData[{ Cell[20654, 498, 243, 4, 30, "Input",ExpressionUUID->"3c1e775c-af2e-4c18-8894-85a24d299c82"], Cell[20900, 504, 1117, 27, 111, "Output",ExpressionUUID->"635785cc-da30-4ada-b0a8-311ff9635453"] }, Open ]] } ] *)