(* 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[      7931,        211]
NotebookOptionsPosition[      7455,        191]
NotebookOutlinePosition[      7794,        206]
CellTagsIndexPosition[      7751,        203]
WindowFrame->Normal
ContainsDynamic->True *)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 StyleBox[
  RowBox[{"The", " ", "MVT", " ", "and", " ", "Approximations", " ", 
   RowBox[{"(", 
    RowBox[{"ie", ",", " ", 
     RowBox[{"Taylor", " ", "Series", " ", "Light"}]}], ")"}]}], 
  "Subtitle"]], "Input",
 CellChangeTimes->{{3.4341008031024427`*^9, 3.4341008174461927`*^9}}],

Cell[BoxData[
 StyleBox["Example", "Subsubtitle"]], "Input",
 CellChangeTimes->{{3.4341008272899427`*^9, 3.4341008280711927`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", "x_", "]"}], " ", ":=", " ", 
   RowBox[{"Sin", "[", "x", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"g", "[", "x_", "]"}], " ", ":=", " ", 
   RowBox[{"Cos", "[", "x", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"h", "[", "x_", "]"}], " ", ":=", " ", 
   RowBox[{"-", 
    RowBox[{"Sin", "[", "x", "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"L", "[", "x_", "]"}], " ", ":=", " ", "x"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Error", "[", "x_", "]"}], " ", ":=", " ", 
   RowBox[{
    RowBox[{"f", "[", "x", "]"}], " ", "-", " ", 
    RowBox[{"L", "[", "x", "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"x1Error", "[", "x_", "]"}], " ", ":=", " ", 
   RowBox[{"ArcCos", "[", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"Error", "[", "x", "]"}], "/", "x"}], ")"}], " ", "+", " ", 
     RowBox[{"g", "[", "0", "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"x2Error", "[", "x_", "]"}], " ", ":=", " ", 
    RowBox[{"-", " ", 
     RowBox[{"ArcSin", "[", " ", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"g", "[", 
          RowBox[{"x1Error", "[", "x", "]"}], "]"}], " ", "-", 
         RowBox[{"g", "[", "0", "]"}]}], ")"}], " ", "/", " ", 
       RowBox[{"x1Error", "[", "x", "]"}]}], "]"}]}]}], ";"}], 
  " "}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.4341008335243177`*^9, 3.4341009603524427`*^9}, {
  3.4341014518055677`*^9, 3.4341015046961927`*^9}, {3.4341018958680677`*^9, 
  3.4341019043368177`*^9}, {3.4341022266649427`*^9, 3.4341022419618177`*^9}}],\


Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
  "Print", "[", 
   "\"\<Comparison of Sin[x] and L[x] = x. Let Error[x] = Sin[x] - L[x]. The \
three points below (in decreasing order) are x, the x from the MVT applied to \
Error (where we evaluate Error'[x]), and the x from the MVT applied to Error' \
(where we evaluate Error''[x].\>\"", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"f", "[", "x", "]"}], ",", "x"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", 
       RowBox[{"Pi", "/", "2"}]}], "}"}], ",", 
     RowBox[{"Epilog", " ", "\[Rule]", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"PointSize", "[", ".025", "]"}], ",", " ", 
        RowBox[{"Point", "[", 
         RowBox[{"{", 
          RowBox[{"x0", ",", "0"}], "}"}], "]"}], ",", "  ", 
        RowBox[{"Point", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"x1Error", "[", "x0", "]"}], ",", "0"}], "}"}], "]"}], ",",
         "  ", 
        RowBox[{"Point", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"x2Error", "[", "x0", "]"}], ",", "0"}], "}"}], "]"}]}], 
       "}"}]}]}], "]"}], ",", " ", 
   RowBox[{"{", 
    RowBox[{"x0", ",", "0.000001", ",", 
     RowBox[{"Pi", "/", "2"}]}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.4341009303993177`*^9, 3.4341009358055677`*^9}, {
   3.4341009818368177`*^9, 3.4341010295711927`*^9}, {3.4341010615399427`*^9, 
   3.4341011514618177`*^9}, {3.4341011878836927`*^9, 
   3.4341011997118177`*^9}, {3.4341013158680677`*^9, 3.4341013338680677`*^9}, 
   3.4341015089461927`*^9, {3.4341015793368177`*^9, 3.4341015982430677`*^9}, {
   3.4341020025399427`*^9, 3.4341020029930677`*^9}, {3.4341021068836927`*^9, 
   3.4341021320086927`*^9}, {3.4341021786961927`*^9, 
   3.4341021856961927`*^9}, {3.4341022555555677`*^9, 
   3.4341022582118177`*^9}, {3.4341022968993177`*^9, 
   3.4341024811649427`*^9}, {3.4341025305555677`*^9, 
   3.4341025921024427`*^9}, {3.4341030540399427`*^9, 3.4341030850868177`*^9}}],\


Cell[BoxData["\<\"Comparison of Sin[x] and L[x] = x. Let Error[x] = Sin[x] - \
L[x]. The three points below (in decreasing order) are x, the x from the MVT \
applied to Error (where we evaluate Error'[x]), and the x from the MVT \
applied to Error' (where we evaluate Error''[x].\"\>"], "Print",
 CellChangeTimes->{
  3.4341025936961927`*^9, {3.4341030623993177`*^9, 3.4341030854149427`*^9}}],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`x0$$ = 1.0838497754884786`, Typeset`show$$ =
     True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", 
    Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = 
    "\"untitled\"", Typeset`specs$$ = {{
      Hold[$CellContext`x0$$], 1.*^-6, Rational[1, 2] Pi}}, Typeset`size$$ = {
    360., {112., 116.}}, Typeset`update$$ = 0, Typeset`initDone$$, 
    Typeset`skipInitDone$$ = True, $CellContext`x0$5644$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, "Variables" :> {$CellContext`x0$$ = 1.*^-6}, 
      "ControllerVariables" :> {
        Hold[$CellContext`x0$$, $CellContext`x0$5644$$, 0]}, 
      "OtherVariables" :> {
       Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
        Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
        Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
         Typeset`skipInitDone$$}, "Body" :> Plot[{
         $CellContext`f[$CellContext`x], $CellContext`x}, {$CellContext`x, 0, 
         Pi/2}, Epilog -> {
          PointSize[0.025], 
          Point[{$CellContext`x0$$, 0}], 
          Point[{
            $CellContext`x1Error[$CellContext`x0$$], 0}], 
          Point[{
            $CellContext`x2Error[$CellContext`x0$$], 0}]}], 
      "Specifications" :> {{$CellContext`x0$$, 1.*^-6, Rational[1, 2] Pi}}, 
      "Options" :> {}, "DefaultOptions" :> {}],
     ImageSizeCache->{408., {156., 161.}},
     SingleEvaluation->True],
    Deinitialization:>None,
    DynamicModuleValues:>{},
    SynchronousInitialization->True,
    UnsavedVariables:>{Typeset`initDone$$},
    UntrackedVariables:>{Typeset`size$$}], "Manipulate",
   Deployed->True,
   StripOnInput->False],
  Manipulate`InterpretManipulate[1]]], "Output",
 CellChangeTimes->{{3.4341021171336927`*^9, 3.4341021861180677`*^9}, {
   3.4341022308055677`*^9, 3.4341022591336927`*^9}, {3.4341023020555677`*^9, 
   3.4341023530399427`*^9}, 3.4341024822430677`*^9, 3.4341025316493177`*^9, 
   3.4341025937118177`*^9, {3.4341030624149427`*^9, 3.4341030854305677`*^9}}]
}, Open  ]]
},
WindowSize->{835, 658},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"6.0 for Microsoft Windows (32-bit) (June 19, 2007)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[568, 21, 301, 7, 44, "Input"],
Cell[872, 30, 130, 2, 34, "Input"],
Cell[1005, 34, 1768, 49, 172, "Input"],
Cell[CellGroupData[{
Cell[2800, 88, 2110, 48, 172, "Input"],
Cell[4915, 139, 392, 5, 62, "Print"],
Cell[5310, 146, 2129, 42, 334, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)