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

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

(* CreatedBy='Mathematica 11.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     27459,        644]
NotebookOptionsPosition[     26611,        615]
NotebookOutlinePosition[     27235,        637]
CellTagsIndexPosition[     27192,        634]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "Resolu\[CCedilla]\[ATilde]o", " ", "do", " ", "Problema", " ", "7"}], " ",
    "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", "n\[OAcute]s", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"xn", "=", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0.25", ",", "0.5", ",", "0.75", ",", "1"}], "}"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"nn", "=", 
     RowBox[{"Length", "[", "xn", "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"valores", " ", "nodais", " ", "tentativa"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"un", "=", 
     RowBox[{"{", 
      RowBox[{"un1", ",", "un2", ",", "un3", ",", "un4", ",", "un5"}], 
      "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"valores", " ", "nodais", " ", "teste"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"wn", "=", 
     RowBox[{"{", 
      RowBox[{"wn1", ",", "wn2", ",", "wn3", ",", "wn4", ",", "wn5"}], 
      "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "fun\[CCedilla]\[OTilde]es", " ", "de", " ", "forma", " ", "lineares"}], 
    " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ll", "[", "xi_", "]"}], ":=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "-", "xi"}], ")"}], "/", "2"}], ",", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "+", "xi"}], ")"}], "/", "2"}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "transforma\[CCedilla]\[ATilde]o", " ", "de", " ", "coordenadas"}], " ", 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"xid", "[", 
      RowBox[{"x1_", ",", "x2_", ",", "x_"}], "]"}], ":=", 
     RowBox[{
      RowBox[{"2", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"x", "-", "x1"}], ")"}], "/", 
        RowBox[{"(", 
         RowBox[{"x2", "-", "x1"}], ")"}]}]}], "-", "1"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "fun\[CCedilla]\[ATilde]o", " ", "de", " ", 
     "interpola\[CCedilla]\[ATilde]o", " ", "por", " ", "tro\[CCedilla]os"}], 
    " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"basic", "[", 
      RowBox[{"i_", ",", "x_", ",", "un_"}], "]"}], ":=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{"ll", "[", 
         RowBox[{"xid", "[", 
          RowBox[{
           RowBox[{"xn", "[", 
            RowBox[{"[", "i", "]"}], "]"}], ",", 
           RowBox[{"xn", "[", 
            RowBox[{"[", 
             RowBox[{"i", "+", "1"}], "]"}], "]"}], ",", "x"}], "]"}], "]"}], 
        ".", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"un", "[", 
           RowBox[{"[", "i", "]"}], "]"}], ",", 
          RowBox[{"un", "[", 
           RowBox[{"[", 
            RowBox[{"i", "+", "1"}], "]"}], "]"}]}], "}"}]}], ",", 
       RowBox[{
        RowBox[{"x", ">", 
         RowBox[{"xn", "[", 
          RowBox[{"[", "i", "]"}], "]"}]}], "&&", 
        RowBox[{"x", "\[LessEqual]", 
         RowBox[{"xn", "[", 
          RowBox[{"[", 
           RowBox[{"i", "+", "1"}], "]"}], "]"}]}]}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"joined", "[", 
      RowBox[{"x_", ",", "un_"}], "]"}], ":=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"basic", "[", 
        RowBox[{"i", ",", "x", ",", "un"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"nn", "-", "1"}]}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"interpl", "[", 
      RowBox[{"x_", ",", "un_"}], "]"}], ":=", 
     RowBox[{"Piecewise", "[", 
      RowBox[{"joined", "[", 
       RowBox[{"x", ",", "un"}], "]"}], "]"}]}], ";"}], "\n", 
   RowBox[{"(*", " ", 
    RowBox[{
    "fun\[CCedilla]\[OTilde]es", " ", "de", " ", "tentativa", " ", "e", " ", 
     "de", " ", "teste"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"u", "[", "x_", "]"}], ":=", 
     RowBox[{"interpl", "[", 
      RowBox[{"x", ",", "un"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"w", "[", "x_", "]"}], ":=", 
     RowBox[{"interpl", "[", 
      RowBox[{"x", ",", "wn"}], "]"}]}], ";"}], "\n", 
   RowBox[{"(*", " ", 
    RowBox[{"integrando", " ", "da", " ", "forma", " ", "fraca"}], " ", 
    "*)"}], "\n", 
   RowBox[{
    RowBox[{"integrando", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"-", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"u", "[", "x", "]"}], ",", "x"}], "]"}]}], "*", 
       RowBox[{"D", "[", 
        RowBox[{
         RowBox[{"w", "[", "x", "]"}], ",", "x"}], "]"}]}], "-", 
      RowBox[{
       RowBox[{"w", "[", "x", "]"}], "*", 
       RowBox[{"u", "[", "x", "]"}]}], "+", 
      RowBox[{"w", "[", "x", "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"integral", " ", "da", " ", "forma", " ", "fraca"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"eq", "=", 
     RowBox[{"Integrate", "[", 
      RowBox[{"integrando", ",", 
       RowBox[{"{", 
        RowBox[{"x", ",", 
         RowBox[{"xn", "[", 
          RowBox[{"[", "1", "]"}], "]"}], ",", 
         RowBox[{"xn", "[", 
          RowBox[{"[", "nn", "]"}], "]"}]}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", "for\[CCedilla]as", " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f", "=", 
     RowBox[{"D", "[", 
      RowBox[{"eq", ",", 
       RowBox[{"{", "wn", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"for\[CCedilla]as", " ", "na", " ", "origem"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"rulezero", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"un", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "\[Rule]", "0"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "nn"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f0", "=", 
     RowBox[{"f", "/.", "rulezero"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"f0", "//", "MatrixForm"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "rigidez", " ", "\[EAcute]", " ", "a", " ", "derivada", " ", "das", " ", 
     "for\[CCedilla]as"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"k", "=", 
     RowBox[{"-", 
      RowBox[{"D", "[", 
       RowBox[{"f", ",", 
        RowBox[{"{", "un", "}"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"k", "//", "MatrixForm"}], "\n", 
   RowBox[{"(*", " ", 
    RowBox[{
    "Matriz", " ", "das", " ", "condi\[CCedilla]\[OTilde]es", " ", "de", " ", 
     "fronteira", " ", "essenciais"}], " ", "*)"}], "\n", 
   RowBox[{
    RowBox[{"transf", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1", ",", "0", ",", "0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "0", ",", "1", ",", "0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "0", ",", "0", ",", "1", ",", "0"}], "}"}]}], 
      "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "Transforma\[CCedilla]\[ATilde]o", " ", "das", " ", "for\[CCedilla]as", 
     " ", "e", " ", "rigidez"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"fred", "=", 
     RowBox[{"transf", ".", "f0"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"kred", "=", 
     RowBox[{"transf", ".", "k", ".", 
      RowBox[{"Transpose", "[", "transf", "]"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "Solu\[CCedilla]\[ATilde]o", " ", "do", " ", "problema", " ", "linear"}], 
    " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ured", "=", 
     RowBox[{
      RowBox[{"Inverse", "[", "kred", "]"}], ".", "fred"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"unod", "=", 
     RowBox[{
      RowBox[{"Transpose", "[", "transf", "]"}], ".", "ured"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", "Interpola\[CCedilla]\[ATilde]o", " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"usol", "[", "x_", "]"}], ":=", 
     RowBox[{"interpl", "[", 
      RowBox[{"x", ",", "unod"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", "Resultados", " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"usol", "[", "x", "]"}], ",", 
       RowBox[{"-", 
        FractionBox[
         RowBox[{
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"-", "x"}]], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", 
            SuperscriptBox["\[ExponentialE]", "x"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "\[ExponentialE]"}], "+", 
            SuperscriptBox["\[ExponentialE]", "x"]}], ")"}]}], 
         RowBox[{"1", "+", "\[ExponentialE]"}]]}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "1"}], "}"}]}], "]"}], "\n", "\n", 
   "\n"}]}]], "Input",
 CellChangeTimes->{{3.8244615993848624`*^9, 3.824461615733876*^9}, {
   3.8244616542877283`*^9, 3.8244618111165605`*^9}, {3.824462280541755*^9, 
   3.8244623214024544`*^9}, {3.824462356198717*^9, 3.824462467607777*^9}, {
   3.8244625115457916`*^9, 3.8244625975848417`*^9}, {3.8244626397347374`*^9, 
   3.824462828101886*^9}, {3.8244628829505644`*^9, 3.8244629522227445`*^9}, {
   3.824463100034727*^9, 3.824463126193757*^9}, {3.8244631671446724`*^9, 
   3.824463407334673*^9}, {3.8244634600622406`*^9, 3.8244634859955797`*^9}, {
   3.8244635306816487`*^9, 3.8244635801797905`*^9}, {3.8244636343582926`*^9, 
   3.8244637041573677`*^9}, {3.824463735170786*^9, 3.824463739013751*^9}, {
   3.824463891273753*^9, 3.8244639662667794`*^9}, {3.824464263168541*^9, 
   3.824464269799592*^9}, {3.8244643883971157`*^9, 3.8244643926086836`*^9}, {
   3.8244644928791103`*^9, 3.824464531981776*^9}, {3.8244657409433355`*^9, 
   3.8244658101337585`*^9}, {3.824465841487799*^9, 3.824466051233199*^9}, {
   3.824466333488858*^9, 3.824466345764758*^9}, {3.8244673879777665`*^9, 
   3.8244674453807645`*^9}, {3.8244674769334364`*^9, 3.824467479207781*^9}, {
   3.8244682468876867`*^9, 3.824468249393775*^9}, {3.8244682840015097`*^9, 
   3.8244684581687064`*^9}, {3.8244684968138385`*^9, 
   3.8244685130257854`*^9}, {3.8244685724017477`*^9, 3.82446858515777*^9}, {
   3.824468626212916*^9, 3.8244686837948236`*^9}, {3.824468809586833*^9, 
   3.8244688139767504`*^9}, {3.824468994488391*^9, 3.8244690880013266`*^9}, {
   3.82446923028374*^9, 3.824469320767507*^9}, {3.824469412207715*^9, 
   3.8244694442285156`*^9}, {3.824469515076744*^9, 3.824469582739748*^9}, {
   3.824469699483801*^9, 3.8244697295514736`*^9}, {3.8244697690216713`*^9, 
   3.8244697819258366`*^9}, {3.8244698802407722`*^9, 
   3.8244699134826775`*^9}, {3.8244701163157787`*^9, 3.824470124290844*^9}, {
   3.824470328362791*^9, 3.8244703351797266`*^9}, {3.824470883490892*^9, 
   3.824471004051675*^9}, {3.824471124840151*^9, 3.8244711290897007`*^9}, 
   3.8244858011482553`*^9, {3.824813076056366*^9, 3.8248131287838235`*^9}},
 ExpressionUUID -> "8c94ab20-d891-4bb7-b431-e14bf4d186c6"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", 
   TagBox[GridBox[{
      {"0.125`"},
      {"0.25`"},
      {"0.25`"},
      {"0.25`"},
      {"0.125`"}
     },
     GridBoxAlignment->{
      "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
        "RowsIndexed" -> {}},
     GridBoxSpacings->{"Columns" -> {
         Offset[0.27999999999999997`], {
          Offset[0.5599999999999999]}, 
         Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
         Offset[0.2], {
          Offset[0.4]}, 
         Offset[0.2]}, "RowsIndexed" -> {}}],
    Column], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.8248131322540407`*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"4.083333333333333`", 
      RowBox[{"-", "3.958333333333333`"}], "0", "0", "0"},
     {
      RowBox[{"-", "3.958333333333333`"}], "8.166666666666666`", 
      RowBox[{"-", "3.958333333333333`"}], "0", "0"},
     {"0", 
      RowBox[{"-", "3.958333333333333`"}], "8.166666666666666`", 
      RowBox[{"-", "3.958333333333333`"}], "0"},
     {"0", "0", 
      RowBox[{"-", "3.958333333333333`"}], "8.166666666666666`", 
      RowBox[{"-", "3.958333333333333`"}]},
     {"0", "0", "0", 
      RowBox[{"-", "3.958333333333333`"}], "4.083333333333333`"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.8248131322540407`*^9}],

Cell[BoxData[
 GraphicsBox[{{{{}, {}, 
     TagBox[
      {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
       Opacity[1.], LineBox[CompressedData["
1:eJwBoQJe/SFib1JlAgAAACkAAAACAAAASmXo0cDpVT4b718I4w4+PgRDGnDf
GTQ/oMadQ5WSGz8zn5icsxlEP+C11yVZkis/S83Xsp0ZVD/YrfQWO5I7P1dk
972SGWQ/1CmDDyySSz/dL4dDjRl0P7xnyoskkls/oBVPhooZhD+wBu7JIJJr
P4EIsyeJGZQ/Ktb/6B6Sez9YWaGv8PGkPwBo/lT2uow/C/+BXEcerz+ig+79
qVeVPyVR7ujDi7Q/ZtDmIs8unD93zjk/mPS5Pxr5DDswzaE/F7QR3gEBvz+Y
kj6WlkOlP7biQ8zPPMI/nZaNCDMEqT/Cap0bWezEPzEZoaNKs6w/9SY9D61t
xz8tUib4fxGwPwL5qxAbJco/Bs0KqnXusT82/2C2U67MPylMuzS6q7M/TE3A
gkW4zD8jlzc9jLKzP2ObH083wsw/HeKzRV65sz+QN97nGtbMPxB4rFYCx7M/
6W9bGeL9zD/3o514SuKzP5zgVXxwTc0/xPt/vNoYtD8BwkpCjezNP1+rRET7
hbQ/GBCqDn/2zT9Z9sBMzYy0Py5eCdtwAM4/U0E9VZ+TtD9a+sdzVBTOP0XX
NWZDobQ/tDJFpRs8zj8tAyeIi7y0P2ajPwiqi84/+loJzBvztD/MhDTOxirP
P5UKzlM8YLU/Og/LApE1zz9kkHjIome1P6iZYTdbQM8/MxYjPQlvtT+Dro6g
71XPP9AheCbWfbU/OtjochiBzz8LOSL5b5u1P6hif6fii88/2b7MbdaitT8W
7RXcrJbPP6hEd+I8qrU/8QFDRUGszz9FUMzLCbm1P6grnRdq188/gGd2nqPW
tT8WtjNMNOLPP0/tIBMK3rU/hEDKgP7szz8ec8uHcOW1P31orA+N9c8/94En
7E7rtT8x8jnp
        "]], LineBox[CompressedData["
1:eJwVj3841HcAx8+P2IgV53L3YJ02O0+eUkxRz97XacxYWy1KhVxzZGuZ2IjN
U1PNaS7VmtJqT6aFjHT3dficU/PwJDzfm/Nb389wnpiZx+oJ+TX74/W8/nn9
8xLLj++JM+fxeP7L/O8Hjz3T/Few+NHzwQn2uQZub18scVzJwvyLNxcGpzXI
1Uu4JwIWV7/2HmfmNXhomTIjc2FhsvCTi5c0mAnRO9wRs9io3Naba8ZAYQwP
TvZi0XQlqPmwFQPZ2OlKKxmL5xWHimxWMUjf0N7yaTAL+MlFqQ4MKk8IR9gw
FrkkPp/yGbgtVggL97EQP0rOUgsZzDo8Ob3xGIsPh747GL1uud++JXx/AYtS
vsaxyo/B6xfGZ0v+YRGbbpcdGcWgPjarpOoZi3c8Lhdci2EQ7eMYWTfDQtQh
Ku+PZXC9K0Dbam6AwUvSFaVg4Oyq/HJSYIDsT5lEfpyBQ6nkmS8McA9Oa008
xWBFY9zfepUBQ/xhfmYxg/FpynV6/wF5pbb4xgSDX/6qa7HP60Cc7mQ3T1kN
e9dfJ00mI7yyUmfPrtPi/FGN4uc3uiAaqclQV2nxQlpmSsjqho3YyX7pgxqk
bG5UC37vwVqzWt4UVwP6W2OuQdyHrdWpFrzMWswlGgs9T/bDyT3xxWZBHdbv
sF5U1A1AuOlOiN3tOgTfEsffW8PBVbBw11JG4Ftg79HkwiHwgq0kdCfB2rw5
U5+YQ7yNqCg/iOBleqfc0otD5aLfVbdQgvLdOdH7pMv906QzW/cSOJpNhS8k
cEjQmqI+UxDQmIadIbUc1Ada7Y1KgtbwcotoPYceY1+O6HsCbei1h8mNHBZ2
jVnGqggubkmRXm/nEBRoNTdxieDd1yTbJ4c49K6XPrX+iaCsXuXzgy3F0vx9
/bYKgq9cY9wHD1G8f7MwqaiDQP6Rd8RALMWO+cQMrpNg17dmym4FhX9kwFlh
D4HHWNFUWxKFp2NvoWqAwKgebajNpnjlnKA5Y4TAOyw55vJdiqbP8133ThO4
ZAVeUt2jqH98WKKaJbCu4jcrGQqNxNunZY6AW1O94VQDRfFQ+3tSng7nh+cW
jhkpsiNWpni9qkOaoG3T0V6KzPv938Tb6nAk5EbcJxzFiVVlObfsdAiokLYd
GKU40hJy09lBB4+h1WYRExQH3xKWfszXYbXTsO/ufyn2ZI+q8wQ6zAerE8Km
l38Hq/WPnHWAb9e5ppcUDbevRO4X6fAf88fB0w==
        "]], LineBox[CompressedData["
1:eJwVkX841HcAx/0sP+bIj6V75OcspzYboYXenjHp6BnzI9EjuWcbYt2K9Xg2
LSRzM8ku/XDMI3nk9Nx65j5d3fcym+nJ/G6ixfZF7kOULZxDp9kfr+f17+t5
Xi6pxz762EBPTy9wnf89NmnLSAxYtDZUHkzgqhDhNunLNWHRZpx8MctOhZ+Z
7t6LZiw6G8tKvW1U8Isn6XYWLPojWvOXrFRwLSmutrJmMX7eOTP/NRUuuQr9
y2xZPN0ZfSTMTAUOk9BvvpnFi6H8eHMTFZaf84w3OrAwdJwIvmCoQl9sV6bO
nYVb0zW76ysMdszGmQ/tZ9H4pyxN+5CBIjayVRTF4oBr7XTpHwxCVO+f2BPD
wjSjPMNlgEFimdfjawks0peFmfxuBiXvmEqzBSx49t5Cya8M1NkM3zqXhTSu
5WTwjwxqda6i/fUsZH2Ks9+UMLCx/MekfYUF6Wirtw9kYNRaJXOTjEHCEZ0+
L1Xi6E7lI9OIcRScDoh8xlFCui+19q5mHFeDWra559zB4gzXSl05AUn8RJ18
8jbkMXsNxaFPIL68asHj38bX4UU/tKif4JNdnfaLCgVcRvnDV4smke2a4lPt
qcBu/zPiLd5qnPL7O8ij/BYOOfYu6z9QQ9BgO2WudwucLcflykE1upzzDKp0
BG02dp/nDKvhUzXp4LlK8KZp0tTUqBr65SQ6XEMwN69+2DulRm1ugvLMLEHB
/bWfqtfUGI2UlOmGCRqz3/7sPU+KuHk337mbBInH+nnzOyhUWaUfnpIRWGRk
TzZ7UbwxtZBucYPg+OE7h5x9KeZHfqvZ3kgQwA+L2BhMUdHxqUlaDUGPU7LH
YBxFz2XpY1ZEsNj53bgwnyIs2KdgQEAwYj6Y9rKQYquY+0p0hOCXSIe54rMU
C1T/q5DDBOU9Tas131LUnevLaUkk8By4Z9t9gWLtr6z0ymiC5Ef6e3lNFPK8
xqhEEIRy93W3NFOU9p/rtg4i2J5UHhMsoxC4n+T/vptAO+KYcqCFYlP3B6GB
fgQVbEBu0V2KLIcJf8e3CDpojnRsgMJdudVpnEvQ7KF6N2uQ4iXHWHLFnuD7
dCOFdojiQeqsfczrBCkzFe2WoxT5Zkrr9k0EK89lo0GUYvRgwoYGEwLWa0lw
f3q97wYKkzcQ3BPueRo7SyHS36a32YhA/KJLc/Rfil3XF5eLX8nxpY9t3tI8
haVu5ItgnRyp2UmGhRoKdVT7vHZFjnB5XQlnef1HvVR4UyuHl2ba8soqhVhb
8SxdI4cTdyZEoKO4xGs70bcgx3/xbOjw
        "]], LineBox[CompressedData["
1:eJwVyn0s1HEAx3Gk3FjaYaMH2TmEPP7hKJtPsSGda2Sx4sjDajcUZZ7pLFIJ
qzwUWZ4uT3XX+G2tna0RyaY8nHB0nqIf7vpmnkrjV/3x3uufNyf6WnCcjpaW
lsO//ju7YCKv1qGxKT6mFqxQkJ0USM+zaHDMGtJqJylUZJa5TevT4K3NmkgU
FGJ3uN6bB2jENEcESnspVOoKbwZb0TDL4zsriinsxpi6DdnQcFV4jP658+/v
Htw4Z0cjwMY6kyum4JLnkxroRCOrf6c3+QaFPsY2w9+DxqyRVMgOo7C1vZbr
xafRVs8uFnAoXFi9e9/+Jo3TPV/Ufe0dMFZRdTrvaSSwLrd0jrdjb7Y45ETS
EnImFg8W7byGVrNQOme0DF/Jc5mBtwzJT6UPjLqW8TCzInWmqw2dxwtZ7MQV
7HOwi3MybgY7INkhz1ANesSyVlXXiCnHwZF0uRqn2ExPGF0LWQ2svgk1MCwq
SBYpqxEh5bHNdzVI+ZWv56hXjlv5teKplh9osSib2ocHeEmbVxnyCczCsj6k
K3NQ+Cw/SiUgKLeY+34gMgexwRrrV0EES4t+LMl8No7I5bLAUILiFOMzw+os
3Cu51FsUTeD6qKXfnsnAVd6TVf10AvNP4wOT3FRY3jbx15MQ1BSKmjKVidj1
yNo/1kTgY8pbejGfgAnN/LCkleCURNteoY5HaWh7hO9rgrnuylYHRgTGPvhG
vpzg827vq6/cK5gaKqnZM0Kgs+wrylUK8aZgK2ZklICvp1oZU4fjsWekXf04
QYVVSoILcxFnG506vFUEk5H112e5oXibNvBRvEQgamU0BUoByh3dSoLUBMyC
46o7E4CkueoQDiFY44Sv01w/2PLjp9+tE/x2SdooU3pBV1vRULpFwJUM/Nni
umOG8hRFbRNU8UJ0O5ROkIvqnV12CFSnD1n+VB5FhYXBBsMQHCaRrHlDL6+/
FZBoBQ==
        "]]},
      Annotation[#, "Charting`Private`Tag$2092#1"]& ], 
     TagBox[
      {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], 
       Opacity[1.], LineBox[CompressedData["
1:eJwV2Hk4Vd0XB3BECSlekmYSKU2ifkW+JSmhzEOvoSghhJQUKRFKlAoVIqIy
ZYi4RyRKMjUa773nXmOmK5SU4bffv87zeZ7z3H32XmuvtfeVczhlclyAj49P
kp+P77/ngeO9H8u/m2tl+BUNFO/S0RLcs9iTLasGgcBXEy7R8tjyMC4+X1YH
MkUbQrNuqkK/rikvQdYMeTohl71uaMPy2HhptOwxDBXvM1gUagLmxtA9N2R9
YCmp0XX2vAOSpA9dvSYbDAEKI9pHvXFZOF8pXPYOuj0ME/9uCIKVdXLO8I9U
2Ch3muZejobozGWbfctzwN4natXPvY/VxkO9b64X4HH674GlNsnwrd+OuLFi
mJ64ouV+/DFkCiRY/PxlePRGa8+G8Sco5O6pbdCswF1zvRN5iZmocE1YYj1Q
iXRjpSbl6Rz8VNw5P8v9HW55L4eLXy4EJjIWLFldA5vGJ/3bJ/KwdeZIoWl2
LRzEGMtZugWoDVarf7u2AfpdS4xHz77Aipv9E0+HGrGvZX/6nvRC9I+zWV83
fYIQ639T+zOKkNrHqBGP/AyN4oW2azNeQnxZOq+z8wsOOIcFbo4oRoTLC6ck
hW/YnRC2861uCX7tyuh0DmzCkVyn+k/jJaBS+x/4BjfBRq9eaOufEgQJq5he
DW9CttHChtjJEog3Zr5OuU3e77aQseRnQNE+O5GV3gSRnYKBlSIMWATkWpk1
NqEryir78DIGXhQX1UKuGZ53i9J8djPgo1pZsPBNM9JeHqyyD2HASW1l0q3q
ZsTdKv0RHcqA1Vb/62L1zYgwNfeqDGdAY4eag0BLM5bOWeAqH8mAgHbq/KGh
ZoxkKmjXxDAQZRzsWrW4BYwUs5vv0xl45qmz0se7BXde+ZZfrWYg3jtJlOfb
AtWP5dbBNQzc8Jn85RLQgi+nn+teqmXA81xBnX1YC67pRy853cjAtksKF/QT
W3Cm6l2uTjMDlVGC3+RrWpD/YN3t2B4G2NmV1z/KtWJuRZn2pCCFxkHBAVWl
VuRKf55HzaZQrrLX4I5KKyYKPDTPC1NIelYlZrmtFZGe+pU/RCk4pL290W7Q
iienRi7WSVLoTqiO6vFthVD800fbVlIYjKi9PVXXCkV56uW37RRYH8TG7D63
4t1V6o+DBoV6EUOz8uZWyJ9PTRrUpJAdVid1paMVOqpbF02Dgmdw/V3hiVb4
6mUXiupSGLvQGPuPQhvaokzyW4wp/HX98kD5fBsOBLFKdzlTMNtqe2tDYBt8
G8fcglwoZPJ3X90S3IYPW7NeVLhSsIsb99oZ2YaF68Sk4E7hddXi/caP2nDW
qFhQyZtC6MqjY+dq2tC992JR/gUK/3wbNKhe3A6rzsxT5TcouCWf1a5b0Y7q
9EyLr5EUKt1mtn1SaEfT/vCm3igKZ2dJrGrf0I7vvew00WgKLZvUJnja7egI
2+y/K4bCw+vnH8ucbEfkwW8CJxIprNs9Z9qJ0Y51V1Svp2dReP4sQ6m5vB3u
T2avT8mmoCZlZKz3th1ON2TWJuRQ0OyJS133sR2W9b/UI3IpHLyxVn+4ux2t
wsw7Fi8oeLcYxvlJMuF/oM4znqLwS3vkdZ8MExnP8C6glMKFzJj+f5cxoaAn
4mf7ikJwIFtLaw0TSwqzh2TKKcSs9uoS0GJCvvWIvN8bCsVed1QjXJiIfViZ
OfCeglbb/2ymPJjoMT1bnF9DoUKHGeLhw8Spkjg5vw8UamQUW4wCmSg1Uhec
qaXQXloUKH2XifaWMcnxBgp8Im21ieVMRK0stYz7SmHfI7kTuTIsWIwZhT5k
k/nGiSu+XcpCaukv6zU0hZWRfztb5Vjg/JV7lUP8x++rg6AKCzlaYholHApZ
xuF2lrtYmLN535xXHSRe/D/Mp5xZeNIwt+9SDxlvnCkl6cHC293OX0eIBwZq
PiueZkGMYWXi2EvhbXOqsVEAC1TMCN+u7xT8nlsZptxiwVFj2ebePgps+3Id
vRIWNvfTzNFBCrXmWbPsylgQbcrhHhyi8FL/foV3JQsC85b4PSGO3uazK76e
BeXY6sOWPAp756/R5HFZOO4Q9SNlmIKqkPTfWb0sGA4eujJKvPwvf8miQRZk
Z+em7f5B4Xd32zbtcRYsmywetRBnvIracleUDa3NWu6/RyjEFfiPPF3AhrLP
ggjNURK/Zy65r6TZGB67qhlIbBezZ2PvCjY0bRxM+MYoSHiMr9VQY0P3J9P4
B/G0Y+f3Q9vZ+OLKnFL+SaHP+uOTY1psnNsRrHGUuHJvhmLkfjY2Hnh/8gOx
7zJ7eY4NG7nWO5Iif5H9brTJov0oGwnHMyZKiA9e4b/W5MSG/Gmh0S5ixe8p
P+o82XB/cKZq2zgFyaVnVr8/w8ZC9YhEO+Kpg7rWlefZWFkfIRNM/KWgt7wk
mI1nnQbfa4jLe4rHXoSzMd0gpj1AnLH4+prcSDYOxxTtEPtNIejS+pvpcWzY
nasX3U/skT/95lECGzXDOzodia27G8YTHrHxdira5iLxJgNv+zuZbFzynb0+
h3hp4J7bUbls/NgoH1JFPCdP6t21QvL7UmsDWolHOrv+hDDIeu9dLj1EzCKN
+nI5G4KcaQO+CZLfB8Ic/KvY6JlVqyBB/CLAOsa3hg2R6pD4lcRJz9fWeDew
4bRXpWADcUTH3yn3L2wwb5Wd1CA+t7Bus0sLGyte440usaNe4vFjLDaqOzNL
jIgP+p+6Z9/BhoHYHBNr4h05u+oO97IxqXMw/AixIleC32KQjHc/6IgTsYR0
h5rxCBshUqlNrsST+wqcDcbJfBh5Q+7EvedD4vdNsmF8O/vxKeI+i3XOyXw0
ZB/GTv7nI8P/CKyYRUON7T7qQfwtfPJBvBANQZPN19yI9Vd1qS8WpvFllFPu
TFxO1TXEitAoexUY60i81aLQRXoejcpnIuK2xJm8xFm359MozwpSMCeWDw9N
WCBJw5vR06pPHCfvuS1SioZ71fbNu4nFKauPojI0lBnnlNSJg813nwyXpTET
nlqqRDwxpCw0ZymNZctKRxcRe4RJPgxeTiPWtuK9MHGn3N//CcjRsFYtwDiJ
V6NZrdvUaho1yraijcS6QwWzL6yhkTVf6mgJMRWakPR7LY1dC4osUojTSzy+
jG6k4dpYLuf1X36YWXp4qdJQ0FYaNyeOHoQwT43MX9XPbTvxxZUSmn3baQwy
WVsnSX6OFU98PaFJw2X7yPU2YldT7qkuLRr8FsOXionNr+an0HtoaE9lHPAi
VhkwF20ypGFWtNi2keyXl2YGZdeMaCz8N3MgiXhPqfZpLVMaUbeZzFPEhyM3
tj22ohFhvaxLmDh809wMH0caYQVzQpXI/pS6N2O3xolG8I3rq/rJ/k7i/yXZ
7kxjS3uGXBbxy0/c83s8aExLlmurEHf7UAck/WiYflBeJEPqhTczb7rqAhlP
P/h1HakvU3uf5vldpHH0uG9tELGUTMxizhUad7F6Ti+pR3uKT/XlRNIIZEbN
f0jqV9KU/DXDVBp6obJuZaQ+qjjJavGn0zg0t2iDDXFR/fyRgqc0vNQ7TMYH
SL9/OGm9NIfGiyedZWuIp3Z/U+4vJuu73FflHKm3h6+Gvw9roJGQr3zofTf5
3qHAAM1PNO4vNF9rQOxleXbz8Bca1RA4UddFIWyN4z3LVjKeDkvwQyeFwhpN
F8UuGqXzrnzK45J+MH9YuPIPDdWsxwVqLAqfSq+F6k2RfPBP+5DOJPF1Wz27
YYaG0PAjP1liyZrDs1oFOVDfdv3+7zbikKopnjgH/pMzvMctFBb8vT+yZBUH
3OwnvdQXkm9P1b2SV3OQ/GHMQ5b4plUjT3ENB53/ro/w+UxhfqHQ4Ob1HNhR
9tcVP5F89/Ls2beNg4a9mfF+pB+K9extP63PwZZjo3PbSX+tvUsfnjDkwOz3
1Wk54gidCy0XjTjIufUz3Kmaguij59/CzTmQXmt/vP8tBRHbJR8f2nOQX6e9
oZP0a+HPw1UfTnPw+rjGnGjS7wXLHuSsiucgZv/vWhY5P7Rq6uxnJ3Lwc6uz
cic5X+SUDND3kzlY8N2D2UPOH4cLd0pKpnOwY/t2g74MCrlZbB++PA46dt3Q
aUsn/SV+1Q7mOw6Eg64muj8k+eWXURkzyoFfd9F3U3I+ihw3tTX5xYHo0y3q
ChGkvp2Z/DlvggO1g0LNI9f+m7+hUsg0B0/MklddD6Nw3JkXfnouF45bAj6l
XyHrbbnlkNEKLjwWvY645UfhpBqjZa4+F4fYpQ6Kx0j9afya+dKQC93jZUuq
HShIuw8HnjDiQtNPWcP5KIWyx6sVq8y5OOm4PivFjsRf5qb3pSNcvNfLtp1j
TYExcUx0/AwXC7f8O3TZkKxnmfjOziQutvfxUXvUSb/Rc0h69YuLjdGz3CZm
GFBTkdPCBBcVjXkqptMMvBKn28r+cpFJnQrLmGSg8bOtzGu+DvjsFJawmmBg
zNY68o1IB0wbXJ2SRhjQ8jbyr17egRn3tBVDneT9+1pWn3Q7MP56dnEAOc//
7F+8oDumAyJaJb1u1xm4k9DzLv5eBxSl4iS1yX1hy6GCQNP4DlgKf4uVJvcJ
rzxDXnlyBz7udSgoCmJg8Nyl+geZHUjcp6/e68dAl1B3hHFFB7LFGkz+nGDg
64rcuaWDHajpv8DmapP7i+m+WXd0OtG60MOgcrQEgftDHhZ0d6I8TVHP5WAJ
5JgHmlNCujCVfOnb2N1i7NgWfEdWtRtYHiCaU/ASNssbJvg/d+OBwti5bzVF
+Flzg+t5uQcj+9bddcsvRJygnY+JQi8K/zIjJWJe4B/Wi0cClb3I8QlP0z1R
AKGAy2bbvb4jsOZrVr9gPvie2uVwJfsgcmdqASM5F973c25IVvShv2pkyZTQ
c5SuCxOW8OjHs+GUvG+ZmZA44K0SJD6Alcze6CCLp2hf3/jZjxqA1bnpyNVB
j/E8EQqddoM42Zdn/o6RDNucrRLLpgehVPchWPv5A1wKSb7c/mwIF1Q9Deo+
RSOrd9kDcQMexBdUVzbpByMsIeQI6yAPMrdKlWInr+CYyeDqbGMeBCequLpZ
V7CUop4bWvLQ4Rhv5yt+Bdei/n0b4cCDtMIWm6WNl+G89d4PET8eXA4H1PQn
XYR8sNT+OWk8pO7XL6i65oPp//nPa3rCw+3HbYutek+jZbDjU1oGD3l/vs/9
sPc0blrm2+rm8rDU+IK/CZ83ZtaanA6heHD76b+yoccD7R+jEmd95kH9wccc
k2fH8PLquOPnrzyMtJ+xCz3qiDsa9sopzTyUNOT/9F/kAP3HGwq0WTw42Kel
m4vao+Rc3fvL33nQqKgQly8xQ8x69SjjAR6sNTeduKdgAi9uvJkcj4fQV/KN
CssOYY2BG7t8jIdbZ88lSFboQJD/S+rNcR6SdZSfLlUE6Bcarkf+8ECpm55Y
/FYNlGvKxk1TPPS5tr3WrFZA7ArRnzMzPAxalwz+93/I/wERs98r
        "]]},
      Annotation[#, 
       "Charting`Private`Tag$2092#2"]& ], {}}, {{}, {}, {}, {}}}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
         Part[#, 1]], 
        (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
         Part[#, 1]], 
        (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 1}, {0., 0.11368324458330914`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.824813132338661*^9, 3.824813399908742*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.8244658228259525`*^9, 3.82446582282695*^9}},
 ExpressionUUID -> "f21a279d-5eb1-4292-bcb8-e69a29b0f7ad"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.82446582384227*^9, 3.8244658238442645`*^9}},
 ExpressionUUID -> "065e843c-91d2-416a-a5ce-1947b9520555"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.8244658253372765`*^9, 3.8244658253382683`*^9}},
 ExpressionUUID -> "4dfd2d26-d6cd-4362-aaf1-fd3257a82e72"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.8244630845037374`*^9, 3.8244630865670843`*^9}},
 ExpressionUUID -> "44368122-897c-47c4-ab2b-d2e2f3b99f2f"]
},
WindowSize->{1540, 825},
WindowMargins->{{Automatic, -1}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PrintingOptions->{"PaperOrientation"->"Portrait",
"PaperSize"->{Automatic, Automatic},
"PostScriptOutputFile"->"/home/pedro/print.pdf"},
Magnification:>3.1,
FrontEndVersion->"11.0 for Microsoft Windows (64-bit) (September 21, 2016)",
StyleDefinitions->"Default.nb",
ExpressionUUID -> "f90bf763-7443-4fa6-b861-1425cd32762f"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 11498, 299, 3027, "Input"],
Cell[12081, 323, 732, 23, 331, "Output"],
Cell[12816, 348, 1152, 29, 331, "Output"],
Cell[13971, 379, 11998, 217, 734, "Output"]
}, Open  ]],
Cell[25984, 599, 152, 2, 132, "Input"],
Cell[26139, 603, 152, 2, 132, "Input"],
Cell[26294, 607, 155, 2, 132, "Input"],
Cell[26452, 611, 155, 2, 132, "Input"]
}
]
*)