(* 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[     32843,        888]
NotebookOptionsPosition[     30840,        821]
NotebookOutlinePosition[     31292,        839]
CellTagsIndexPosition[     31249,        836]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
Math 481A
California State University, Northridge
Newton's Interpolation Method: Backward Differences\
\>", "Section"],

Cell["\[Copyright]2007 BE Shapiro", "Subsection",
 CellChangeTimes->{{3.398601216828899*^9, 3.398601217595446*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["BackwardInterpolate", "Section"],

Cell["\<\
Input: List of points of the form {{x0,f0}, {x1, f1}, ....} and a point x
Output: Interpolated value of x to the highest order possible, given the data

To print out the interpolation table, use NewtonInterpolationTable, which is \
otherwise identical to this function but has lots of printing.\
\>", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"BackwardInterpolate", "[", 
     RowBox[{"points_", ",", " ", "X_"}], "]"}], ":=", " ", 
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "n", ",", " ", "F", ",", " ", "i", ",", " ", "j", ",", 
        "\[CapitalDelta]F", ",", " ", "P", ",", "s", ",", "h", ",", " ", 
        "x"}], "}"}], ",", "\[IndentingNewLine]", " ", "\[IndentingNewLine]", 
      
      StyleBox[
       RowBox[{"(*", " ", 
        RowBox[{
        "points", " ", "should", " ", "have", " ", "the", " ", "form", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"x0", ",", " ", 
             RowBox[{"f", "[", "x0", "]"}]}], "}"}], ",", " ", 
           RowBox[{"{", 
            RowBox[{"x1", ",", " ", 
             RowBox[{"f", "[", "x1", "]"}]}], "}"}], ",", " ", "..."}], 
          "}"}]}], "*)"}],
       CellFrame->True,
       Background->RGBColor[1, 1, 0]], 
      StyleBox["\[IndentingNewLine]",
       CellFrame->True,
       Background->RGBColor[1, 1, 0]], 
      StyleBox["\[IndentingNewLine]",
       CellFrame->True,
       Background->RGBColor[1, 1, 0]], 
      StyleBox[
       RowBox[{"(*", " ", 
        RowBox[{
         RowBox[{"NOTE", ":", " ", 
          RowBox[{
          "if", " ", "the", " ", "xi", " ", "are", " ", "not", " ", "equally",
            " ", "spaced"}]}], ",", " ", 
         RowBox[{
          RowBox[{
          "this", " ", "will", " ", "not", " ", "calculate", " ", "the", " ", 
           "correct", " ", "answer"}], " ", "-", " ", 
          RowBox[{
          "ERROR", " ", "CHECKING", " ", "IS", " ", "NOT", " ", 
           "PERFORMED"}]}]}], " ", "*)"}],
       CellFrame->True,
       FontSize->18,
       Background->RGBColor[1, 1, 0]], 
      StyleBox[" ",
       CellFrame->True,
       FontSize->18,
       Background->RGBColor[1, 1, 0]], "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"n", "=", 
        RowBox[{"Length", "[", "points", "]"}]}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"x", ",", " ", "F"}], "}"}], " ", "=", " ", 
        RowBox[{"Transpose", "[", "points", "]"}]}], ";", " ", 
       "\[IndentingNewLine]", "\[IndentingNewLine]", 
       StyleBox[
        RowBox[{"(*", " ", 
         RowBox[{
          RowBox[{"The", " ", "variables", " ", "ddtable"}], ",", " ", 
          RowBox[{
           RowBox[{
           "nextline", " ", "are", " ", "purely", " ", "for", " ", "show"}], 
           " ", "-", " ", 
           RowBox[{
           "it", " ", "is", " ", "used", " ", "to", " ", "print", " ", "the", 
            " ", "table"}]}], ",", " ", 
          RowBox[{
          "but", " ", "is", " ", "not", " ", "otherwise", " ", 
           "referenced"}]}], " ", "*)"}],
        Background->RGBColor[1, 1, 0]], 
       StyleBox[" ",
        Background->GrayLevel[0.666667]], 
       StyleBox["\[IndentingNewLine]",
        Background->GrayLevel[0.666667]], "\[IndentingNewLine]", 
       StyleBox[
        RowBox[{"(*", " ", 
         RowBox[{
          RowBox[{
           RowBox[{"\[CapitalDelta]F", "[", 
            RowBox[{"0", ",", " ", "i"}], "]"}], " ", "represents", " ", 
           SubscriptBox["f", "i"]}], " ", "=", " ", 
          RowBox[{"f", 
           RowBox[{"(", 
            SubscriptBox["x", "i"], ")"}]}]}], "*)"}],
        Background->RGBColor[1, 1, 0]], 
       StyleBox[" ",
        Background->RGBColor[1, 1, 0]], 
       StyleBox["\[IndentingNewLine]",
        Background->RGBColor[1, 1, 0]], 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"i", "=", "1"}], ",", " ", 
         RowBox[{"i", "\[LessEqual]", " ", "n"}], ",", " ", 
         RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           RowBox[{"\[CapitalDelta]F", "[", 
            RowBox[{"0", ",", "i"}], "]"}], "=", 
           RowBox[{"F", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], ";"}]}], "\[IndentingNewLine]", 
        "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
       StyleBox[
        RowBox[{"(*", " ", 
         RowBox[{
          RowBox[{"\[CapitalDelta]F", "[", 
           RowBox[{"i", ",", " ", "j"}], "]"}], " ", "represents", " ", 
          SuperscriptBox["\[CapitalDelta]", "i"], 
          RowBox[{"(", 
           SubscriptBox["x", "j"], ")"}]}], "*)"}],
        Background->RGBColor[1, 1, 0]], 
       StyleBox[" ",
        Background->RGBColor[1, 1, 0]], 
       StyleBox["\[IndentingNewLine]",
        Background->RGBColor[1, 1, 0]], 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"i", "=", "1"}], ",", " ", 
         RowBox[{"i", "\[LessEqual]", " ", 
          RowBox[{"n", "-", "1"}]}], ",", 
         RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"For", "[", 
           RowBox[{
            RowBox[{"j", "=", "1"}], ",", " ", 
            RowBox[{"j", "\[LessEqual]", " ", 
             RowBox[{"n", "-", "i"}]}], ",", " ", 
            RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", " ", 
            RowBox[{
             RowBox[{
              RowBox[{"\[CapitalDelta]F", "[", 
               RowBox[{"i", ",", "j"}], "]"}], "=", " ", 
              RowBox[{
               RowBox[{"\[CapitalDelta]F", "[", 
                RowBox[{
                 RowBox[{"i", "-", "1"}], ",", " ", 
                 RowBox[{"j", "+", "1"}]}], "]"}], "-", 
               RowBox[{"\[CapitalDelta]F", "[", 
                RowBox[{
                 RowBox[{"i", "-", "1"}], ",", "j"}], "]"}]}]}], ";"}]}], 
           "\[IndentingNewLine]", "]"}], ";"}]}], " ", "\[IndentingNewLine]", 
        "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
       StyleBox[
        RowBox[{"(*", " ", 
         RowBox[{
          RowBox[{
          "Mathematica", " ", "indices", " ", "start", " ", "at", " ", "1"}], 
          ",", " ", 
          RowBox[{"not", " ", "zero"}]}], " ", "*)"}],
        Background->RGBColor[1, 1, 0]], 
       StyleBox["\[IndentingNewLine]",
        Background->RGBColor[1, 1, 0]], 
       RowBox[{"h", "=", 
        RowBox[{
         RowBox[{"x", "[", 
          RowBox[{"[", "2", "]"}], "]"}], "-", 
         RowBox[{"x", "[", 
          RowBox[{"[", "1", "]"}], "]"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"s", "=", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"X", "-", 
           RowBox[{"x", "[", 
            RowBox[{"[", "n", "]"}], "]"}]}], ")"}], "/", "h"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"P", "[", 
         RowBox[{"u_", ",", "j_"}], "]"}], ":=", 
        RowBox[{
         RowBox[{"\[CapitalDelta]F", "[", 
          RowBox[{"0", ",", "j"}], "]"}], "+", 
         RowBox[{
          UnderoverscriptBox["\[Sum]", 
           RowBox[{"i", "=", "1"}], 
           RowBox[{"j", "-", "1"}]], 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"-", "1"}], ")"}], "^", "i"}], " ", "*", " ", 
           RowBox[{"Binomial", "[", 
            RowBox[{
             RowBox[{"-", "u"}], ",", "i"}], "]"}], 
           RowBox[{"\[CapitalDelta]F", "[", 
            RowBox[{"i", ",", 
             RowBox[{"j", "-", "i"}]}], "]"}]}]}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"Return", "[", 
        RowBox[{"P", "[", 
         RowBox[{"s", ",", "n"}], "]"}], "]"}], ";"}]}], " ", 
     "\[IndentingNewLine]", "\[IndentingNewLine]", "]"}]}], ";"}], 
  " "}]], "Input"]
}, Closed]],

Cell[CellGroupData[{

Cell["BackwardInterpolationTable", "Section"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"BackwardInterpolationTable", "[", 
     RowBox[{"points_", ",", " ", "X_"}], "]"}], ":=", " ", 
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "n", ",", " ", "F", ",", " ", "i", ",", " ", "j", ",", 
        "\[CapitalDelta]F", ",", " ", "P", ",", "t", ",", "h", ",", " ", "x", 
        ",", " ", "ddtable", ",", " ", "nextline"}], "}"}], ",", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      StyleBox[
       RowBox[{"(*", " ", 
        RowBox[{
        "points", " ", "should", " ", "have", " ", "the", " ", "form", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"x0", ",", " ", 
             RowBox[{"f", "[", "x0", "]"}]}], "}"}], ",", " ", 
           RowBox[{"{", 
            RowBox[{"x1", ",", " ", 
             RowBox[{"f", "[", "x1", "]"}]}], "}"}], ",", " ", "..."}], 
          "}"}]}], "*)"}],
       Background->RGBColor[1, 1, 0]], 
      StyleBox["\[IndentingNewLine]",
       Background->RGBColor[1, 1, 0]], "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"n", "=", 
        RowBox[{"Length", "[", "points", "]"}]}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"x", ",", " ", "F"}], "}"}], " ", "=", " ", 
        RowBox[{"Transpose", "[", "points", "]"}]}], ";", " ", 
       "\[IndentingNewLine]", "\[IndentingNewLine]", 
       StyleBox[
        RowBox[{"(*", " ", 
         RowBox[{
          RowBox[{"The", " ", "variables", " ", "ddtable"}], ",", " ", 
          RowBox[{
           RowBox[{
           "nextline", " ", "are", " ", "purely", " ", "for", " ", "show"}], 
           " ", "-", " ", 
           RowBox[{
           "it", " ", "is", " ", "used", " ", "to", " ", "print", " ", "the", 
            " ", "table"}]}], ",", " ", 
          RowBox[{
          "but", " ", "is", " ", "not", " ", "otherwise", " ", 
           "referenced"}]}], " ", "*)"}],
        Background->RGBColor[1, 1, 0]], 
       StyleBox[" ",
        Background->RGBColor[1, 1, 0]], 
       StyleBox["\[IndentingNewLine]",
        Background->RGBColor[1, 1, 0]], 
       RowBox[{"ddtable", " ", "=", " ", 
        RowBox[{"{", 
         RowBox[{"x", ",", " ", "F"}], "}"}]}], ";", " ", 
       "\[IndentingNewLine]", "\[IndentingNewLine]", 
       StyleBox[
        RowBox[{"(*", " ", 
         RowBox[{
          RowBox[{
           RowBox[{"\[CapitalDelta]F", "[", 
            RowBox[{"0", ",", " ", "i"}], "]"}], " ", "represents", " ", 
           SubscriptBox["f", "i"]}], " ", "=", " ", 
          RowBox[{"f", 
           RowBox[{"(", 
            SubscriptBox["x", "i"], ")"}]}]}], "*)"}],
        Background->RGBColor[1, 1, 0]], 
       StyleBox[" ",
        Background->RGBColor[1, 1, 0]], 
       StyleBox["\[IndentingNewLine]",
        Background->RGBColor[1, 1, 0]], 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"i", "=", "1"}], ",", " ", 
         RowBox[{"i", "\[LessEqual]", " ", "n"}], ",", " ", 
         RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           RowBox[{"\[CapitalDelta]F", "[", 
            RowBox[{"0", ",", "i"}], "]"}], "=", 
           RowBox[{"F", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
          RowBox[{"Print", "[", 
           RowBox[{
           "\"\<\[CapitalDelta]F[0,\>\"", ",", "i", ",", "\"\<]=\>\"", ",", 
            RowBox[{"\[CapitalDelta]F", "[", 
             RowBox[{"0", ",", "i"}], "]"}]}], "]"}], ";"}]}], 
        "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       StyleBox[
        RowBox[{"(*", " ", 
         RowBox[{
          RowBox[{"\[CapitalDelta]F", "[", 
           RowBox[{"i", ",", " ", "j"}], "]"}], " ", "represents", " ", 
          SuperscriptBox["\[CapitalDelta]", "i"], 
          RowBox[{"(", 
           SubscriptBox["x", "j"], ")"}]}], "*)"}],
        Background->RGBColor[1, 1, 0]], 
       StyleBox[" ",
        Background->RGBColor[1, 1, 0]], 
       StyleBox["\[IndentingNewLine]",
        Background->RGBColor[1, 1, 0]], 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"i", "=", "1"}], ",", " ", 
         RowBox[{"i", "\[LessEqual]", " ", 
          RowBox[{"n", "-", "1"}]}], ",", 
         RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"nextline", "=", 
           RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
          RowBox[{"For", "[", 
           RowBox[{
            RowBox[{"j", "=", "1"}], ",", " ", 
            RowBox[{"j", "\[LessEqual]", " ", 
             RowBox[{"n", "-", "i"}]}], ",", " ", 
            RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", " ", 
            RowBox[{
             RowBox[{
              RowBox[{"\[CapitalDelta]F", "[", 
               RowBox[{"i", ",", "j"}], "]"}], "=", " ", 
              RowBox[{
               RowBox[{"\[CapitalDelta]F", "[", 
                RowBox[{
                 RowBox[{"i", "-", "1"}], ",", " ", 
                 RowBox[{"j", "+", "1"}]}], "]"}], "-", 
               RowBox[{"\[CapitalDelta]F", "[", 
                RowBox[{
                 RowBox[{"i", "-", "1"}], ",", "j"}], "]"}]}]}], ";", 
             "\[IndentingNewLine]", 
             RowBox[{"nextline", " ", "=", " ", 
              RowBox[{"Append", "[", 
               RowBox[{"nextline", ",", " ", 
                RowBox[{"\[CapitalDelta]F", "[", 
                 RowBox[{"i", ",", "j"}], "]"}]}], "]"}]}], ";", " ", 
             "\[IndentingNewLine]", 
             RowBox[{"Print", "[", 
              RowBox[{
              "\"\<\[CapitalDelta]F[\>\"", ",", "i", ",", "\"\<,\>\"", ",", 
               "j", ",", "\"\<]=\>\"", ",", 
               RowBox[{"\[CapitalDelta]F", "[", 
                RowBox[{"i", ",", "j"}], "]"}]}], "]"}], ";"}]}], 
           "\[IndentingNewLine]", "]"}], ";", " ", "\[IndentingNewLine]", 
          RowBox[{"ddtable", " ", "=", " ", 
           RowBox[{"Append", "[", 
            RowBox[{"ddtable", ",", " ", 
             RowBox[{"PadRight", "[", 
              RowBox[{"nextline", ",", "n", ",", "\"\< \>\""}], "]"}]}], 
            "]"}]}], ";"}]}], " ", "\[IndentingNewLine]", "]"}], ";", 
       "\[IndentingNewLine]", "\[IndentingNewLine]", 
       StyleBox[
        RowBox[{"(*", " ", 
         RowBox[{
          RowBox[{
          "format", " ", "the", " ", "table", " ", "for", " ", "printing"}], 
          " ", "-", " ", 
          RowBox[{
          "this", " ", "part", " ", "is", " ", "not", " ", "needed", " ", 
           "for", " ", "the", " ", "algorithm"}]}], " ", "*)"}],
        Background->RGBColor[1, 1, 0]], 
       StyleBox["\[IndentingNewLine]",
        Background->RGBColor[1, 1, 0]], 
       RowBox[{"ddtable", " ", "=", " ", 
        RowBox[{"Transpose", "[", "ddtable", "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"nextline", " ", "=", 
        RowBox[{"Join", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
           "\"\<\!\(\*SubscriptBox[\(x\), \(k\)]\)\>\"", ",", 
            "\"\<f(\!\(\*SubscriptBox[\(x\), \(k\)]\))\>\""}], "}"}], ",", 
          "  ", 
          RowBox[{"\"\<\[CapitalDelta]\>\"", "^", 
           RowBox[{"ToString", "/@", 
            RowBox[{"Range", "[", 
             RowBox[{"1", ",", 
              RowBox[{"n", "-", "1"}]}], "]"}]}]}]}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"ddtable", " ", "=", " ", 
        RowBox[{"Prepend", "[", 
         RowBox[{"ddtable", ",", " ", "nextline"}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"Print", "[", 
        RowBox[{"ddtable", "//", "TableForm"}], "]"}], ";", 
       "\[IndentingNewLine]", "\[IndentingNewLine]", 
       StyleBox[
        RowBox[{"(*", " ", 
         RowBox[{
          RowBox[{
          "Mathematica", " ", "indices", " ", "start", " ", "at", " ", "1"}], 
          ",", " ", 
          RowBox[{"not", " ", "zero"}]}], " ", "*)"}],
        Background->RGBColor[1, 1, 0]], 
       StyleBox["\[IndentingNewLine]",
        Background->RGBColor[1, 1, 0]], 
       RowBox[{"h", "=", 
        RowBox[{
         RowBox[{"x", "[", 
          RowBox[{"[", "2", "]"}], "]"}], "-", 
         RowBox[{"x", "[", 
          RowBox[{"[", "1", "]"}], "]"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"t", "=", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"X", "-", 
           RowBox[{"x", "[", 
            RowBox[{"[", "n", "]"}], "]"}]}], ")"}], "/", "h"}]}], ";", 
       "\[IndentingNewLine]", "\[IndentingNewLine]", 
       RowBox[{"Print", "[", 
        RowBox[{"\"\<h=\>\"", ",", "h", ",", "\"\< t=\>\"", ",", "t"}], "]"}],
        ";", " ", "\[IndentingNewLine]", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"P", "[", 
         RowBox[{"u_", ",", "j_"}], "]"}], ":=", 
        RowBox[{
         RowBox[{"\[CapitalDelta]F", "[", 
          RowBox[{"0", ",", "n"}], "]"}], "+", 
         RowBox[{
          UnderoverscriptBox["\[Sum]", 
           RowBox[{"i", "=", "1"}], 
           RowBox[{"j", "-", "1"}]], 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"-", "1"}], ")"}], "^", "i"}], " ", "*", " ", 
           RowBox[{"Binomial", "[", 
            RowBox[{
             RowBox[{"-", "u"}], ",", "i"}], "]"}], 
           RowBox[{"\[CapitalDelta]F", "[", 
            RowBox[{"i", ",", 
             RowBox[{"n", "-", "i"}]}], "]"}]}]}]}]}], ";", 
       "\[IndentingNewLine]", "\[IndentingNewLine]", 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"i", "=", "1"}], ",", " ", 
         RowBox[{"i", "\[LessEqual]", " ", "n"}], ",", 
         RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"Print", "[", 
           RowBox[{
           "\"\<P[\>\"", ",", "X", ",", "\"\<] to \>\"", ",", "i", ",", 
            RowBox[{"Switch", "[", 
             RowBox[{
             "i", ",", "1", ",", "\"\<st\>\"", ",", "2", ",", "\"\<d\>\"", 
              ",", "3", ",", "\"\<rd\>\"", ",", "_", ",", " ", "\"\<th\>\""}],
              "]"}], ",", "\"\< order = \>\"", ",", " ", 
            RowBox[{"P", "[", 
             RowBox[{"t", ",", "i"}], "]"}]}], "]"}], ";"}]}], 
        "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{"Return", "[", 
        RowBox[{"P", "[", 
         RowBox[{"t", ",", "n"}], "]"}], "]"}], ";"}]}], " ", 
     "\[IndentingNewLine]", "\[IndentingNewLine]", "]"}]}], ";"}], 
  " "}]], "Input"]
}, Closed]],

Cell[CellGroupData[{

Cell["Examples", "Section"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"points", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"1", ",", "0.7651977"}], "}"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"1.3", ",", " ", "0.6200860"}], "}"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"1.6", ",", " ", "0.4554022"}], "}"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"1.9", ",", " ", "0.2818186"}], "}"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"2.2", ",", " ", "0.1103623"}], "}"}]}], "}"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"1", ",", "0.7651977`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1.3`", ",", "0.620086`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1.6`", ",", "0.4554022`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1.9`", ",", "0.2818186`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"2.2`", ",", "0.1103623`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.398601172117259*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"BackwardInterpolate", "[", 
  RowBox[{"points", ",", " ", "2.0"}], "]"}]], "Input"],

Cell[BoxData["0.22387536460905355`"], "Output",
 CellChangeTimes->{3.398601174883074*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"BackwardInterpolationTable", "[", 
  RowBox[{"points", ",", " ", "2."}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[0,\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", "0.7651977`"}],
  SequenceForm["\[CapitalDelta]F[0,", 1, "]=", 0.7651977],
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178398138*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[0,\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", "0.620086`"}],
  SequenceForm["\[CapitalDelta]F[0,", 2, "]=", 0.620086],
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178400476*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[0,\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", "0.4554022`"}],
  SequenceForm["\[CapitalDelta]F[0,", 3, "]=", 0.4554022],
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178403534*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[0,\"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", "0.2818186`"}],
  SequenceForm["\[CapitalDelta]F[0,", 4, "]=", 0.2818186],
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178404935*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[0,\"\>", "\[InvisibleSpace]", "5", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", "0.1103623`"}],
  SequenceForm["\[CapitalDelta]F[0,", 5, "]=", 0.1103623],
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178406321*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\",\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "0.14511169999999995`"}]}],
  SequenceForm["\[CapitalDelta]F[", 1, ",", 1, "]=", -0.14511169999999995`],
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178407702*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\",\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "0.16468380000000005`"}]}],
  SequenceForm["\[CapitalDelta]F[", 1, ",", 2, "]=", -0.16468380000000005`],
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178409047*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\",\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "0.1735836`"}]}],
  SequenceForm["\[CapitalDelta]F[", 1, ",", 3, "]=", -0.1735836],
  Editable->False]], "Print",
 CellChangeTimes->{3.3986011784102783`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\",\"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "0.17145629999999998`"}]}],
  SequenceForm["\[CapitalDelta]F[", 1, ",", 4, "]=", -0.17145629999999998`],
  Editable->False]], "Print",
 CellChangeTimes->{3.3986011784115*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\",\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "0.019572100000000092`"}]}],
  SequenceForm["\[CapitalDelta]F[", 2, ",", 1, "]=", -0.019572100000000092`],
  Editable->False]], "Print",
 CellChangeTimes->{3.3986011784126673`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\",\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "0.008899799999999958`"}]}],
  SequenceForm["\[CapitalDelta]F[", 2, ",", 2, "]=", -0.008899799999999958],
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178413765*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\",\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", 
   "0.0021273000000000264`"}],
  SequenceForm["\[CapitalDelta]F[", 2, ",", 3, "]=", 0.0021273000000000264`],
  Editable->False]], "Print",
 CellChangeTimes->{3.3986011784148684`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\",\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", 
   "0.010672300000000134`"}],
  SequenceForm["\[CapitalDelta]F[", 3, ",", 1, "]=", 0.010672300000000134`],
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178415955*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\",\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", 
   "0.011027099999999984`"}],
  SequenceForm["\[CapitalDelta]F[", 3, ",", 2, "]=", 0.011027099999999984`],
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178417038*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[CapitalDelta]F[\"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\",\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\"]=\"\>", "\[InvisibleSpace]", 
   "0.0003547999999998497`"}],
  SequenceForm["\[CapitalDelta]F[", 4, ",", 1, "]=", 0.0003547999999998497],
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178418126*^9}],

Cell[BoxData[
 TagBox[GridBox[{
    {"\<\"\\!\\(\\*SubscriptBox[\\(x\\), \\(k\\)]\\)\"\>", \
"\<\"f(\\!\\(\\*SubscriptBox[\\(x\\), \\(k\\)]\\))\"\>", 
     SuperscriptBox["\<\"\[CapitalDelta]\"\>", "\<\"1\"\>"], 
     SuperscriptBox["\<\"\[CapitalDelta]\"\>", "\<\"2\"\>"], 
     SuperscriptBox["\<\"\[CapitalDelta]\"\>", "\<\"3\"\>"], 
     SuperscriptBox["\<\"\[CapitalDelta]\"\>", "\<\"4\"\>"]},
    {"1", "0.7651977`", 
     RowBox[{"-", "0.14511169999999995`"}], 
     RowBox[{"-", "0.019572100000000092`"}], "0.010672300000000134`", 
     "0.0003547999999998497`"},
    {"1.3`", "0.620086`", 
     RowBox[{"-", "0.16468380000000005`"}], 
     RowBox[{"-", "0.008899799999999958`"}], 
     "0.011027099999999984`", "\<\" \"\>"},
    {"1.6`", "0.4554022`", 
     RowBox[{"-", "0.1735836`"}], 
     "0.0021273000000000264`", "\<\" \"\>", "\<\" \"\>"},
    {"1.9`", "0.2818186`", 
     RowBox[{
     "-", "0.17145629999999998`"}], "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
    {"2.2`", 
     "0.1103623`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"}
   },
   GridBoxAlignment->{
    "Columns" -> {{Left}}, "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" -> {}}],
  Function[BoxForm`e$, 
   TableForm[BoxForm`e$]]]], "Print",
 CellChangeTimes->{3.3986011784195023`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"h=\"\>", "\[InvisibleSpace]", "0.30000000000000004`", 
   "\[InvisibleSpace]", "\<\" t=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "0.6666666666666672`"}]}],
  SequenceForm["h=", 0.30000000000000004`, " t=", -0.6666666666666672],
  Editable->False]], "Print",
 CellChangeTimes->{3.3986011784214973`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"P[\"\>", "\[InvisibleSpace]", "2.`", 
   "\[InvisibleSpace]", "\<\"] to \"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\"st\"\>", 
   "\[InvisibleSpace]", "\<\" order = \"\>", "\[InvisibleSpace]", 
   "0.1103623`"}],
  SequenceForm["P[", 2., "] to ", 1, "st", " order = ", 0.1103623],
  Editable->False]], "Print",
 CellChangeTimes->{3.3986011784225893`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"P[\"\>", "\[InvisibleSpace]", "2.`", 
   "\[InvisibleSpace]", "\<\"] to \"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\"d\"\>", "\[InvisibleSpace]", "\<\" order = \"\>",
    "\[InvisibleSpace]", "0.22466650000000005`"}],
  SequenceForm["P[", 2., "] to ", 2, "d", " order = ", 0.22466650000000005`],
  Editable->False]], "Print",
 CellChangeTimes->{3.3986011784237213`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"P[\"\>", "\[InvisibleSpace]", "2.`", 
   "\[InvisibleSpace]", "\<\"] to \"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\"rd\"\>", 
   "\[InvisibleSpace]", "\<\" order = \"\>", "\[InvisibleSpace]", 
   "0.2244301333333334`"}],
  SequenceForm["P[", 2., "] to ", 3, "rd", " order = ", 0.2244301333333334],
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178424843*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"P[\"\>", "\[InvisibleSpace]", "2.`", 
   "\[InvisibleSpace]", "\<\"] to \"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\"th\"\>", 
   "\[InvisibleSpace]", "\<\" order = \"\>", "\[InvisibleSpace]", 
   "0.22388558518518523`"}],
  SequenceForm["P[", 2., "] to ", 4, "th", " order = ", 0.22388558518518523`],
  
  Editable->False]], "Print",
 CellChangeTimes->{3.39860117842597*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"P[\"\>", "\[InvisibleSpace]", "2.`", 
   "\[InvisibleSpace]", "\<\"] to \"\>", "\[InvisibleSpace]", "5", 
   "\[InvisibleSpace]", "\<\"th\"\>", 
   "\[InvisibleSpace]", "\<\" order = \"\>", "\[InvisibleSpace]", 
   "0.22387536460905355`"}],
  SequenceForm["P[", 2., "] to ", 5, "th", " order = ", 0.22387536460905355`],
  
  Editable->False]], "Print",
 CellChangeTimes->{3.398601178427105*^9}]
}, Open  ]],

Cell[BoxData["0.22387536460905355`"], "Output",
 CellChangeTimes->{3.398601178428193*^9}]
}, Open  ]]
}, Open  ]]
},
WindowToolbars->"RulerBar",
WindowSize->{811, 955},
WindowMargins->{{34, Automatic}, {Automatic, 11}},
Magnification->1,
FrontEndVersion->"6.0 for Mac OS X x86 (32-bit) (June 19, 2007)",
StyleDefinitions->FrontEnd`FileName[{"Creative"}, "PastelColor.nb", 
  CharacterEncoding -> "UTF-8"]
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[590, 23, 128, 4, 112, "Section"],
Cell[721, 29, 115, 1, 34, "Subsection"]
}, Open  ]],
Cell[CellGroupData[{
Cell[873, 35, 38, 0, 72, "Section"],
Cell[914, 37, 318, 6, 71, "Text"],
Cell[1235, 45, 7565, 200, 614, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[8837, 250, 45, 0, 42, "Section"],
Cell[8885, 252, 10623, 261, 866, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[19545, 518, 27, 0, 42, "Section"],
Cell[CellGroupData[{
Cell[19597, 522, 474, 13, 53, "Input"],
Cell[20074, 537, 446, 13, 37, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[20557, 555, 107, 2, 37, "Input"],
Cell[20667, 559, 89, 1, 37, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[20793, 565, 113, 2, 37, "Input"],
Cell[CellGroupData[{
Cell[20931, 571, 307, 6, 20, "Print"],
Cell[21241, 579, 305, 6, 20, "Print"],
Cell[21549, 587, 307, 6, 20, "Print"],
Cell[21859, 595, 307, 6, 20, "Print"],
Cell[22169, 603, 307, 6, 20, "Print"],
Cell[22479, 611, 416, 8, 20, "Print"],
Cell[22898, 621, 416, 8, 20, "Print"],
Cell[23317, 631, 397, 8, 20, "Print"],
Cell[23717, 641, 414, 8, 20, "Print"],
Cell[24134, 651, 420, 8, 20, "Print"],
Cell[24557, 661, 417, 8, 20, "Print"],
Cell[24977, 671, 406, 8, 20, "Print"],
Cell[25386, 681, 402, 8, 20, "Print"],
Cell[25791, 691, 402, 8, 20, "Print"],
Cell[26196, 701, 403, 8, 20, "Print"],
Cell[26602, 711, 1556, 37, 105, "Print"],
Cell[28161, 750, 353, 7, 20, "Print"],
Cell[28517, 759, 421, 9, 20, "Print"],
Cell[28941, 770, 436, 8, 20, "Print"],
Cell[29380, 780, 437, 9, 20, "Print"],
Cell[29820, 791, 442, 10, 20, "Print"],
Cell[30265, 803, 443, 10, 20, "Print"]
}, Open  ]],
Cell[30723, 816, 89, 1, 37, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)