(* 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[     21141,        695]
NotebookOptionsPosition[     18346,        598]
NotebookOutlinePosition[     18832,        617]
CellTagsIndexPosition[     18789,        614]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{
      "Compute", " ", "the", " ", "densities", " ", "and", " ", "nodes", " ", 
       "\[IndentingNewLine]", "for", " ", "Gaussian", " ", "quadrature", " ", 
       "\[IndentingNewLine]", "of", " ", "degree", " ", 
       RowBox[{"(", "precision", ")"}], "2", "n"}], "+", 
      RowBox[{"1", " ", 
       RowBox[{"over", " ", "[", 
        RowBox[{"a", ",", "b"}], "]"}]}]}], ";", "\[IndentingNewLine]", "  ", 
     
     RowBox[{
      RowBox[{"generates", " ", "n"}], "+", 
      RowBox[{"1", " ", "nodes", " ", "and", " ", 
       RowBox[{"densities", ".", "\[IndentingNewLine]", "Based"}], " ", "on", 
       " ", 
       RowBox[{"R", ".", "Curto"}], " ", "and", " ", 
       RowBox[{"L", ".", 
        RowBox[{"Fialkow", "\[IndentingNewLine]", "[", 
         RowBox[{
          RowBox[{"Houston", " ", 
           RowBox[{"J", ".", " ", "Math", ".", " ", "17"}], 
           RowBox[{"(", "1991", ")"}]}], ",", " ", 
          RowBox[{"603", "-", "635"}]}], "]"}]}]}]}]}], ",", 
    "\[IndentingNewLine]", "      ", 
    RowBox[{"MR", " ", "93", 
     RowBox[{"a", ":", "47016"}]}]}], "*)"}], "        ", 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Gaussian", "[", 
    RowBox[{"n_", ",", "a_", ",", "b_"}], "]"}], ":=", 
   RowBox[{"(", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "0"}], ",", 
       RowBox[{"i", "\[LessEqual]", 
        RowBox[{
         RowBox[{"2", "n"}], "+", "1"}]}], ",", 
       RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"bb", "[", "i", "]"}], " ", "=", " ", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"b", "^", 
            RowBox[{"(", 
             RowBox[{"i", "+", "1"}], ")"}]}], "-", 
           RowBox[{"a", "^", 
            RowBox[{"(", 
             RowBox[{"i", "+", "1"}], ")"}]}]}], ")"}], "/", 
         RowBox[{"(", 
          RowBox[{"i", "+", "1"}], ")"}]}]}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "0"}], ",", 
       RowBox[{"i", "\[LessEqual]", "n"}], ",", 
       RowBox[{"i", "=", 
        RowBox[{"i", "+", "1"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"j", "=", "0"}], ",", 
         RowBox[{"j", "\[LessEqual]", "n"}], ",", 
         RowBox[{"j", "=", 
          RowBox[{"j", "+", "1"}]}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"Mtab", "[", 
           RowBox[{"i", ",", "j"}], "]"}], "=", 
          RowBox[{"bb", "[", 
           RowBox[{"i", "+", "j"}], "]"}]}]}], "]"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"H", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Mtab", "[", 
         RowBox[{"r", ",", "s"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"r", ",", "0", ",", "n"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"s", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"J", " ", "=", " ", 
      RowBox[{"N", "[", 
       RowBox[{"Inverse", "[", "H", "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "0"}], ",", " ", 
       RowBox[{"i", "\[LessEqual]", "n"}], ",", 
       RowBox[{"i", "=", 
        RowBox[{"i", "+", "1"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"v", "[", "i", "]"}], " ", "=", " ", 
        RowBox[{"bb", "[", 
         RowBox[{"n", "+", "1", "+", "i"}], "]"}]}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"w", " ", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"v", "[", "i", "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"z", " ", "=", " ", 
      RowBox[{"J", ".", "w"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"p", " ", "=", " ", "0"}], ";", "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "0"}], ",", 
       RowBox[{"i", "\[LessEqual]", "n"}], ",", 
       RowBox[{"i", "=", 
        RowBox[{"i", "+", "1"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{"p", " ", "=", " ", 
        RowBox[{"p", " ", "+", " ", 
         RowBox[{
          RowBox[{"z", "[", 
           RowBox[{"[", 
            RowBox[{"i", "+", "1"}], "]"}], "]"}], " ", 
          RowBox[{"t", "^", "i"}]}]}]}]}], "]"}], ";", "\[IndentingNewLine]", 
     
     RowBox[{"q", " ", "=", " ", 
      RowBox[{
       RowBox[{"t", "^", 
        RowBox[{"(", 
         RowBox[{"n", "+", "1"}], ")"}]}], "-", "p"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"roots", " ", "=", " ", 
      RowBox[{"Solve", "[", 
       RowBox[{
        RowBox[{"q", "\[Equal]", "0"}], ",", "t"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "1"}], ",", 
       RowBox[{"i", "\[LessEqual]", 
        RowBox[{"n", "+", "1"}]}], ",", 
       RowBox[{"i", "=", 
        RowBox[{"i", "+", "1"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"x", "[", 
         RowBox[{"i", "-", "1"}], "]"}], " ", "=", " ", 
        RowBox[{
         RowBox[{
          RowBox[{"roots", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "[", 
          RowBox[{"[", "1", "]"}], "]"}], "[", 
         RowBox[{"[", "2", "]"}], "]"}]}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", " ", 
      RowBox[{
       RowBox[{"i", "=", "0"}], ",", 
       RowBox[{"i", "\[LessEqual]", "n"}], ",", 
       RowBox[{"i", "=", 
        RowBox[{"i", "+", "1"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"j", "=", "0"}], ",", 
         RowBox[{"j", "\[LessEqual]", "n"}], ",", 
         RowBox[{"j", "=", 
          RowBox[{"j", "+", "1"}]}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"TTab", "[", 
           RowBox[{"i", ",", "j"}], "]"}], "=", " ", 
          RowBox[{
           RowBox[{"x", "[", "j", "]"}], "^", "i"}]}]}], "]"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"V", "=", " ", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"TTab", "[", 
         RowBox[{"rr", ",", "ss"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"rr", ",", "0", ",", "n"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"ss", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "0"}], ",", " ", 
       RowBox[{"i", "\[LessEqual]", "n"}], ",", 
       RowBox[{"i", "=", 
        RowBox[{"i", "+", "1"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"vvv", "[", "i", "]"}], " ", "=", " ", 
        RowBox[{"bb", "[", "i", "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
     
     RowBox[{"w", " ", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"vvv", "[", "i", "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"weights", "=", 
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{"Inverse", "[", "V", "]"}], ".", "w"}], "]"}]}], ";"}], 
    ")"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"General", "::", "\<\"spell1\"\>"}], ":", 
  " ", "\<\"Possible spelling error: new symbol name \\\"\\!\\(roots\\)\\\" \
is similar to existing symbol \\\"\\!\\(Roots\\)\\\".\"\>"}]], "Message"],

Cell[BoxData[
 RowBox[{
  RowBox[{"General", "::", "\<\"spell1\"\>"}], ":", 
  " ", "\<\"Possible spelling error: new symbol name \\\"\\!\\(TTab\\)\\\" is \
similar to existing symbol \\\"\\!\\(Tab\\)\\\".\"\>"}]], "Message"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{
     "Iterated", " ", "Gaussian", " ", "quadrature", " ", "estimate", " ", 
      "of", " ", "\[IndentingNewLine]", "degree", " ", "2", "n"}], "+", 
     RowBox[{
     "1", " ", "for", " ", "integral", " ", "of", " ", "predefined", 
      "\[IndentingNewLine]", "function", " ", 
      RowBox[{"f", "[", 
       RowBox[{"x", ",", "y"}], "]"}], " ", "over", " ", 
      RowBox[{"rectangle", " ", "[", 
       RowBox[{"a", ",", "b"}], "]"}], " ", 
      RowBox[{"x", " ", "[", 
       RowBox[{"c", ",", "d"}], "]"}]}]}], ";", "\[IndentingNewLine]", "  ", 
    RowBox[{
     RowBox[{
     "exact", " ", "if", " ", "f", " ", "is", " ", "a", " ", "polynomial", 
      " ", "whose", " ", "degrees", " ", "in", " ", "x", " ", "and", " ", "y",
       "\[IndentingNewLine]", " ", "are", " ", "each"}], " ", "\[LessEqual]", 
     " ", 
     RowBox[{
      RowBox[{"2", "n"}], "+", "1"}]}], ";", " ", 
    RowBox[{
    "the", " ", "value", " ", "is", " ", "stored", " ", "in", " ", "II"}]}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"G2dim", "[", 
     RowBox[{"n_", ",", "a_", ",", "b_", ",", "c_", ",", "d_"}], "]"}], ":=", 
    
    RowBox[{"(", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Gaussian", "[", 
       RowBox[{"n", ",", "a", ",", "b"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"i", "=", "0"}], ",", 
        RowBox[{"i", "\[LessEqual]", "n"}], ",", 
        RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{
          RowBox[{"xx", "[", "i", "]"}], "=", 
          RowBox[{"x", "[", "i", "]"}]}], ";", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"xweight", "[", "i", "]"}], "=", " ", 
          RowBox[{"weights", "[", 
           RowBox[{"[", 
            RowBox[{"i", "+", "1"}], "]"}], "]"}]}]}]}], "]"}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Gaussian", "[", 
       RowBox[{"n", ",", "c", ",", "d"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"i", "=", "0"}], ",", 
        RowBox[{"i", "\[LessEqual]", "n"}], ",", 
        RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{
          RowBox[{"yy", "[", "i", "]"}], "=", 
          RowBox[{"x", "[", "i", "]"}]}], ";", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"yweight", "[", "i", "]"}], " ", "=", " ", 
          RowBox[{"weights", "[", 
           RowBox[{"[", 
            RowBox[{"i", "+", "1"}], "]"}], "]"}]}]}]}], "]"}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"II", "=", "0"}], ";", "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"j", "=", "0"}], ",", 
        RowBox[{"j", "\[LessEqual]", "n"}], ",", 
        RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
        RowBox[{"For", "[", 
         RowBox[{
          RowBox[{"i", "=", "0"}], ",", 
          RowBox[{"i", "\[LessEqual]", "n"}], ",", 
          RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
          RowBox[{"II", " ", "=", " ", 
           RowBox[{"II", " ", "+", " ", 
            RowBox[{
             RowBox[{"xweight", "[", "j", "]"}], " ", 
             RowBox[{"yweight", "[", "i", "]"}], " ", 
             RowBox[{"f", "[", 
              RowBox[{
               RowBox[{"xx", "[", "j", "]"}], ",", 
               RowBox[{"yy", "[", "i", "]"}]}], "]"}]}]}]}]}], "]"}]}], "]"}],
       ";"}], "\[IndentingNewLine]", ")"}]}], 
   "\[IndentingNewLine]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", 
    RowBox[{"x_", ",", "y_"}], "]"}], ":=", " ", 
   RowBox[{
    RowBox[{"x", "^", "7"}], " ", 
    RowBox[{"y", "^", "7"}]}]}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"a", "=", "1"}]], "Input"],

Cell[BoxData["1"], "Output"]
}, Open  ]],

Cell[CellGroupData[{

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

Cell[BoxData["2"], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"c", "=", "3"}]], "Input"],

Cell[BoxData["3"], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"d", "=", "5"}]], "Input"],

Cell[BoxData["5"], "Output"]
}, Open  ]],

Cell[CellGroupData[{

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

Cell[BoxData["3"], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NIntegrate", "[", 
  RowBox[{
   RowBox[{"f", "[", 
    RowBox[{"x", ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "a", ",", "b"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "c", ",", "d"}], "}"}]}], "]"}]], "Input"],

Cell[BoxData["1.5302550004199378`*^6"], "Output"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"G2dim", "[", 
  RowBox[{"n", ",", "a", ",", "b", ",", "c", ",", "d"}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["II"], "Input"],

Cell[BoxData["1.5302550145218275`*^6"], "Output"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
   "An", " ", "example", " ", "where", " ", "the", " ", "result", " ", "is", 
    " ", "exact", " ", "even", "\[IndentingNewLine]", "though", " ", "the", 
    " ", "degree", " ", "in", " ", "x", " ", "is", " ", "greater", " ", 
    "than", " ", "2", "n"}], "+", "1"}], "*)"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", 
    RowBox[{"x_", ",", "y_"}], "]"}], ":=", " ", 
   RowBox[{
    RowBox[{"x", "^", "9"}], " ", 
    RowBox[{"y", "^", "3"}]}]}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NIntegrate", "[", 
  RowBox[{
   RowBox[{"f", "[", 
    RowBox[{"x", ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "a", ",", "b"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "c", ",", "d"}], "}"}]}], "]"}]], "Input"],

Cell[BoxData["13912.799998833125`"], "Output"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"G2dim", "[", 
  RowBox[{"n", ",", "a", ",", "b", ",", "c", ",", "d"}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["II"], "Input"],

Cell[BoxData["13912.758373598572`"], "Output"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
     "An", " ", "example", " ", "where", " ", "degree", " ", "in", " ", "x"}],
      " ", "+", " ", 
     RowBox[{"degree", " ", "in", " ", "y"}]}], " ", "\[LessEqual]", " ", 
    RowBox[{"2", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], ",", "\[IndentingNewLine]", 
   " ", 
   RowBox[{
   "but", " ", "the", " ", "result", " ", "is", " ", "not", " ", "exact"}]}], 
  "*)"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"f", "[", 
   RowBox[{"x_", ",", "y_"}], "]"}], ":=", " ", 
  RowBox[{
   RowBox[{"x", "^", "10"}], " ", 
   RowBox[{"y", "^", "4"}]}]}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NIntegrate", "[", 
  RowBox[{
   RowBox[{"f", "[", 
    RowBox[{"x", ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "a", ",", "b"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "c", ",", "d"}], "}"}]}], "]"}]], "Input"],

Cell[BoxData["107262.80002135695`"], "Output"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"G2dim", "[", 
  RowBox[{"n", ",", "a", ",", "b", ",", "c", ",", "d"}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["II"], "Input"],

Cell[BoxData["107261.46947900583`"], "Output"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"A", " ", "non"}], "-", 
   RowBox[{"polynomial", " ", "example"}]}], "*)"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", 
    RowBox[{"x_", ",", "y_"}], "]"}], ":=", " ", 
   RowBox[{"Exp", "[", 
    RowBox[{
     RowBox[{"Sin", "[", "x", "]"}], " ", "+", "y"}], "]"}]}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NIntegrate", "[", 
  RowBox[{
   RowBox[{"f", "[", 
    RowBox[{"x", ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "a", ",", "b"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "c", ",", "d"}], "}"}]}], "]"}]], "Input"],

Cell[BoxData["334.2500266297184`"], "Output"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"G2dim", "[", 
  RowBox[{"n", ",", "a", ",", "b", ",", "c", ",", "d"}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["II"], "Input"],

Cell[BoxData["334.24991826644174`"], "Output"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
  "Example", " ", "where", " ", "a", " ", "difficult", " ", "integral", " ", 
   "leads", "\[IndentingNewLine]", "to", " ", "a", " ", 
   RowBox[{"(", "probably", ")"}], " ", "not", " ", "very", " ", "accurate", 
   " ", "estimate"}], "*)"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", 
    RowBox[{"x_", ",", "y_"}], "]"}], ":=", " ", 
   RowBox[{"Sin", "[", 
    RowBox[{"Exp", "[", 
     RowBox[{
      RowBox[{"x", "^", "2"}], " ", "+", " ", 
      RowBox[{"y", "^", "2"}]}], "]"}], "]"}]}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NIntegrate", "[", 
  RowBox[{
   RowBox[{"f", "[", 
    RowBox[{"x", ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "a", ",", "b"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "c", ",", "d"}], "}"}]}], "]"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"NIntegrate", "::", "\<\"slwcon\"\>"}], ":", 
  " ", "\<\"Numerical integration converging too slowly; suspect one of the \
following: singularity, value of the integration being 0, oscillatory \
integrand, or insufficient WorkingPrecision. If your integrand is oscillatory \
try using the option Method->Oscillatory in NIntegrate.\"\>"}]], "Message"],

Cell[BoxData["1.9808972282177625`*^-6"], "Output"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"G2dim", "[", 
  RowBox[{"n", ",", "a", ",", "b", ",", "c", ",", "d"}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["II"], "Input"],

Cell[BoxData[
 RowBox[{"-", "0.28675739436512604`"}]], "Output"]
}, Open  ]]
},
WindowSize->{1051, 893},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
DockedCells->FEPrivate`FrontEndResource[
 "FEExpressions", "CompatibilityToolbar"],
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
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[CellGroupData[{
Cell[590, 23, 7370, 202, 692, "Input"],
Cell[7963, 227, 228, 4, 22, "Message"],
Cell[8194, 233, 225, 4, 22, "Message"]
}, Open  ]],
Cell[8434, 240, 3602, 90, 412, "Input"],
Cell[12039, 332, 203, 7, 31, "Input"],
Cell[CellGroupData[{
Cell[12267, 343, 49, 1, 31, "Input"],
Cell[12319, 346, 28, 0, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12384, 351, 49, 1, 31, "Input"],
Cell[12436, 354, 28, 0, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12501, 359, 49, 1, 31, "Input"],
Cell[12553, 362, 28, 0, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12618, 367, 49, 1, 31, "Input"],
Cell[12670, 370, 28, 0, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12735, 375, 49, 1, 31, "Input"],
Cell[12787, 378, 28, 0, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12852, 383, 268, 8, 31, "Input"],
Cell[13123, 393, 49, 0, 30, "Output"]
}, Open  ]],
Cell[13187, 396, 111, 2, 31, "Input"],
Cell[CellGroupData[{
Cell[13323, 402, 28, 0, 31, "Input"],
Cell[13354, 404, 49, 0, 30, "Output"]
}, Open  ]],
Cell[13418, 407, 339, 7, 52, "Input"],
Cell[13760, 416, 203, 7, 31, "Input"],
Cell[CellGroupData[{
Cell[13988, 427, 268, 8, 31, "Input"],
Cell[14259, 437, 46, 0, 30, "Output"]
}, Open  ]],
Cell[14320, 440, 111, 2, 31, "Input"],
Cell[CellGroupData[{
Cell[14456, 446, 28, 0, 31, "Input"],
Cell[14487, 448, 46, 0, 30, "Output"]
}, Open  ]],
Cell[14548, 451, 506, 16, 52, "Input"],
Cell[15057, 469, 181, 6, 31, "Input"],
Cell[CellGroupData[{
Cell[15263, 479, 268, 8, 31, "Input"],
Cell[15534, 489, 46, 0, 30, "Output"]
}, Open  ]],
Cell[15595, 492, 111, 2, 31, "Input"],
Cell[CellGroupData[{
Cell[15731, 498, 28, 0, 31, "Input"],
Cell[15762, 500, 46, 0, 30, "Output"]
}, Open  ]],
Cell[15823, 503, 139, 4, 31, "Input"],
Cell[15965, 509, 223, 7, 31, "Input"],
Cell[CellGroupData[{
Cell[16213, 520, 268, 8, 31, "Input"],
Cell[16484, 530, 45, 0, 30, "Output"]
}, Open  ]],
Cell[16544, 533, 111, 2, 31, "Input"],
Cell[CellGroupData[{
Cell[16680, 539, 28, 0, 31, "Input"],
Cell[16711, 541, 46, 0, 30, "Output"]
}, Open  ]],
Cell[16772, 544, 293, 6, 52, "Input"],
Cell[17068, 552, 282, 9, 31, "Input"],
Cell[CellGroupData[{
Cell[17375, 565, 268, 8, 31, "Input"],
Cell[17646, 575, 385, 6, 77, "Message"],
Cell[18034, 583, 50, 0, 30, "Output"]
}, Open  ]],
Cell[18099, 586, 111, 2, 31, "Input"],
Cell[CellGroupData[{
Cell[18235, 592, 28, 0, 31, "Input"],
Cell[18266, 594, 64, 1, 30, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)