(* Content-type: application/mathematica *)

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

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     25207,        706]
NotebookOptionsPosition[     24288,        673]
NotebookOutlinePosition[     24778,        691]
CellTagsIndexPosition[     24735,        688]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
1\.aa Practica: Sistemas Lineales (M\[EAcute]todos directos)\
\>", "Title"],

Cell[CellGroupData[{

Cell[TextData[{
 StyleBox["Ejercicio 1\.ba",
  FontVariations->{"Underline"->True}],
 ".\nSea el siguiente sistema:\n\n",
 Cell[BoxData[
  InterpretationBox[GridBox[{
     {
      RowBox[{
       RowBox[{
        RowBox[{"2.1754`", " ", "x1"}], "+", 
        RowBox[{"4.0231`", " ", "x2"}], "-", 
        RowBox[{"2.1732`", " ", "x3"}], "+", 
        RowBox[{"5.1967`", " ", "x4"}]}], "==", "17.108`"}]},
     {
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"-", "4.0231`"}], " ", "x1"}], "+", 
        RowBox[{"6", " ", "x2"}], "+", 
        RowBox[{"1.1973`", " ", "x4"}]}], "==", 
       RowBox[{"-", "6.1593`"}]}]},
     {
      RowBox[{
       RowBox[{
        RowBox[{"-", "x1"}], "-", 
        RowBox[{"5.2107`", " ", "x2"}], "+", 
        RowBox[{"1.1111`", " ", "x3"}]}], "==", "3.0004`"}]},
     {
      RowBox[{
       RowBox[{
        RowBox[{"6.0235`", " ", "x1"}], "+", 
        RowBox[{"7", " ", "x2"}], "-", 
        RowBox[{"4.1561`", " ", "x4"}]}], "==", "0"}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[2.0999999999999996`]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}],
   TableForm[{
    2.1754 x1 + 4.0231 x2 - 2.1732 x3 + 5.1967 x4 == 
     17.108, (-4.0231) x1 + 6 x2 + 1.1973 x4 == -6.1593, -x1 - 5.2107 x2 + 
      1.1111 x3 == 3.0004, 6.0235 x1 + 7 x2 - 4.1561 x4 == 0}]]]],
 "\n\nResolverlo utilizando los siguientes procedimientos de ",
 StyleBox["Mathematica",
  FontSlant->"Italic"],
 " :\n\t- Solve\n\t- LinearSolve\n\t- Mediante las funciones preparadas de \
eliminaci\[OAcute]n gaussiana (GaussElim []) y de pivotaje parcial escalado \
(PivParEsc[])\n\t- LUDecomposition"
}], "Section",
 CellChangeTimes->{{3.5075623928697805`*^9, 3.5075624159635305`*^9}}],

Cell["1\.ba.- Comando Solve", "Subsection",
 CellChangeTimes->{{3.506342019546875*^9, 3.50634202253125*^9}}],

Cell[TextData[{
 "2\.ba.- El comando ",
 StyleBox["LinearSolve", "MR"],
 " sirve para sistemas en que la matriz de coeficientes es tanto \
num\[EAcute]rica como simb\[OAcute]lica. La matriz puede ser cuadrada o \
rectangular. en el caso de sistemas indeterminados dara como resultado una de \
las posibles soluciones mientras que solve da la soluci\[OAcute]n general.\n\
La salida que se obtiene tambi\[EAcute]n es distinta.\nCuando es necesario \
resolver el mismo sistema con distintos vectores de t\[EAcute]rminos \
independientes es mucho m\[AAcute]s eficiente utilizar el comandos \
LUDecomposition.\n"
}], "Subsection",
 TextAlignment->Left,
 TextJustification->1],

Cell[CellGroupData[{

Cell["\<\
3\.ba.- Se utilizan dos funciones creadas por nosotros para ilustrar el \
procedimiento de resoluci\[OAcute]n de istemas de ecuaciones lineales \
mediante la eliminaci\[OAcute]n gaussiana y mediante el pivotaje parcial \
escalado.
Se ejecutan las funciones para que mathematica las identifique y sepa como \
operar y a continuaci\[OAcute]n se aplican a las matrices de coeficientes y \
termino independiente que se desee.\
\>", "Subsection"],

Cell[BoxData[
 RowBox[{
  RowBox[{"GaussElim", " ", "[", 
   RowBox[{"AA_", ",", "bb_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "x", ",", "c", ",", "m", ",", "n", ",", "p", ",", "i", ",", "j", ",", 
      "k", ",", "A", ",", "NoUniqueSoln"}], "}"}], ",", 
    RowBox[{
     RowBox[{"A", "=", 
      RowBox[{"Transpose", "[", 
       RowBox[{"Append", "[", 
        RowBox[{
         RowBox[{"Transpose", "[", "AA", "]"}], ",", "b"}], "]"}], "]"}]}], 
     ";", 
     RowBox[{"n", "=", 
      RowBox[{"Length", "[", 
       RowBox[{"A", "[", 
        RowBox[{"[", 
         RowBox[{"All", ",", "1"}], "]"}], "]"}], "]"}]}], ";", 
     RowBox[{"m", "=", "A"}], ";", 
     RowBox[{"x", "=", 
      RowBox[{"A", "[", 
       RowBox[{"[", 
        RowBox[{"All", ",", "1"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"p", "=", "1"}], "  ", ";", 
     RowBox[{"NoUniqueSoln", "=", "0"}], ";", " ", 
     RowBox[{"Print", "[", 
      RowBox[{"\"\<Sistema Inicial\>\"", ",", 
       RowBox[{"MatrixForm", "[", "A", "]"}]}], "]"}], ";", 
     "\[IndentingNewLine]", "  ", 
     RowBox[{"For", "\[IndentingNewLine]", "[", 
      RowBox[{
       RowBox[{"i", "=", "1"}], ",", 
       RowBox[{"i", "<", "n"}], ",", 
       RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"Print", "[", 
         RowBox[{"\"\<Etapa\>\"", ",", "i"}], "]"}], ";", 
        "\[IndentingNewLine]", "\t", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"A", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "i"}], "]"}], "]"}], "\[NotEqual]", "0"}], ",", 
          RowBox[{"p", "=", "i"}], ",", "\[IndentingNewLine]", "\t", 
          RowBox[{
           RowBox[{"p", "=", "i"}], ";", 
           RowBox[{"While", "[", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{"p", "\[LessEqual]", "n"}], ")"}], "&&", 
              RowBox[{"(", 
               RowBox[{"p", "\[GreaterEqual]", "i"}], ")"}], "&&", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"Abs", "[", 
                 RowBox[{"A", "[", 
                  RowBox[{"[", 
                   RowBox[{"p", ",", "i"}], "]"}], "]"}], "]"}], "<", 
                SuperscriptBox["10", 
                 RowBox[{"-", "15"}]]}], ")"}]}], ",", 
             RowBox[{"p", "=", 
              RowBox[{"p", "+", "1"}]}]}], "]"}]}]}], "]"}], ";", 
        "\[IndentingNewLine]", "\t", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{"p", ">", "n"}], ",", 
          RowBox[{
           RowBox[{"NoUniqueSoln", "=", "1"}], ";", "Abort"}], ",", 
          RowBox[{
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{"p", "\[NotEqual]", "i"}], ",", 
             RowBox[{
              RowBox[{"For", "[", 
               RowBox[{
                RowBox[{"j", "=", "1"}], ",", 
                RowBox[{"j", "\[LessEqual]", 
                 RowBox[{"n", "+", "1"}]}], ",", 
                RowBox[{"j", "++"}], ",", 
                RowBox[{
                 RowBox[{"c", "=", 
                  RowBox[{"A", "[", 
                   RowBox[{"[", 
                    RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ";", 
                 RowBox[{
                  RowBox[{"A", "[", 
                   RowBox[{"[", 
                    RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", 
                  RowBox[{"A", "[", 
                   RowBox[{"[", 
                    RowBox[{"p", ",", "j"}], "]"}], "]"}]}], ";", 
                 RowBox[{
                  RowBox[{"A", "[", 
                   RowBox[{"[", 
                    RowBox[{"p", ",", "j"}], "]"}], "]"}], "=", "c"}]}]}], 
               "]"}], ";", 
              RowBox[{"Print", "[", 
               RowBox[{
               "\"\<F\>\"", ",", "p", ",", "\"\<<->F\>\"", ",", "i", ",", 
                RowBox[{"MatrixForm", "[", "A", "]"}]}], "]"}]}], ",", 
             "Continue"}], "]"}], 
           RowBox[{"For", "[", 
            RowBox[{
             RowBox[{"j", "=", 
              RowBox[{"i", "+", "1"}]}], ",", 
             RowBox[{"j", "\[LessEqual]", "n"}], ",", 
             RowBox[{"j", "++"}], ",", 
             RowBox[{
              RowBox[{
               RowBox[{"m", "[", 
                RowBox[{"[", 
                 RowBox[{"j", ",", "i"}], "]"}], "]"}], "=", 
               RowBox[{
                RowBox[{"A", "[", 
                 RowBox[{"[", 
                  RowBox[{"j", ",", "i"}], "]"}], "]"}], "/", 
                RowBox[{"A", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "i"}], "]"}], "]"}]}]}], ";", " ", 
              "\[IndentingNewLine]", 
              RowBox[{"A", "=", 
               RowBox[{"ReplacePart", "[", 
                RowBox[{"A", ",", 
                 RowBox[{
                  RowBox[{"A", "[", 
                   RowBox[{"[", "j", "]"}], "]"}], "-", 
                  RowBox[{
                   RowBox[{"m", "[", 
                    RowBox[{"[", 
                    RowBox[{"j", ",", "i"}], "]"}], "]"}], "*", 
                   RowBox[{"A", "[", 
                    RowBox[{"[", "i", "]"}], "]"}]}]}], " ", ",", "j"}], 
                "]"}]}], ";", "\[IndentingNewLine]", 
              RowBox[{"Print", "[", 
               RowBox[{"\"\<F\>\"", ",", "j", ",", "\"\<-(\>\"", ",", 
                RowBox[{"m", "[", 
                 RowBox[{"[", 
                  RowBox[{"j", ",", "i"}], "]"}], "]"}], ",", "\"\<)F\>\"", 
                ",", "i", ",", 
                RowBox[{"MatrixForm", "[", 
                 RowBox[{"Chop", "[", "A", "]"}], "]"}]}], "]"}]}]}], 
            "\[IndentingNewLine]", "]"}]}]}], "  ", "]"}]}]}], " ", "]"}], 
     "  ", ";", "  ", "\[IndentingNewLine]", 
     RowBox[{"If", "\[IndentingNewLine]", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"A", "[", 
           RowBox[{"[", 
            RowBox[{"n", ",", "n"}], "]"}], "]"}], "==", "0"}], ")"}], " ", "||", 
        RowBox[{"(", " ", 
         RowBox[{"NoUniqueSoln", "==", "1"}], ")"}]}], ",", "Abort", ",", 
       "\[IndentingNewLine]", "\t\t", 
       RowBox[{
        RowBox[{
         RowBox[{"x", "[", 
          RowBox[{"[", "n", "]"}], "]"}], "=", 
         RowBox[{
          RowBox[{"A", "[", 
           RowBox[{"[", 
            RowBox[{"n", ",", 
             RowBox[{"n", "+", "1"}]}], "]"}], "]"}], "/", 
          RowBox[{"A", "[", 
           RowBox[{"[", 
            RowBox[{"n", ",", "n"}], "]"}], "]"}]}]}], " ", ";", 
        "\[IndentingNewLine]", 
        RowBox[{"For", "[", 
         RowBox[{
          RowBox[{"i", "=", 
           RowBox[{"n", "-", "1"}]}], ",", 
          RowBox[{"i", "\[GreaterEqual]", "1"}], ",", 
          RowBox[{"i", "--"}], ",", 
          RowBox[{
           RowBox[{"x", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"A", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", 
                 RowBox[{"n", "+", "1"}]}], "]"}], "]"}], "-", 
              RowBox[{"Sum", "[", 
               RowBox[{
                RowBox[{
                 RowBox[{"A", "[", 
                  RowBox[{"[", 
                   RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", 
                 RowBox[{"x", "[", 
                  RowBox[{"[", "j", "]"}], "]"}]}], ",", 
                RowBox[{"{", 
                 RowBox[{"j", ",", 
                  RowBox[{"i", "+", "1"}], ",", "n"}], "}"}]}], "]"}]}], 
             ")"}], "/", 
            RowBox[{"A", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "i"}], "]"}], "]"}]}]}]}], "  ", "]"}]}]}], 
      "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"Print", "[", 
      RowBox[{"\"\<Sistema Triangularizado\>\"", ",", 
       RowBox[{"MatrixForm", "[", 
        RowBox[{"Chop", "[", "A", "]"}], "]"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"NoUniqueSoln", "\[Equal]", "1"}], " ", ",", 
       RowBox[{
       "Print", "[", "\"\<No existe soluci\[OAcute]n \[UAcute]nica\>\"", 
        "]"}], ",", 
       RowBox[{"Print", "[", 
        RowBox[{"\"\<Soluci\[OAcute]n=\>\"", ",", "x"}], "]"}]}], "]"}]}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.431773896140625*^9, 3.431773912203125*^9}, {
  3.431774004421875*^9, 3.43177401909375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"PivParEsc", "[", 
   RowBox[{"AA_", ",", "B_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "s", ",", "per", ",", "i", ",", "j", ",", "k", ",", "calpiv", ",", 
      "pospiv", ",", "el", ",", "m", ",", "a", ",", "b", ",", "sol"}], "}"}], 
    ",", 
    RowBox[{
     RowBox[{"a", "=", "AA"}], ";", 
     RowBox[{"b", "=", "B"}], ";", "\[IndentingNewLine]", 
     RowBox[{"s", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Abs", "[", 
         RowBox[{"Max", "[", 
          RowBox[{"Abs", "[", 
           RowBox[{"a", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", 
          RowBox[{"Length", "[", "a", "]"}]}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"per", "=", 
      RowBox[{"Table", "[", 
       RowBox[{"i", ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", 
          RowBox[{"Length", "[", "a", "]"}]}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "1"}], ",", 
       RowBox[{"i", "<", 
        RowBox[{"Length", "[", "a", "]"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"calpiv", "=", 
         RowBox[{"Table", "[", 
          RowBox[{
           RowBox[{"Abs", "[", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{
                RowBox[{"per", "[", 
                 RowBox[{"[", "j", "]"}], "]"}], ",", "i"}], "]"}], "]"}], 
             "/", 
             RowBox[{"s", "[", 
              RowBox[{"[", 
               RowBox[{"per", "[", 
                RowBox[{"[", "j", "]"}], "]"}], "]"}], "]"}]}], "]"}], ",", 
           RowBox[{"{", 
            RowBox[{"j", ",", "i", ",", 
             RowBox[{"Length", "[", "per", "]"}]}], "}"}]}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"pospiv", "=", 
         RowBox[{"Flatten", "[", 
          RowBox[{"Position", "[", 
           RowBox[{"calpiv", ",", 
            RowBox[{"Max", "[", "calpiv", "]"}]}], "]"}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"pospiv", "=", 
         RowBox[{
          RowBox[{"pospiv", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "+", 
          RowBox[{"(", 
           RowBox[{"i", "-", "1"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"el", "=", 
         RowBox[{"per", "[", 
          RowBox[{"[", "pospiv", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"per", "=", 
         RowBox[{"Delete", "[", 
          RowBox[{"per", ",", "pospiv"}], "]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"per", "=", 
         RowBox[{"Insert", "[", 
          RowBox[{"per", ",", "el", ",", "i"}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"Print", "[", 
         RowBox[{"\"\<Etapa \>\"", ",", "i", ",", "\"\<   Pivote \>\"", ",", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"el", ",", "i"}], "]"}], "]"}], ",", 
          "\"\< situado en la posici\[OAcute]n (\>\"", ",", "el", ",", 
          "\"\<,\>\"", ",", "i", ",", "\"\<)\>\""}], "]"}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"Print", "[", 
         RowBox[{"\"\<Vector de permutaci\[OAcute]n:\>\"", ",", "per"}], 
         "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
        RowBox[{"For", "[", 
         RowBox[{
          RowBox[{"j", "=", 
           RowBox[{"i", "+", "1"}]}], ",", " ", 
          RowBox[{"j", "\[LessEqual]", 
           RowBox[{"Length", "[", "a", "]"}]}], ",", 
          RowBox[{
           RowBox[{"m", "=", 
            FractionBox[
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{
                RowBox[{"per", "[", 
                 RowBox[{"[", "j", "]"}], "]"}], ",", "i"}], "]"}], "]"}], 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"el", ",", "i"}], "]"}], "]"}]]}], ";", 
           "\[IndentingNewLine]", "\[IndentingNewLine]", 
           RowBox[{"For", "[", 
            RowBox[{
             RowBox[{"k", "=", 
              RowBox[{"i", "+", "1"}]}], ",", 
             RowBox[{"k", "\[LessEqual]", 
              RowBox[{"Length", "[", "a", "]"}]}], ",", "\[IndentingNewLine]", 
             RowBox[{
              RowBox[{
               RowBox[{"a", "[", 
                RowBox[{"[", 
                 RowBox[{
                  RowBox[{"per", "[", 
                   RowBox[{"[", "j", "]"}], "]"}], ",", "k"}], "]"}], "]"}], 
               "=", 
               RowBox[{
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{
                   RowBox[{"per", "[", 
                    RowBox[{"[", "j", "]"}], "]"}], ",", "k"}], "]"}], "]"}], 
                "-", 
                RowBox[{"m", "*", 
                 RowBox[{"a", "[", 
                  RowBox[{"[", 
                   RowBox[{"el", ",", "k"}], "]"}], "]"}]}]}]}], ";", 
              RowBox[{"k", "++"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
           "\[IndentingNewLine]", "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{
               RowBox[{"per", "[", 
                RowBox[{"[", "j", "]"}], "]"}], ",", "i"}], "]"}], "]"}], "=",
             "0"}], ";", 
           RowBox[{
            RowBox[{"b", "[", 
             RowBox[{"[", 
              RowBox[{"per", "[", 
               RowBox[{"[", "j", "]"}], "]"}], "]"}], "]"}], "=", 
            RowBox[{
             RowBox[{"b", "[", 
              RowBox[{"[", 
               RowBox[{"per", "[", 
                RowBox[{"[", "j", "]"}], "]"}], "]"}], "]"}], "-", 
             RowBox[{"m", "*", 
              RowBox[{"b", "[", 
               RowBox[{"[", "el", "]"}], "]"}]}]}]}], ";", 
           RowBox[{"j", "++"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
        RowBox[{"Print", "[", 
         RowBox[{"MatrixForm", "[", "a", "]"}], "]"}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"Print", "[", 
         RowBox[{"\"\<T\[EAcute]rmino independiente=\>\"", ",", "b"}], "]"}], 
        ";", "\[IndentingNewLine]", 
        RowBox[{"i", "++"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
     "\[IndentingNewLine]", 
     RowBox[{
     "Print", "[", "\"\<Por sustituci\[OAcute]n regresiva se obtiene: \>\"", 
      "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
     RowBox[{"sol", "=", 
      RowBox[{"Table", "[", 
       RowBox[{"0", ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", 
          RowBox[{"Length", "[", "a", "]"}]}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"sol", "[", 
       RowBox[{"[", 
        RowBox[{"Length", "[", "a", "]"}], "]"}], "]"}], "=", 
      FractionBox[
       RowBox[{"b", "[", 
        RowBox[{"[", 
         RowBox[{"per", "[", 
          RowBox[{"[", 
           RowBox[{"Length", "[", "a", "]"}], "]"}], "]"}], "]"}], "]"}], 
       RowBox[{"a", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"per", "[", 
           RowBox[{"[", 
            RowBox[{"Length", "[", "a", "]"}], "]"}], "]"}], ",", 
          RowBox[{"Length", "[", "a", "]"}]}], "]"}], "]"}]]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", 
        RowBox[{
         RowBox[{"Length", "[", "a", "]"}], "-", "1"}]}], ",", 
       RowBox[{"i", "\[GreaterEqual]", "1"}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"sol", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "=", 
         FractionBox[
          RowBox[{
           RowBox[{"b", "[", 
            RowBox[{"[", 
             RowBox[{"per", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}], "-", 
           RowBox[{
            UnderoverscriptBox["\[Sum]", 
             RowBox[{"j", "=", 
              RowBox[{"i", "+", "1"}]}], 
             RowBox[{"Length", "[", "a", "]"}]], 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{
                RowBox[{"per", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], ",", "j"}], "]"}], "]"}], 
             "*", 
             RowBox[{"sol", "[", 
              RowBox[{"[", "j", "]"}], "]"}]}]}]}], 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{
             RowBox[{"per", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", "i"}], "]"}], "]"}]]}], 
        ";", "\[IndentingNewLine]", 
        RowBox[{"i", "--"}], ";"}]}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"Print", "[", 
      RowBox[{"\"\<Solucion= \>\"", ",", "sol"}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Return", "[", 
      RowBox[{"{", 
       RowBox[{"a", ",", "b", ",", "sol"}], "}"}], "]"}]}]}], 
   "]"}]}]], "Input"]
}, Open  ]],

Cell[TextData[{
 "4\.ba.-  El comando LUDecomposition realiza una descomposicion de A en \
producto de una triangular inferior por una triangular superior con \
intercambio de filas.",
 StyleBox[" LUDecomposition,", "MR"],
 " devuelve una lista tres elementos. El primero es una combinacion de las \
matrices triangular inferior y superior, el segundo es un vector que \
especifica las filas que se han usado en el pivoteo, y el tercero es una \
estimaci\[OAcute]n del n\.ba de condici\[OAcute]n de la matriz. El comando \
LUMatrices[ ] transforma el resultado de la matriz conjunta que devuelve el \
comando LUDecomposition, en dos matrices L y U, de modo que ya se puede \
realizar el proceso de resolucion del sistema en dos etapas. Se debe tener en \
cuenta que si el comando LUDecomposition ha realizado intercambio de filas, \
tambien se debera aplicar dicho intercambio a el vector b , termino \
independiente del sistema que se pretende resolver."
}], "Subsection",
 TextAlignment->Left,
 TextJustification->1]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 StyleBox["Ejercicio 2\.ba",
  FontVariations->{"Underline"->True}],
 "\nResolver el siguiente sistema usando la factorizaci\[OAcute]n de \
Choleski:\n6",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "1"], TraditionalForm]]],
 "+2",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "2"], TraditionalForm]]],
 "+",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "3"], TraditionalForm]]],
 "-",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "4"], TraditionalForm]]],
 "=0\n2",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "1"], TraditionalForm]]],
 "+4",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "2"], TraditionalForm]]],
 "+",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "3"], TraditionalForm]]],
 "     =7\n",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "1"], TraditionalForm]]],
 "+",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "2"], TraditionalForm]]],
 "+4",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "3"], TraditionalForm]]],
 "-",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "4"], TraditionalForm]]],
 "  =-1\n-",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "1"], TraditionalForm]]],
 "-",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "3"], TraditionalForm]]],
 "+3",
 Cell[BoxData[
  FormBox[
   SubscriptBox["x", "4"], TraditionalForm]]],
 "      =-2"
}], "Section"],

Cell["\<\
La descomposici\[OAcute]n de  Cholesky se aplica a  una matriz \
sim\[EAcute]trica definida positiva y da como resultado el producto de una \
matriz triangular superior por una triangular inferior que es la traspuesta \
de la anterior. \
\>", "Subsection",
 CellChangeTimes->{3.5075623317760305`*^9},
 TextAlignment->Left,
 TextJustification->1]
}, Open  ]]
}, Open  ]]
},
WindowToolbars->{"RulerBar", "EditBar"},
WindowSize->{1028, 746},
WindowMargins->{{3, Automatic}, {Automatic, 30}},
CellLabelAutoDelete->True,
FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (February 18, 2009)",
StyleDefinitions->FrontEnd`FileName[{"Report"}, "StandardReport.nb", 
  CharacterEncoding -> "WindowsANSI"]
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[567, 22, 85, 2, 81, "Title"],
Cell[CellGroupData[{
Cell[677, 28, 2019, 55, 401, "Section"],
Cell[2699, 85, 108, 1, 24, "Subsection"],
Cell[2810, 88, 670, 13, 127, "Subsection"],
Cell[CellGroupData[{
Cell[3505, 105, 451, 8, 75, "Subsection"],
Cell[3959, 115, 8546, 220, 456, "Input"],
Cell[12508, 337, 9003, 236, 792, "Input"]
}, Open  ]],
Cell[21526, 576, 1018, 17, 127, "Subsection"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22581, 598, 1321, 61, 203, "Section"],
Cell[23905, 661, 355, 8, 41, "Subsection"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)