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

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

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     12331,        358]
NotebookOptionsPosition[     11511,        325]
NotebookOutlinePosition[     12007,        345]
CellTagsIndexPosition[     11964,        342]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Exercise 4.49(a)", "Subsection",
 CellChangeTimes->{{3.608792515645569*^9, 3.6087925295123796`*^9}, {
  3.608816519920402*^9, 3.608816522571476*^9}, {3.608913275869972*^9, 
  3.6089132760338707`*^9}, {3.608913363716382*^9, 3.608913364226076*^9}, {
  3.609412065906448*^9, 3.609412078752801*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"R", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "\[Lambda]"}], ",", "\[Lambda]"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"\[Mu]", ",", 
       RowBox[{"-", "\[Mu]"}]}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "R", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dim", "=", 
   RowBox[{
    RowBox[{"MatrixRank", "[", "R", "]"}], "+", "1"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"um", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"1", ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "dim"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vec", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"i", "\[Equal]", "dim"}], ",", "1", ",", "0"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "dim"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", 
  RowBox[{"LinearSolve", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"Transpose", "[", "R", "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}], ",", "um"}], "}"}], ",", "vec"}], 
   "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.609412097962874*^9, 3.609412159731246*^9}, {
  3.6094122171699743`*^9, 3.609412225754033*^9}, {3.6094123118834333`*^9, 
  3.609412391835557*^9}, {3.609412642844348*^9, 3.609412720273533*^9}, {
  3.609412849657915*^9, 3.609412875813037*^9}, {3.609412945538601*^9, 
  3.609412949641547*^9}, {3.609413047019086*^9, 3.609413047938262*^9}, {
  3.609413095551879*^9, 3.6094131188512373`*^9}, {3.6094136994218483`*^9, 
  3.609413701411078*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", "\[Lambda]"}], "\[Lambda]"},
     {"\[Mu]", 
      RowBox[{"-", "\[Mu]"}]}
    },
    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.6088165701246023`*^9, 3.608816579301999*^9}, {
   3.60881706018918*^9, 3.608817084376479*^9}, 3.608817507300185*^9, {
   3.6088180484721327`*^9, 3.608818107518712*^9}, 3.608821047314049*^9, 
   3.608913246999919*^9, {3.608913422852075*^9, 3.608913448478587*^9}, 
   3.6089135477062683`*^9, {3.60941212449616*^9, 3.609412178624639*^9}, {
   3.6094123216073008`*^9, 3.6094123512240257`*^9}, {3.6094123872727337`*^9, 
   3.609412392252173*^9}, {3.609412709588524*^9, 3.609412721165134*^9}, {
   3.6094128550646*^9, 3.609412876670982*^9}, 3.609412997841967*^9, 
   3.609413048461776*^9, {3.609413096786017*^9, 3.60941311994621*^9}, 
   3.6094137019862413`*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", 
   TagBox[GridBox[{
      {
       FractionBox["\[Mu]", 
        RowBox[{"\[Lambda]", "+", "\[Mu]"}]]},
      {
       FractionBox["\[Lambda]", 
        RowBox[{"\[Lambda]", "+", "\[Mu]"}]]}
     },
     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.6088165701246023`*^9, 3.608816579301999*^9}, {
   3.60881706018918*^9, 3.608817084376479*^9}, 3.608817507300185*^9, {
   3.6088180484721327`*^9, 3.608818107518712*^9}, 3.608821047314049*^9, 
   3.608913246999919*^9, {3.608913422852075*^9, 3.608913448478587*^9}, 
   3.6089135477062683`*^9, {3.60941212449616*^9, 3.609412178624639*^9}, {
   3.6094123216073008`*^9, 3.6094123512240257`*^9}, {3.6094123872727337`*^9, 
   3.609412392252173*^9}, {3.609412709588524*^9, 3.609412721165134*^9}, {
   3.6094128550646*^9, 3.609412876670982*^9}, 3.609412997841967*^9, 
   3.609413048461776*^9, {3.609413096786017*^9, 3.60941311994621*^9}, 
   3.609413701989644*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Exercise 4.49(b)", "Subsection",
 CellChangeTimes->{{3.608792515645569*^9, 3.6087925295123796`*^9}, {
  3.608816519920402*^9, 3.608816522571476*^9}, {3.608913275869972*^9, 
  3.6089132760338707`*^9}, {3.608913363716382*^9, 3.608913364226076*^9}, {
  3.609412065906448*^9, 3.609412066752741*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"R", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["\[Lambda]", "1"]}], ",", "0", ",", 
       SubscriptBox["\[Lambda]", "1"], ",", "0", ",", "0"}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"-", 
        SubscriptBox["\[Lambda]", "2"]}], ",", "0", ",", 
       SubscriptBox["\[Lambda]", "2"], ",", "0"}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       SubscriptBox["\[Mu]", "1"], ",", 
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{
          SubscriptBox["\[Mu]", "1"], "+", 
          SubscriptBox["\[Lambda]", "2"]}], ")"}]}], ",", " ", "0", ",", " ", 
       
       SubscriptBox["\[Lambda]", "2"]}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["\[Mu]", "2"], ",", "0", ",", "0", ",", 
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{
          SubscriptBox["\[Mu]", "2"], "+", 
          SubscriptBox["\[Lambda]", "1"]}], ")"}]}], ",", 
       SubscriptBox["\[Lambda]", "1"]}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", 
       SubscriptBox["\[Mu]", "2"], ",", 
       SubscriptBox["\[Mu]", "1"], ",", 
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{
          SubscriptBox["\[Mu]", "1"], "+", 
          SubscriptBox["\[Mu]", "2"]}], ")"}]}]}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "R", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["\[Lambda]", "1"], "=", "1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["\[Lambda]", "2"], "=", "2"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["\[Mu]", "1"], "=", "0.1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["\[Mu]", "2"], "=", "0.15"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dim", "=", 
   RowBox[{
    RowBox[{"MatrixRank", "[", "R", "]"}], "+", "1"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"um", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"1", ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "dim"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vec", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"i", "\[Equal]", "dim"}], ",", "1", ",", "0"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "dim"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", 
  RowBox[{"LinearSolve", "[", 
   RowBox[{
    RowBox[{"Append", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"Transpose", "[", "R", "]"}], "[", 
         RowBox[{"[", "i", "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", 
          RowBox[{"dim", "-", "1"}]}], "}"}]}], "]"}], ",", "um"}], "]"}], 
    ",", "vec"}], "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.609412097962874*^9, 3.609412159731246*^9}, {
  3.6094122171699743`*^9, 3.609412225754033*^9}, {3.6094123118834333`*^9, 
  3.609412391835557*^9}, {3.609412642844348*^9, 3.609412720273533*^9}, {
  3.609412849657915*^9, 3.609412875813037*^9}, {3.609412945538601*^9, 
  3.609412949641547*^9}, {3.609413047019086*^9, 3.609413047938262*^9}, {
  3.609413095551879*^9, 3.609413142179937*^9}, {3.609413179357265*^9, 
  3.609413446073924*^9}, {3.609413519122067*^9, 3.609413604706319*^9}, {
  3.609413646108212*^9, 3.609413692218253*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", "1"}], "0", "1", "0", "0"},
     {"0", 
      RowBox[{"-", "2"}], "0", "2", "0"},
     {"0", "0.1`", 
      RowBox[{"-", "2.1`"}], "0", "2"},
     {"0.15`", "0", "0", 
      RowBox[{"-", "1.15`"}], "1"},
     {"0", "0", "0.15`", "0.1`", 
      RowBox[{"-", "0.25`"}]}
    },
    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.609413396318708*^9, {3.609413429617772*^9, 3.609413446602715*^9}, {
   3.60941351971349*^9, 3.609413563166078*^9}, 3.609413605701542*^9, {
   3.609413663175281*^9, 3.609413705310433*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", 
   TagBox[GridBox[{
      {"0.011818778726198308`"},
      {"0.0032829940906106366`"},
      {"0.06565988181221273`"},
      {"0.0787918581746554`"},
      {"0.840446487196323`"}
     },
     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.609413396318708*^9, {3.609413429617772*^9, 3.609413446602715*^9}, {
   3.60941351971349*^9, 3.609413563166078*^9}, 3.609413605701542*^9, {
   3.609413663175281*^9, 3.6094137053338223`*^9}}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{912, 773},
WindowMargins->{{451, Automatic}, {Automatic, 22}},
PrintingCopies->1,
PrintingPageRange->{1, Automatic},
PrivateNotebookOptions -> {
 "VersionedStylesheet" -> {"Default.nb"[8.] -> True}},
Magnification->2,
FrontEndVersion->"7.0 for Mac OS X x86 (32-bit) (November 10, 2008)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 301, 4, 68, "Subsection"],
Cell[CellGroupData[{
Cell[905, 30, 1739, 50, 257, "Input"],
Cell[2647, 82, 1332, 29, 105, "Output"],
Cell[3982, 113, 1450, 33, 137, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[5481, 152, 301, 4, 68, "Subsection"],
Cell[CellGroupData[{
Cell[5807, 160, 3637, 102, 557, "Input"],
Cell[9447, 264, 1058, 29, 199, "Output"],
Cell[10508, 295, 975, 26, 199, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)