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

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

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     23689,        881]
NotebookOptionsPosition[     19025,        712]
NotebookOutlinePosition[     19682,        739]
CellTagsIndexPosition[     19617,        734]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Instrucciones de control de programa", "Title"],

Cell[TextData[StyleBox["Para poder realizar la mayor\[IAcute]a de los \
programas que representan m\[EAcute]todos numer\[IAcute]cos de resoluci\
\[OAcute]n de problemas matem\[AAcute]ticos, se deben emplear sentencias de \
control.\nEste tipo de sentencias permiten realizar fundamentalmente \
operaciones de forma c\[IAcute]clica. La sustituci\[OAcute]n regresiva de un \
sistema triangular superior consta de n etapas. Todas ellas siguen una ley de \
recurrecia com\[UAcute]n, basta con indicarla una vez y pedir que se realice \
n veces. Los m\[EAcute]todos iterativos, necesitan de un n\[UAcute]mero \
indeterminado, a priori, de ciclos para conseguir una buena aproximaci\
\[OAcute]n a la soluci\[OAcute]n del problema.\nTodos estos problemas se \
resuelven con sentencias de control.",
 FontFamily->"Arial"]], "Text",
 TextAlignment->Left,
 TextJustification->1],

Cell[CellGroupData[{

Cell["Sentencia Do", "Subtitle"],

Cell[TextData[StyleBox["Su estructura es la siguiente: \
Do[expresi\[OAcute]n,{i,imin, imax,di}]\nRealiza la expresi\[OAcute]n que se \
indica dentro del comando desde imin hasta imax con un incremento di.\nSi no \
se indica el incremento di se considera di=1",
 FontFamily->"Arial"]], "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"sum", "=", "0"}], ";"}]], "Input",
 CellTags->"Do"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Do", "[", 
   RowBox[{
    RowBox[{"sum", "=", 
     RowBox[{"sum", "+", "i"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1000"}], "}"}]}], "]"}], ";"}]], "Input",
 CellTags->"Do"],

Cell[CellGroupData[{

Cell[BoxData["sum"], "Input"],

Cell[BoxData["500500"], "Output",
 CellChangeTimes->{3.5060574164375*^9, 3.506340797546875*^9, 
  3.50634084815625*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
   "En", "el", "ejemplo", "que", "se", "ha", "realizado", "se", "calcula", 
    "el", "sumatorio", "de", "i", "desde", "i"}], "=", 
   RowBox[{"1", "hasta", "1000"}]}], ",", 
  RowBox[{
   RowBox[{"como", "ya", "es", "sabido", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"i", "=", "1"}], "n"], "i"}]}], "=", 
   RowBox[{
    FractionBox[
     SuperscriptBox["n", "2"], "2"], "+", 
    FractionBox["n", "2"]}]}]}]], "Text",
 FontFamily->"Arial"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  FractionBox[
   SuperscriptBox["1000", "2"], "2"], "+", 
  FractionBox["1000", "2"]}]], "Input"],

Cell[BoxData["500500"], "Output",
 CellChangeTimes->{3.506057422453125*^9, 3.5063407975625*^9, 
  3.506340848171875*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  UnderoverscriptBox["\[Sum]", 
   RowBox[{"i", "=", "1"}], "1000"], "i"}]], "Input"],

Cell[BoxData["500500"], "Output",
 CellChangeTimes->{3.506057423234375*^9, 3.50634079759375*^9, 
  3.506340848203125*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Sentencia While", "Subtitle"],

Cell["\<\
Su estructura es la siguiente: While[test,ejecuci\[OAcute]n]
Realiza todo lo indicado en la parte de ejecuci\[OAcute]n, tantas sentencias \
como se desee, mientras se cumpla la condici\[OAcute]n del test.
Si no se indica el incremento di se considera di=1\
\>", "Text",
 FontFamily->"Arial"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"sum", "=", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"i", "=", "0"}], ";"}]}], "Input"],

Cell[BoxData[
 RowBox[{"While", "[", 
  RowBox[{
   RowBox[{"i", "\[LessEqual]", "1000"}], ",", 
   RowBox[{
    RowBox[{"sum", "=", 
     RowBox[{"sum", "+", "i"}]}], ";", 
    RowBox[{"i", "=", 
     RowBox[{"i", "+", "1"}]}], ";"}]}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["sum"], "Input"],

Cell[BoxData["500500"], "Output",
 CellChangeTimes->{3.506057428875*^9, 3.506340797625*^9, 3.50634084825*^9}]
}, Open  ]],

Cell["\<\
Se puede terminar el ciclo antes de que deje de cumplirse la \
condici\[OAcute]n del test con el comando Break[ ]\
\>", "Text",
 FontFamily->"Arial"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"sum", "=", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"i", "=", "0"}], ";"}]}], "Input"],

Cell[BoxData[
 RowBox[{"While", "[", 
  RowBox[{
   RowBox[{"i", "\[LessEqual]", "1000"}], ",", 
   RowBox[{
    RowBox[{"sum", "=", 
     RowBox[{"sum", "+", "i"}]}], ";", 
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{"sum", ">", "10000"}], ",", 
      RowBox[{"Break", "[", "]"}]}], "]"}], ";", 
    RowBox[{"i", "=", 
     RowBox[{"i", "+", "1"}]}], ";"}]}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["sum"], "Input"],

Cell[BoxData["10011"], "Output",
 CellChangeTimes->{3.5060574323125*^9, 3.506340797671875*^9, 
  3.50634084828125*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["i"], "Input"],

Cell[BoxData["141"], "Output",
 CellChangeTimes->{3.506057432953125*^9, 3.5063407976875*^9, 
  3.5063408483125*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
  "En", " ", "este", " ", "ejemplo", " ", "se", " ", "controla", " ", "que", 
   " ", "el", " ", "sumatorio", " ", "supere", " ", "10000"}], ",", " ", 
  RowBox[{
  "y", " ", "en", " ", "ese", " ", "momento", " ", "se", " ", "sale", " ", 
   "del", " ", "ciclo", " ", "while"}], ",", " ", 
  RowBox[{
   RowBox[{
   "almacenando", " ", "el", " ", "valor", " ", "de", " ", "i", " ", 
    "as\[IAcute]", " ", "como", " ", "la", " ", "suma", " ", 
    RowBox[{"realizada", ".", " ", "Por"}], " ", "tanto", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"i", "=", "1"}], "141"], "i"}]}], "=", "10011"}]}]], "Text",
 FontFamily->"Arial"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Sentencia For", "Subtitle"],

Cell[TextData[StyleBox["Su estructura es la siguiente: \
For[inicio,test,incremento,ejecuci\[OAcute]n]\nSe trata de un ciclo que tiene \
un inicio, se realizan todas las instrucciones del cuerpo ejecuci\[OAcute]n \
mientras se cumpla el test y se puede incrementar la variable inicio una \
cantidad en el apartado incremento. En ejecucion se pueden realizar todas las \
intrucciones que se desee mientras est\[EAcute]n separados por ; Se puede \
salir del ciclo con la sentencia Break[ ]\nEl incremento es opcional.",
 FontFamily->"Arial"]], "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"sum", "=", "0"}], ";"}]], "Input"],

Cell[BoxData[
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "\[LessEqual]", "1000"}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{"sum", "=", 
    RowBox[{"sum", "+", "i"}]}]}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["sum"], "Input"],

Cell[BoxData["500500"], "Output",
 CellChangeTimes->{3.5060574393125*^9, 3.506340797734375*^9, 
  3.50634084834375*^9}]
}, Open  ]],

Cell[TextData[StyleBox["Se puede eliminar la parte de incremento y \
posteriormente incrementar la variable de control dentro del cuerpo de \
ejecuci\[OAcute]n",
 FontFamily->"Arial"]], "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"sum", "=", "0"}], ";"}]], "Input"],

Cell[BoxData[
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "\[LessEqual]", "1000"}], ",", 
   RowBox[{
    RowBox[{"sum", "=", 
     RowBox[{"sum", "+", "i"}]}], ";", 
    RowBox[{"i", "=", 
     RowBox[{"i", "+", "1"}]}], ";"}]}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["sum"], "Input"],

Cell[BoxData["500500"], "Output",
 CellChangeTimes->{3.5060574466875*^9, 3.5063407978125*^9, 3.506340848375*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ejemplo: Sustituci\[OAcute]n regresiva", "Subtitle"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
  "El", " ", "proceso", " ", "de", " ", "sustituci\[OAcute]n", " ", 
   "regresiva", " ", "se", " ", "puede", " ", "programar", " ", "mediante", 
   " ", "una", " ", "sentencia", " ", "de", " ", "tipo", " ", "Do"}], ",", 
  " ", "While", ",", " ", "For", ",", " ", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{"sentencias", " ", "ciclicas"}], ")"}], ".", 
     "\[IndentingNewLine]", "La"}], " ", "ley", " ", "de", " ", "recurrencia",
     " ", 
    RowBox[{"es", ":", " ", 
     SubscriptBox["x", "k"]}]}], "=", 
   FractionBox[
    RowBox[{
     SubscriptBox["b", "k"], "-", 
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"j", "=", 
        RowBox[{"k", "+", "1"}]}], "n"], 
      RowBox[{
       SubscriptBox["a", "kj"], "*", 
       SubscriptBox["x", "j"]}]}]}], 
    SubscriptBox["a", "kk"]]}], ",", " ", 
  RowBox[{
   RowBox[{"v\[AAcute]lido", " ", "desde", " ", "k"}], "=", 
   RowBox[{
    RowBox[{
     RowBox[{"n", "..."}], "..."}], "1."}]}]}], "\[IndentingNewLine]", 
 RowBox[{"Implementaremos", " ", "este", " ", "algoritmo", " ", "con", " ", 
  "un", " ", 
  RowBox[{"ejemplo", "."}]}]}], "Text",
 FontFamily->"Arial"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"a", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"1", ",", "1", ",", "2"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0", ",", "1", ",", 
      RowBox[{"-", "1"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0", ",", "0", ",", 
      RowBox[{"-", "4"}]}], "}"}]}], "}"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"1", ",", "1", ",", "2"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", 
     RowBox[{"-", "1"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", 
     RowBox[{"-", "4"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.506057449140625*^9, 3.506340797828125*^9, 
  3.5063408484375*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"a", "//", "TableForm"}]], "Input"],

Cell[BoxData[
 TagBox[GridBox[{
    {"1", "1", "2"},
    {"0", "1", 
     RowBox[{"-", "1"}]},
    {"0", "0", 
     RowBox[{"-", "4"}]}
   },
   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" -> {}}],
  Function[BoxForm`e$, 
   TableForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.506057449875*^9, 3.506340797859375*^9, 
  3.506340848453125*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"b", "=", 
  RowBox[{"{", 
   RowBox[{"1", ",", "0", ",", 
    RowBox[{"-", "1"}]}], "}"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", "0", ",", 
   RowBox[{"-", "1"}]}], "}"}]], "Output",
 CellChangeTimes->{3.506057450859375*^9, 3.506340797890625*^9, 
  3.506340848484375*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n", "=", 
  RowBox[{"Length", "[", "b", "]"}]}]], "Input"],

Cell[BoxData["3"], "Output",
 CellChangeTimes->{3.506057451421875*^9, 3.506340797921875*^9, 
  3.506340848515625*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"x", "=", 
  RowBox[{"Table", "[", 
   RowBox[{"0", ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0", ",", "0", ",", "0"}], "}"}]], "Output",
 CellChangeTimes->{3.506057452125*^9, 3.506340797953125*^9, 
  3.506340848546875*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Con el comando For:", "Subsection"],

Cell[BoxData[
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"k", "=", "n"}], ",", 
   RowBox[{"k", "\[GreaterEqual]", "1"}], ",", 
   RowBox[{"k", "--"}], ",", 
   RowBox[{
    RowBox[{
     RowBox[{"x", "[", 
      RowBox[{"[", "k", "]"}], "]"}], "=", 
     FractionBox[
      RowBox[{
       RowBox[{"b", "[", 
        RowBox[{"[", "k", "]"}], "]"}], "-", 
       RowBox[{
        UnderoverscriptBox["\[Sum]", 
         RowBox[{"j", "=", 
          RowBox[{"k", "+", "1"}]}], "n"], 
        RowBox[{
         RowBox[{"a", "[", 
          RowBox[{"[", 
           RowBox[{"k", ",", "j"}], "]"}], "]"}], "*", 
         RowBox[{"x", "[", 
          RowBox[{"[", "j", "]"}], "]"}]}]}]}], 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"k", ",", "k"}], "]"}], "]"}]]}], ";"}]}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["x"], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox["1", "4"], ",", 
   FractionBox["1", "4"], ",", 
   FractionBox["1", "4"]}], "}"}]], "Output",
 CellChangeTimes->{3.50605745428125*^9, 3.506340797984375*^9, 
  3.506340848578125*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"a", ".", "x"}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", "0", ",", 
   RowBox[{"-", "1"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5060574551875*^9, 3.506340798015625*^9, 
  3.506340848609375*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Con el comando While", "Subsection"],

Cell["\<\
No es necesario utilizar un vector soluci\[OAcute]n x, se pueden ir \
almacenandio los resultados en el propio vector de t\[EAcute]rminos \
independientes b,as\[IAcute] a medida que se va resolviendo el problema \
pasamos de tener el vector b a tener la soluci\[OAcute]n del problema. Esto \
permite disminuir los requerimientos de memoria del problema, dado que el \
sistema de ecuaciones puede ser de gr\[AAcute]n dimensi\[OAcute]n.\
\>", "Text",
 FontFamily->"Arial"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"a", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"1", ",", "1", ",", "2"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0", ",", "1", ",", 
      RowBox[{"-", "1"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0", ",", "0", ",", 
      RowBox[{"-", "4"}]}], "}"}]}], "}"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"1", ",", "1", ",", "2"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", 
     RowBox[{"-", "1"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", 
     RowBox[{"-", "4"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.506057457171875*^9, 3.50634079803125*^9, 
  3.506340848625*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"b", "=", 
  RowBox[{"{", 
   RowBox[{"1", ",", "0", ",", 
    RowBox[{"-", "1"}]}], "}"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", "0", ",", 
   RowBox[{"-", "1"}]}], "}"}]], "Output",
 CellChangeTimes->{3.506057457765625*^9, 3.5063407980625*^9, 
  3.50634084865625*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n", "=", 
  RowBox[{"Length", "[", "b", "]"}]}]], "Input"],

Cell[BoxData["3"], "Output",
 CellChangeTimes->{3.50605745840625*^9, 3.50634079809375*^9, 
  3.506340848703125*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"k", "=", "n"}]], "Input"],

Cell[BoxData["3"], "Output",
 CellChangeTimes->{3.5060574589375*^9, 3.506340798109375*^9, 3.50634084875*^9}
  ]
}, Open  ]],

Cell[BoxData[
 RowBox[{"While", "[", 
  RowBox[{
   RowBox[{"k", "\[GreaterEqual]", "1"}], ",", 
   RowBox[{
    RowBox[{
     RowBox[{"b", "[", 
      RowBox[{"[", "k", "]"}], "]"}], "=", 
     FractionBox[
      RowBox[{
       RowBox[{"b", "[", 
        RowBox[{"[", "k", "]"}], "]"}], "-", 
       RowBox[{
        UnderoverscriptBox["\[Sum]", 
         RowBox[{"j", "=", 
          RowBox[{"k", "+", "1"}]}], "n"], 
        RowBox[{
         RowBox[{"a", "[", 
          RowBox[{"[", 
           RowBox[{"k", ",", "j"}], "]"}], "]"}], "*", 
         RowBox[{"b", "[", 
          RowBox[{"[", "j", "]"}], "]"}]}]}]}], 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"k", ",", "k"}], "]"}], "]"}]]}], ";", 
    RowBox[{"k", "=", 
     RowBox[{"k", "-", "1"}]}], ";"}]}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["b"], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox["1", "4"], ",", 
   FractionBox["1", "4"], ",", 
   FractionBox["1", "4"]}], "}"}]], "Output",
 CellChangeTimes->{3.50605746009375*^9, 3.506340798140625*^9, 
  3.50634084878125*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Con el comando Do", "Subsection"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"a", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"1", ",", "1", ",", "2"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0", ",", "1", ",", 
      RowBox[{"-", "1"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0", ",", "0", ",", 
      RowBox[{"-", "4"}]}], "}"}]}], "}"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"1", ",", "1", ",", "2"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", 
     RowBox[{"-", "1"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", 
     RowBox[{"-", "4"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5060574613125*^9, 3.506340798171875*^9, 
  3.5063408488125*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"b", "=", 
  RowBox[{"{", 
   RowBox[{"1", ",", "0", ",", 
    RowBox[{"-", "1"}]}], "}"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", "0", ",", 
   RowBox[{"-", "1"}]}], "}"}]], "Output",
 CellChangeTimes->{3.50605746184375*^9, 3.506340798203125*^9, 
  3.506340848828125*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n", "=", 
  RowBox[{"Length", "[", "b", "]"}]}]], "Input"],

Cell[BoxData["3"], "Output",
 CellChangeTimes->{3.506057462296875*^9, 3.50634079821875*^9, 
  3.506340848859375*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"Do", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"b", "[", 
     RowBox[{"[", "k", "]"}], "]"}], "=", 
    FractionBox[
     RowBox[{
      RowBox[{"b", "[", 
       RowBox[{"[", "k", "]"}], "]"}], "-", 
      RowBox[{
       UnderoverscriptBox["\[Sum]", 
        RowBox[{"j", "=", 
         RowBox[{"k", "+", "1"}]}], "n"], 
       RowBox[{
        RowBox[{"a", "[", 
         RowBox[{"[", 
          RowBox[{"k", ",", "j"}], "]"}], "]"}], "*", 
        RowBox[{"b", "[", 
         RowBox[{"[", "j", "]"}], "]"}]}]}]}], 
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"k", ",", "k"}], "]"}], "]"}]]}], ",", 
   RowBox[{"{", 
    RowBox[{"k", ",", "n", ",", "1", ",", 
     RowBox[{"-", "1"}]}], "}"}]}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["b"], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox["1", "4"], ",", 
   FractionBox["1", "4"], ",", 
   FractionBox["1", "4"]}], "}"}]], "Output",
 CellChangeTimes->{3.506057463375*^9, 3.506340798265625*^9, 
  3.506340848890625*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{1672, 933},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
FrontEndVersion->"8.0 for Microsoft Windows (32-bit) (November 7, 2010)",
StyleDefinitions->FrontEnd`FileName[{"Report"}, "StandardReport.nb", 
  CharacterEncoding -> "WindowsANSI"]
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{
 "Do"->{
  Cell[1859, 49, 86, 3, 40, "Input",
   CellTags->"Do"],
  Cell[1948, 54, 223, 8, 40, "Input",
   CellTags->"Do"]}
 }
*)
(*CellTagsIndex
CellTagsIndex->{
 {"Do", 19478, 726}
 }
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 53, 0, 81, "Title"],
Cell[635, 24, 868, 13, 83, "Text"],
Cell[CellGroupData[{
Cell[1528, 41, 32, 0, 29, "Subtitle"],
Cell[1563, 43, 293, 4, 65, "Text"],
Cell[1859, 49, 86, 3, 40, "Input",
 CellTags->"Do"],
Cell[1948, 54, 223, 8, 40, "Input",
 CellTags->"Do"],
Cell[CellGroupData[{
Cell[2196, 66, 29, 0, 40, "Input"],
Cell[2228, 68, 119, 2, 39, "Output"]
}, Open  ]],
Cell[2362, 73, 510, 16, 57, "Text"],
Cell[CellGroupData[{
Cell[2897, 93, 122, 4, 59, "Input"],
Cell[3022, 99, 120, 2, 39, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3179, 106, 109, 3, 64, "Input"],
Cell[3291, 111, 121, 2, 39, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[3461, 119, 35, 0, 29, "Subtitle"],
Cell[3499, 121, 301, 6, 65, "Text"],
Cell[3803, 129, 139, 4, 61, "Input"],
Cell[3945, 135, 254, 8, 40, "Input"],
Cell[CellGroupData[{
Cell[4224, 147, 29, 0, 40, "Input"],
Cell[4256, 149, 109, 1, 39, "Output"]
}, Open  ]],
Cell[4380, 153, 159, 4, 29, "Text"],
Cell[4542, 159, 139, 4, 61, "Input"],
Cell[4684, 165, 385, 12, 40, "Input"],
Cell[CellGroupData[{
Cell[5094, 181, 29, 0, 40, "Input"],
Cell[5126, 183, 118, 2, 39, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5281, 190, 27, 0, 40, "Input"],
Cell[5311, 192, 115, 2, 39, "Output"]
}, Open  ]],
Cell[5441, 197, 683, 16, 56, "Text"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6161, 218, 33, 0, 29, "Subtitle"],
Cell[6197, 220, 549, 7, 83, "Text"],
Cell[6749, 229, 69, 2, 40, "Input"],
Cell[6821, 233, 234, 7, 40, "Input"],
Cell[CellGroupData[{
Cell[7080, 244, 29, 0, 40, "Input"],
Cell[7112, 246, 119, 2, 39, "Output"]
}, Open  ]],
Cell[7246, 251, 193, 3, 29, "Text"],
Cell[7442, 256, 69, 2, 40, "Input"],
Cell[7514, 260, 286, 9, 40, "Input"],
Cell[CellGroupData[{
Cell[7825, 273, 29, 0, 40, "Input"],
Cell[7857, 275, 112, 1, 39, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[8018, 282, 58, 0, 29, "Subtitle"],
Cell[8079, 284, 1212, 35, 98, "Text"],
Cell[CellGroupData[{
Cell[9316, 323, 333, 11, 40, "Input"],
Cell[9652, 336, 389, 12, 39, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10078, 353, 58, 1, 40, "Input"],
Cell[10139, 356, 693, 21, 81, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10869, 382, 126, 4, 40, "Input"],
Cell[10998, 388, 191, 5, 39, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11226, 398, 82, 2, 40, "Input"],
Cell[11311, 402, 117, 2, 39, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11465, 409, 166, 5, 40, "Input"],
Cell[11634, 416, 169, 4, 39, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11840, 425, 41, 0, 24, "Subsection"],
Cell[11884, 427, 806, 26, 60, "Input"],
Cell[CellGroupData[{
Cell[12715, 457, 27, 0, 40, "Input"],
Cell[12745, 459, 237, 7, 56, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13019, 471, 49, 1, 40, "Input"],
Cell[13071, 474, 189, 5, 39, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[13309, 485, 42, 0, 24, "Subsection"],
Cell[13354, 487, 480, 8, 47, "Text"],
Cell[CellGroupData[{
Cell[13859, 499, 333, 11, 40, "Input"],
Cell[14195, 512, 387, 12, 39, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14619, 529, 126, 4, 40, "Input"],
Cell[14748, 535, 188, 5, 39, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14973, 545, 82, 2, 40, "Input"],
Cell[15058, 549, 115, 2, 39, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15210, 556, 49, 1, 40, "Input"],
Cell[15262, 559, 111, 2, 39, "Output"]
}, Open  ]],
Cell[15388, 564, 805, 26, 60, "Input"],
Cell[CellGroupData[{
Cell[16218, 594, 27, 0, 40, "Input"],
Cell[16248, 596, 236, 7, 56, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[16533, 609, 39, 0, 24, "Subsection"],
Cell[CellGroupData[{
Cell[16597, 613, 333, 11, 40, "Input"],
Cell[16933, 626, 387, 12, 39, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17357, 643, 126, 4, 40, "Input"],
Cell[17486, 649, 190, 5, 39, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17713, 659, 82, 2, 40, "Input"],
Cell[17798, 663, 116, 2, 39, "Output"]
}, Open  ]],
Cell[17929, 668, 754, 25, 60, "Input"],
Cell[CellGroupData[{
Cell[18708, 697, 27, 0, 40, "Input"],
Cell[18738, 699, 235, 7, 56, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)