(* 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[ 253012, 6372] NotebookOptionsPosition[ 229396, 5761] NotebookOutlinePosition[ 229791, 5778] CellTagsIndexPosition[ 229748, 5775] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"PyraGaussRuleInfo", "[", " ", RowBox[{ RowBox[{"{", RowBox[{"rule_", ",", "numer_"}], "}"}], ",", "point_"}], "]"}], " ", ":=", " ", RowBox[{"Module", " ", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{ "g1", ",", "g2", ",", "g3", ",", "g4", ",", "g5", ",", "w1", ",", "w2", ",", "w3", ",", "w4", ",", RowBox[{"jk", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}], ",", "\[IndentingNewLine]", RowBox[{"jk4", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", RowBox[{"jk9", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "1"}], "}"}]}], "}"}]}], ",", "\[IndentingNewLine]", RowBox[{"wg9", "=", RowBox[{ RowBox[{"{", RowBox[{"64", ",", "40", ",", "25"}], "}"}], "/", "81"}]}], ",", "\[IndentingNewLine]", RowBox[{"eps", "=", RowBox[{"10", "^", "16"}]}], ",", "\[IndentingNewLine]", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"Null", ",", "Null", ",", "Null"}], "}"}], ",", "0"}], "}"}]}], ",", "\[IndentingNewLine]", "j", ",", "k", ",", "m", ",", "\[IndentingNewLine]", RowBox[{"p", "=", "rule"}], ",", "\[IndentingNewLine]", RowBox[{"i", "=", "point"}]}], "}"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"If", " ", "[", " ", RowBox[{ RowBox[{"p", " ", "\[Equal]", " ", "1"}], ",", " ", "\[IndentingNewLine]", RowBox[{"info", " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{ RowBox[{"-", "1"}], "/", "2"}]}], "}"}], ",", RowBox[{"128", "/", "27"}]}], "}"}]}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"If", " ", "[", " ", RowBox[{ RowBox[{"p", " ", "==", " ", "5"}], ",", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"g1", " ", "=", " ", RowBox[{"8", "*", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{"2", "/", "15"}], "]"}], "/", "5"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", " ", RowBox[{ RowBox[{"i", "<", "5"}], ",", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk", "[", RowBox[{"[", "i", "]"}], "]"}]}], ";", " ", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g1"}], ",", RowBox[{"k", "*", "g1"}], ",", RowBox[{ RowBox[{"-", "2"}], "/", "3"}]}], "}"}], ",", RowBox[{"81", "/", "100"}]}], "}"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[Equal]", "5"}], ",", " ", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{"2", "/", "5"}]}], "}"}], ",", RowBox[{"125", "/", "27"}]}], "}"}]}]}], "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"If", " ", "[", " ", RowBox[{ RowBox[{"p", "\[Equal]", "6"}], ",", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"g1", "=", RowBox[{"Sqrt", "[", RowBox[{"12", "/", "35"}], "]"}]}], ";", "\[IndentingNewLine]", " ", RowBox[{"g2", "=", RowBox[{"{", RowBox[{ RowBox[{"1", "/", "6"}], ",", RowBox[{"1", "/", "2"}]}], "}"}]}], ";", " ", "\[IndentingNewLine]", RowBox[{"w2", "=", RowBox[{"{", RowBox[{ RowBox[{"576", "/", "625"}], ",", RowBox[{"64", "/", "15"}]}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "<", "5"}], ",", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk", "[", RowBox[{"[", "i", "]"}], "]"}]}], ";", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g1"}], ",", RowBox[{"k", "*", "g1"}], ",", RowBox[{ RowBox[{"-", "2"}], "/", "3"}]}], "}"}], ",", RowBox[{"504", "/", "625"}]}], "}"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", ">", "4"}], ",", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{"g2", "[", RowBox[{"[", RowBox[{"i", "-", "4"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"w2", "[", RowBox[{"[", RowBox[{"i", "-", "4"}], "]"}], "]"}]}], "}"}]}]}], "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"If", " ", "[", " ", RowBox[{ RowBox[{"p", " ", "\[Equal]", " ", "8"}], ",", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"g1", "=", RowBox[{"Sqrt", "[", RowBox[{"1", "/", "3"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"g2", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"2", "*", RowBox[{"Sqrt", "[", "10", "]"}]}], "-", "5"}], ")"}], "/", "15"}]}], ";", "\[IndentingNewLine]", RowBox[{"g3", "=", RowBox[{ RowBox[{ RowBox[{"-", "2"}], "/", "3"}], "-", "g2"}]}], ";", "\[IndentingNewLine]", RowBox[{"w1", "=", RowBox[{"5", "*", RowBox[{ RowBox[{"(", RowBox[{"68", "+", RowBox[{"5", "*", RowBox[{"Sqrt", "[", "10", "]"}]}]}], ")"}], "/", "432"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"w2", "=", RowBox[{ RowBox[{"85", "/", "54"}], "-", "w1"}]}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "<", "5"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk", "[", RowBox[{"[", "i", "]"}], "]"}]}], ";", " ", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g1"}], ",", RowBox[{"k", "*", "g1"}], ",", "g2"}], "}"}], ",", "w1"}], "}"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", ">", "4"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk", "[", RowBox[{"[", RowBox[{"i", "-", "4"}], "]"}], "]"}]}], ";", " ", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g1"}], ",", RowBox[{"k", "*", "g1"}], ",", "g3"}], "}"}], ",", "w2"}], "}"}]}]}]}], "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"If", "[", " ", RowBox[{ RowBox[{"p", " ", "\[Equal]", " ", RowBox[{"-", "8"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"g1", " ", "=", " ", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"(", RowBox[{"2", "/", "15"}], ")"}], "*", RowBox[{"(", RowBox[{"573", "-", RowBox[{"2", "*", RowBox[{"Sqrt", "[", "51", "]"}]}]}], ")"}]}], "]"}], "/", "15"}]}], ";", "\[IndentingNewLine]", RowBox[{"g2", " ", "=", " ", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"(", RowBox[{"2", "/", "15"}], ")"}], "*", RowBox[{"(", RowBox[{"573", "+", RowBox[{"2", "*", RowBox[{"Sqrt", "[", "51", "]"}]}]}], ")"}]}], "]"}], "/", "15"}]}], ";", "\[IndentingNewLine]", RowBox[{"g3", "=", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ RowBox[{"2", "*", RowBox[{"Sqrt", "[", "51", "]"}]}], "+", "13"}], ")"}]}], "/", "35"}]}], ";", "\[IndentingNewLine]", RowBox[{"g4", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"2", "*", RowBox[{"Sqrt", "[", "51", "]"}]}], "-", "13"}], ")"}], "/", "35"}]}], ";", "\[IndentingNewLine]", RowBox[{"w1", "=", RowBox[{ RowBox[{"(", RowBox[{"11764", "-", RowBox[{"461", "*", RowBox[{"Sqrt", "[", "51", "]"}]}]}], ")"}], "/", "15300"}]}], ";", "\[IndentingNewLine]", RowBox[{"w2", " ", "=", " ", RowBox[{ RowBox[{"346", "/", "225"}], "-", "w1"}]}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "<", "5"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk", "[", RowBox[{"[", "i", "]"}], "]"}]}], ";", " ", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g1"}], ",", RowBox[{"k", "*", "g1"}], ",", "g3"}], "}"}], ",", "w1"}], "}"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", ">", "4"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk", "[", RowBox[{"[", RowBox[{"i", "-", "4"}], "]"}], "]"}]}], ";", " ", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g2"}], ",", RowBox[{"k", "*", "g2"}], ",", "g4"}], "}"}], ",", "w2"}], "}"}]}]}]}], "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"p", "\[Equal]", "9"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"g1", "=", RowBox[{"8", "*", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"(", RowBox[{"573", "+", RowBox[{"5", "*", RowBox[{"Sqrt", "[", "2865", "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{"109825", "+", RowBox[{"969", "*", RowBox[{"Sqrt", "[", "2865", "]"}]}]}], ")"}]}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"g2", "=", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"(", RowBox[{"2", "*", RowBox[{"(", RowBox[{"8025", "+", RowBox[{"Sqrt", "[", "2865", "]"}]}], ")"}]}], ")"}], "/", "35"}], "]"}], "/", "37"}]}], ";", "\[IndentingNewLine]", RowBox[{"g3", "=", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"87", "+", RowBox[{"Sqrt", "[", "2865", "]"}]}], ")"}]}], "/", "168"}]}], ";", "\[IndentingNewLine]", RowBox[{"g4", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "87"}], "+", RowBox[{"Sqrt", "[", "2865", "]"}]}], ")"}], "/", "168"}]}], ";", "\[IndentingNewLine]", RowBox[{"w1", "=", RowBox[{"7", "*", RowBox[{ RowBox[{"(", RowBox[{"11472415", "-", RowBox[{"70057", "*", RowBox[{"Sqrt", "[", "2865", "]"}]}]}], ")"}], "/", "130739500"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"w2", "=", RowBox[{ RowBox[{"84091", "/", "68450"}], "-", "w1"}]}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "<", "5"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk", "[", RowBox[{"[", "i", "]"}], "]"}]}], ";", " ", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g1"}], ",", RowBox[{"k", "*", "g1"}], ",", "g3"}], "}"}], ",", "w1"}], "}"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"i", ">", "4"}], "&&", RowBox[{"i", "<", "9"}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk", "[", RowBox[{"[", RowBox[{"i", "-", "4"}], "]"}], "]"}]}], ";", " ", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g2"}], ",", RowBox[{"k", "*", "g2"}], ",", "g4"}], "}"}], ",", "w2"}], "}"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[Equal]", "9"}], ",", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{"2", "/", "3"}]}], "}"}], ",", RowBox[{"18", "/", "5"}]}], "}"}]}]}], "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"If", "[", " ", RowBox[{ RowBox[{"p", " ", "\[Equal]", " ", "13"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"g1", "=", RowBox[{"7", "*", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{"35", "/", "59"}], "]"}], "/", "8"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"g2", "=", RowBox[{"224", "*", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{"336633710", "/", "33088740423"}], "]"}], "/", "37"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"g3", "=", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{"37043", "/", "35"}], "]"}], "/", "56"}]}], ";", "\[IndentingNewLine]", RowBox[{"g4", "=", RowBox[{ RowBox[{"-", "127"}], "/", "153"}]}], ";", "\[IndentingNewLine]", RowBox[{"g5", "=", RowBox[{"1490761", "/", "2842826"}]}], ";", "\[IndentingNewLine]", RowBox[{"w1", "=", RowBox[{"170569", "/", "331200"}]}], ";", "\[IndentingNewLine]", RowBox[{"w2", "=", RowBox[{"276710106577408", "/", "1075923777052725"}]}], ";", "\[IndentingNewLine]", RowBox[{"w3", "=", RowBox[{"12827693806929", "/", "30577384040000"}]}], ";", "\[IndentingNewLine]", RowBox[{"w4", "=", RowBox[{ "10663383340655070643544192", "/", "4310170528879365193704375"}]}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "<", "5"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk", "[", RowBox[{"[", "i", "]"}], "]"}]}], ";", " ", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g1"}], ",", RowBox[{"k", "*", "g1"}], ",", RowBox[{ RowBox[{"-", "1"}], "/", "7"}]}], "}"}], ",", "w1"}], "}"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"i", ">", "4"}], "&&", RowBox[{"i", "<", "9"}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk4", "[", RowBox[{"[", RowBox[{"i", "-", "4"}], "]"}], "]"}]}], ";", " ", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g2"}], ",", RowBox[{"k", "*", "g2"}], ",", RowBox[{ RowBox[{"-", "9"}], "/", "28"}]}], "}"}], ",", "w2"}], "}"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"i", ">", "8"}], "&&", RowBox[{"i", "<", "13"}]}], ",", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk", "[", RowBox[{"[", RowBox[{"i", "-", "8"}], "]"}], "]"}]}], ";", " ", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g3"}], ",", RowBox[{"k", "*", "g3"}], ",", "g4"}], "}"}], ",", "w3"}], "}"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[Equal]", "13"}], ",", " ", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "g5"}], "}"}], ",", "w4"}], "}"}]}]}], "]"}]}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"p", "\[Equal]", "18"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"g1", "=", RowBox[{"Sqrt", "[", RowBox[{"3", "/", "5"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"g2", "=", RowBox[{"1", "-", RowBox[{"2", "*", RowBox[{ RowBox[{"(", RowBox[{"10", "-", RowBox[{"Sqrt", "[", "10", "]"}]}], ")"}], "/", "15"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"g3", "=", RowBox[{ RowBox[{ RowBox[{"-", "2"}], "/", "3"}], "-", "g2"}]}], ";", "\[IndentingNewLine]", RowBox[{"w1", "=", RowBox[{"5", "*", RowBox[{ RowBox[{"(", RowBox[{"68", "+", RowBox[{"5", "*", RowBox[{"Sqrt", "[", "10", "]"}]}]}], ")"}], "/", "432"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"w2", "=", RowBox[{ RowBox[{"85", "/", "54"}], "-", "w1"}]}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "<", "10"}], ",", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk9", "[", RowBox[{"[", "i", "]"}], "]"}]}], ";", " ", RowBox[{"m", "=", RowBox[{ RowBox[{"Abs", "[", "j", "]"}], "+", RowBox[{"Abs", "[", "k", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g1"}], ",", RowBox[{"k", "*", "g1"}], ",", "g2"}], "}"}], ",", RowBox[{"w1", "*", RowBox[{"wg9", "[", RowBox[{"[", RowBox[{"m", "+", "1"}], "]"}], "]"}]}]}], "}"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"i", ">", "9"}], "&&", RowBox[{"i", "<", "19"}]}], ",", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk9", "[", RowBox[{"[", RowBox[{"i", "-", "9"}], "]"}], "]"}]}], ";", RowBox[{"m", "=", RowBox[{ RowBox[{"Abs", "[", "j", "]"}], "+", RowBox[{"Abs", "[", "k", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g1"}], ",", RowBox[{"k", "*", "g1"}], ",", "g3"}], "}"}], ",", RowBox[{"w2", "*", RowBox[{"wg9", "[", RowBox[{"[", RowBox[{"m", "+", "1"}], "]"}], "]"}]}]}], "}"}]}]}]}], "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"p", "\[Equal]", "27"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"g1", "=", RowBox[{"Sqrt", "[", RowBox[{"3", "/", "5"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"g3", "=", RowBox[{"-", "0.854011951853700535688324041975993416"}]}], ";", "\[IndentingNewLine]", RowBox[{"g4", "=", RowBox[{"-", "0.305992467923296230556472913192103090"}]}], ";", "\[IndentingNewLine]", RowBox[{"g5", "=", RowBox[{"+", "0.410004419776996766244796955168096505"}]}], ";", "\[IndentingNewLine]", RowBox[{"If", " ", "[", RowBox[{ RowBox[{"!", "numer"}], ",", " ", RowBox[{ RowBox[{"{", RowBox[{"g3", ",", "g4", ",", "g5"}], "}"}], "=", RowBox[{"Rationalize", "[", RowBox[{ RowBox[{"{", RowBox[{"g3", ",", "g4", ",", "g5"}], "}"}], ",", "eps"}], "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"w1", "=", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "15"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"4", "+", RowBox[{"5", "*", RowBox[{"(", RowBox[{"g4", "+", "g5"}], ")"}]}], "+", RowBox[{"10", "*", "g4", "*", "g5"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"g3", "-", "g4"}], ")"}], "*", RowBox[{"(", RowBox[{"g3", "-", "g5"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "g3"}], ")"}], "^", "2"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"w2", "=", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "15"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"4", "+", RowBox[{"5", "*", RowBox[{"(", RowBox[{"g3", "+", "g5"}], ")"}]}], "+", RowBox[{"10", "*", "g3", "*", "g5"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"g3", "-", "g4"}], ")"}], "*", RowBox[{"(", RowBox[{"g5", "-", "g4"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "g4"}], ")"}], "^", "2"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"w3", "=", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "15"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"4", "+", RowBox[{"5", "*", RowBox[{"(", RowBox[{"g3", "+", "g4"}], ")"}]}], "+", RowBox[{"10", "*", "g3", "*", "g4"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"g3", "-", "g5"}], ")"}], "*", RowBox[{"(", RowBox[{"g4", "-", "g5"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "g5"}], ")"}], "^", "2"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "<", "10"}], ",", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk9", "[", RowBox[{"[", "i", "]"}], "]"}]}], ";", RowBox[{"m", "=", RowBox[{ RowBox[{"Abs", "[", "j", "]"}], "+", RowBox[{"Abs", "[", "k", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g1"}], ",", RowBox[{"k", "*", "g1"}], ",", "g3"}], "}"}], ",", RowBox[{"w1", "*", RowBox[{"wg9", "[", RowBox[{"[", RowBox[{"m", "+", "1"}], "]"}], "]"}]}]}], "}"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"i", ">", "9"}], "&&", RowBox[{"i", "<", "19"}]}], ",", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk9", "[", RowBox[{"[", RowBox[{"i", "-", "9"}], "]"}], "]"}]}], ";", RowBox[{"m", "=", RowBox[{ RowBox[{"Abs", "[", "j", "]"}], "+", RowBox[{"Abs", "[", "k", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g1"}], ",", RowBox[{"k", "*", "g1"}], ",", "g4"}], "}"}], ",", RowBox[{"w2", "*", RowBox[{"wg9", "[", RowBox[{"[", RowBox[{"m", "+", "1"}], "]"}], "]"}]}]}], "}"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"i", ">", "18"}], ",", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}], "=", RowBox[{"jk9", "[", RowBox[{"[", RowBox[{"i", "-", "18"}], "]"}], "]"}]}], ";", RowBox[{"m", "=", RowBox[{ RowBox[{"Abs", "[", "j", "]"}], "+", RowBox[{"Abs", "[", "k", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"info", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"j", "*", "g1"}], ",", RowBox[{"k", "*", "g1"}], ",", "g5"}], "}"}], ",", RowBox[{"w3", "*", RowBox[{"wg9", "[", RowBox[{"[", RowBox[{"m", "+", "1"}], "]"}], "]"}]}]}], "}"}]}]}]}], "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"If", " ", "[", " ", RowBox[{"numer", ",", " ", RowBox[{"Return", "[", RowBox[{"N", "[", RowBox[{"info", ",", "20"}], "]"}], "]"}], ",", " ", RowBox[{"Return", "[", "info", "]"}]}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]"}]], "Input",\ CellChangeTimes->{{3.413376886297916*^9, 3.4133773256082172`*^9}, { 3.413378742795382*^9, 3.413378888399621*^9}, {3.413378923615098*^9, 3.4133789328638353`*^9}, {3.413396985311747*^9, 3.413397147592699*^9}, { 3.413653443541477*^9, 3.413653597200617*^9}, {3.41365363703336*^9, 3.413653900925736*^9}, {3.413654239162406*^9, 3.413654567881105*^9}, { 3.413654917203484*^9, 3.413654925192285*^9}, {3.4136549840012903`*^9, 3.413654985096965*^9}, {3.413655112059662*^9, 3.413655115833828*^9}, { 3.413655182544448*^9, 3.413655403531612*^9}, {3.413655645788496*^9, 3.413655858704184*^9}, {3.4136558968602753`*^9, 3.4136559696428537`*^9}, { 3.413657095602586*^9, 3.4136573081847973`*^9}, 3.413819685894218*^9, { 3.4138197550324383`*^9, 3.413819756256186*^9}, 3.413819882135359*^9, { 3.414673591270568*^9, 3.414673592847736*^9}, {3.414746509322014*^9, 3.414746510392467*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< Mathematica version\>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ "Print", "[", "\"\< Test the quadrature rules for a pyramid, returned\>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ "Print", "[", "\"\< by {{x,y,z},w} = PyraGaussRuleInfo[{rule,numer},i]\>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< INPUT:\>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ "Print", "[", "\"\< RULE is 1, 5, 6, 8, -8, 9, 13, 18 or 27,\>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ "Print", "[", "\"\< the order of the rule. There are two rules of order 8,\>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ "Print", "[", "\"\< so the value of -8 returns the second one.\>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ "Print", "[", "\"\< NUMER is True to return numerical values for the output, or\>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ "Print", "[", "\"\< False to return symbolic (exact, except for rule 27) values.\>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"Print", "[", "\"\< I is the index, between 1 and |RULE|, of the point desired.\>\"", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< OUTPUT:\>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"Print", "[", "\"\< {X,Y,Z} are the coordinates of the quadrature point,\>\"", "]"}], "\[IndentingNewLine]", RowBox[{"Print", "[", "\"\< with -1 <= X, Y, Z <= +1. X and Y represent the base\>\"", "]"}], "\[IndentingNewLine]", RowBox[{"Print", "[", "\"\< of the pyramid, which Z measures the height.\>\"", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"Print", "[", "\"\< W is the weight of the quadrature point, and should be\>\"", "]"}], "\[IndentingNewLine]", RowBox[{"Print", "[", "\"\< multiplied by (1/4) (1-Z)^2 when used in a sum.\>\"", "]"}], "\[IndentingNewLine]", RowBox[{"Print", "[", "\"\< \>\"", "]"}], "\[IndentingNewLine]", RowBox[{"Print", "[", "\"\< The following cells print out each rule, and the sum\>\"", "]"}], "\[IndentingNewLine]", RowBox[{"Print", "[", "\"\< of the weights, scaled by (1/4) (1-Z)^2. The scaled\>\"", "]"}], "\[IndentingNewLine]", RowBox[{"Print", "[", "\"\< weight should sum to 8/3 or 2.666..., the volume of\>\"", "]"}], "\[IndentingNewLine]", RowBox[{"Print", "[", "\"\< the pyramid.\>\"", "]"}]}], "Input", CellChangeTimes->{{3.413820383448031*^9, 3.413820964218034*^9}, { 3.4138210017963247`*^9, 3.413821002539048*^9}}], Cell[CellGroupData[{ Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413820967703148*^9}], Cell[BoxData["\<\"PyraGaussRuleInfoTest.nb\"\>"], "Print", CellChangeTimes->{3.413820967765345*^9}], Cell[BoxData["\<\" Mathematica version\"\>"], "Print", CellChangeTimes->{3.413820967811675*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413820967862602*^9}], Cell[BoxData["\<\" Test the quadrature rules for a pyramid, returned\"\>"], \ "Print", CellChangeTimes->{3.413820967902644*^9}], Cell[BoxData["\<\" by {{x,y,z},w} = PyraGaussRuleInfo[{rule,numer},i]\"\>"], \ "Print", CellChangeTimes->{3.413820967958683*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413820967986137*^9}], Cell[BoxData["\<\" INPUT:\"\>"], "Print", CellChangeTimes->{3.4138209680176373`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413820968050599*^9}], Cell[BoxData["\<\" RULE is 1, 5, 6, 8, -8, 9, 13, 18 or 27,\"\>"], "Print", CellChangeTimes->{3.413820968083971*^9}], Cell[BoxData["\<\" the order of the rule. There are two rules of order 8,\"\ \>"], "Print", CellChangeTimes->{3.413820968117506*^9}], Cell[BoxData["\<\" so the value of -8 returns the second one.\"\>"], "Print", CellChangeTimes->{3.4138209681669493`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413820968193915*^9}], Cell[BoxData["\<\" NUMER is True to return numerical values for the output, \ or\"\>"], "Print", CellChangeTimes->{3.4138209682328*^9}], Cell[BoxData["\<\" False to return symbolic (exact, except for rule 27) \ values.\"\>"], "Print", CellChangeTimes->{3.4138209682678947`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413820968301022*^9}], Cell[BoxData["\<\" I is the index, between 1 and |RULE|, of the point \ desired.\"\>"], "Print", CellChangeTimes->{3.413820968337451*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413820968368371*^9}], Cell[BoxData["\<\" OUTPUT:\"\>"], "Print", CellChangeTimes->{3.413820968401188*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.4138209684343357`*^9}], Cell[BoxData["\<\" {X,Y,Z} are the coordinates of the quadrature \ point,\"\>"], "Print", CellChangeTimes->{3.413820968469121*^9}], Cell[BoxData["\<\" with -1 <= X, Y, Z <= +1. X and Y represent the \ base\"\>"], "Print", CellChangeTimes->{3.413820968506566*^9}], Cell[BoxData["\<\" of the pyramid, which Z measures the height.\"\>"], \ "Print", CellChangeTimes->{3.41382096853677*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413820968568121*^9}], Cell[BoxData["\<\" W is the weight of the quadrature point, and should \ be\"\>"], "Print", CellChangeTimes->{3.413820968601959*^9}], Cell[BoxData["\<\" multiplied by (1/4) (1-Z)^2 when used in a sum.\"\>"], \ "Print", CellChangeTimes->{3.413820968652543*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413820968682118*^9}], Cell[BoxData["\<\" The following cells print out each rule, and the \ sum\"\>"], "Print", CellChangeTimes->{3.4138209687183237`*^9}], Cell[BoxData["\<\" of the weights (scaled by (1/4) (1-Z)^2. The \ scaled\"\>"], "Print", CellChangeTimes->{3.4138209687513113`*^9}], Cell[BoxData["\<\" weight should sum to 8/3 or 2.666..., the volume \ of\"\>"], "Print", CellChangeTimes->{3.413820968784669*^9}], Cell[BoxData["\<\" the pyramid.\"\>"], "Print", CellChangeTimes->{3.413820968818873*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"rule", " ", "=", "1"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"numer", " ", "=", " ", "True"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"wsum", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "rule"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"w", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"3", "/", "8"}], " ", ")"}], " ", "*", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], ")"}], "^", "2"}], " ", "*", " ", "w"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "w", " ", "]"}], ";", "\[IndentingNewLine]", RowBox[{"wsum", " ", "=", " ", RowBox[{"wsum", " ", "+", "w"}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "wsum"}], " ", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"x", " ", "=", " ", RowBox[{"a", " ", "*", " ", RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], " ", ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "x", " ", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", " ", "=", " ", RowBox[{"b", " ", "*", " ", RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], " ", ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "y", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "z", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.413654633924918*^9, 3.413654799625798*^9}, { 3.41365484902993*^9, 3.413654878240939*^9}, {3.413654957362526*^9, 3.413654963399206*^9}, {3.413655009030931*^9, 3.413655015807214*^9}, { 3.413655056153854*^9, 3.4136550564604473`*^9}, {3.41365560416474*^9, 3.4136556090032673`*^9}, {3.4138136887636223`*^9, 3.413813699871728*^9}, { 3.414673677845633*^9, 3.4146737578261414`*^9}, {3.4146737973797894`*^9, 3.414673820186462*^9}, {3.4152015668451138`*^9, 3.4152015975359373`*^9}}], Cell[CellGroupData[{ Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598441349*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Pyramid rule number \"\>", "\[InvisibleSpace]", "1"}], SequenceForm["Pyramid rule number ", 1], Editable->False]], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.4152015984790087`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598519987*^9}], Cell[BoxData["\<\"W\"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.4152015985537853`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598587816*^9}], Cell[BoxData["0.99999999999999999999999999999684455638`19.778151250383644"], \ "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598621944*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"sum(w) = \"\>", "\[InvisibleSpace]", "0.99999999999999999999999999999684455638`19.778151250383644"}], SequenceForm["sum(w) = ", 1.`19.778151250383644], Editable->False]], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598656389*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598680682*^9}], Cell[BoxData["\<\"X\"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598710535*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.4152015987449417`*^9}], Cell[BoxData["0"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598778454*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598811008*^9}], Cell[BoxData["\<\"Y\"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598844962*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598900466*^9}], Cell[BoxData["0"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.4152015989253387`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598949967*^9}], Cell[BoxData["\<\"Z\"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201598978764*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201599011983*^9}], Cell[BoxData["0.25`20."], "Print", CellChangeTimes->{3.4136550570040913`*^9, 3.413655609904941*^9, 3.413813706748934*^9, 3.41467376248269*^9, 3.414673828767861*^9, 3.415201599045574*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"rule", " ", "=", " ", "5"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"numer", " ", "=", " ", "True"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"wsum", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "rule"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"w", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"3", "/", "8"}], " ", ")"}], " ", "*", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], ")"}], "^", "2"}], " ", "*", " ", "w"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "w", " ", "]"}], ";", "\[IndentingNewLine]", RowBox[{"wsum", " ", "=", " ", RowBox[{"wsum", " ", "+", "w"}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "wsum"}], " ", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"x", " ", "=", " ", RowBox[{"a", " ", "*", " ", RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], " ", ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "x", " ", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", " ", "=", " ", RowBox[{"b", " ", "*", " ", RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], " ", ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "y", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "z", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.413654633924918*^9, 3.413654799625798*^9}, { 3.41365484902993*^9, 3.413654878240939*^9}, {3.413654957362526*^9, 3.413654963399206*^9}, {3.413655009030931*^9, 3.413655015807214*^9}, { 3.4136555956892033`*^9, 3.413655599497733*^9}, {3.4138137133463373`*^9, 3.413813722344843*^9}, {3.414673847465712*^9, 3.414673847940092*^9}, { 3.415201081249057*^9, 3.4152013249591084`*^9}, {3.415201392820709*^9, 3.415201396805786*^9}}], Cell[CellGroupData[{ Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201398584159*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Pyramid rule number \"\>", "\[InvisibleSpace]", "5"}], SequenceForm["Pyramid rule number ", 5], Editable->False]], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.4152013986182833`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201398648987*^9}], Cell[BoxData["\<\"W\"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201398682036*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201398715253*^9}], Cell[BoxData["0.2109375`19.744727494896697"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.4152013987487383`*^9}], Cell[BoxData["0.2109375`19.744727494896697"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.4152013987828903`*^9}], Cell[BoxData["0.2109375`19.744727494896697"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201398813858*^9}], Cell[BoxData["0.2109375`19.744727494896697"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.41520139884013*^9}], Cell[BoxData["0.15624999999999999999999999999987378226`19.63202321470541"], \ "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201398874689*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"sum(w) = \"\>", "\[InvisibleSpace]", "0.99999999999999999999999999999987378226`19.725072806900226"}], SequenceForm["sum(w) = ", 1.`19.725072806900226], Editable->False]], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.4152013989079247`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201398941025*^9}], Cell[BoxData["\<\"X\"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201398975898*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.4152013990086813`*^9}], Cell[BoxData[ RowBox[{"-", "0.48686449556014765640619536248960189685`19.85387196432177"}]], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.41520139906431*^9}], Cell[BoxData["0.48686449556014765640619536248960189685`19.85387196432177"], \ "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399087089*^9}], Cell[BoxData["0.48686449556014765640619536248960189685`19.85387196432177"], \ "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399109806*^9}], Cell[BoxData[ RowBox[{"-", "0.48686449556014765640619536248960189685`19.85387196432177"}]], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399142003*^9}], Cell[BoxData["0"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.4152013991754723`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399209103*^9}], Cell[BoxData["\<\"Y\"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399242753*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399276143*^9}], Cell[BoxData[ RowBox[{"-", "0.48686449556014765640619536248960189685`19.85387196432177"}]], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399309084*^9}], Cell[BoxData[ RowBox[{"-", "0.48686449556014765640619536248960189685`19.85387196432177"}]], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399342329*^9}], Cell[BoxData["0.48686449556014765640619536248960189685`19.85387196432177"], \ "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399375963*^9}], Cell[BoxData["0.48686449556014765640619536248960189685`19.85387196432177"], \ "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.4152013994095697`*^9}], Cell[BoxData["0"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.41520139944326*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.4152013994761553`*^9}], Cell[BoxData["\<\"Z\"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399529716*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.4152013995541983`*^9}], Cell[BoxData["0.16666666666666666666666666666666666666`19.698970004336026"], \ "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399578861*^9}], Cell[BoxData["0.16666666666666666666666666666666666666`19.698970004336026"], \ "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399610002*^9}], Cell[BoxData["0.16666666666666666666666666666666666666`19.698970004336026"], \ "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399643691*^9}], Cell[BoxData["0.16666666666666666666666666666666666666`19.698970004336026"], \ "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399677368*^9}], Cell[BoxData["0.7`20.54406804435028"], "Print", CellChangeTimes->{ 3.413654888980755*^9, {3.4136549957687607`*^9, 3.4136550179711027`*^9}, 3.4136556003409863`*^9, 3.413813725817009*^9, 3.4146738522661257`*^9, { 3.415201328468487*^9, 3.415201341158718*^9}, 3.415201399709393*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"rule", " ", "=", " ", "6"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"numer", " ", "=", " ", "True"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"wsum", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "rule"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"w", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"3", "/", "8"}], " ", ")"}], " ", "*", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], ")"}], "^", "2"}], " ", "*", " ", "w"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "w", " ", "]"}], ";", "\[IndentingNewLine]", RowBox[{"wsum", " ", "=", " ", RowBox[{"wsum", " ", "+", "w"}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "wsum"}], " ", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"x", " ", "=", " ", RowBox[{"a", " ", "*", " ", RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], " ", ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "x", " ", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", " ", "=", " ", RowBox[{"b", " ", "*", " ", RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], " ", ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "y", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "z", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}]}], "Input", CellChangeTimes->{{3.413654633924918*^9, 3.413654799625798*^9}, { 3.41365484902993*^9, 3.413654878240939*^9}, {3.413654957362526*^9, 3.413654963399206*^9}, {3.413655009030931*^9, 3.413655015807214*^9}, { 3.413655074777342*^9, 3.413655075108171*^9}, {3.4136555839228897`*^9, 3.413655587706215*^9}, {3.4138136414533052`*^9, 3.413813652649139*^9}, { 3.413813735131393*^9, 3.413813741542747*^9}, {3.4146738704055853`*^9, 3.4146738708597116`*^9}, {3.415201436066044*^9, 3.4152014387062798`*^9}}], Cell[CellGroupData[{ Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439249247*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Pyramid rule number \"\>", "\[InvisibleSpace]", "6"}], SequenceForm["Pyramid rule number ", 6], Editable->False]], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439280696*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439312674*^9}], Cell[BoxData["\<\"W\"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.4152014393442917`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.4152014393750353`*^9}], Cell[BoxData["0.21`19.744727494896697"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.4152014394008713`*^9}], Cell[BoxData["0.21`19.744727494896697"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439434623*^9}], Cell[BoxData["0.21`19.744727494896697"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439469058*^9}], Cell[BoxData["0.21`19.744727494896697"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.4152014395019197`*^9}], Cell[BoxData["0.06`19.85387196432176"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.4152014395351143`*^9}], Cell[BoxData["0.09999999999999999999999999999992111391`19.52287874528034"], \ "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439568434*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"sum(w) = \"\>", "\[InvisibleSpace]", "0.99999999999999999999999999999992111391`19.72216166699796"}], SequenceForm["sum(w) = ", 1.`19.72216166699796], Editable->False]], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.4152014396021748`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439636169*^9}], Cell[BoxData["\<\"X\"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439669175*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.41520143972404*^9}], Cell[BoxData[ RowBox[{"-", "0.4879500364742665896771923181200500948`19.85387196432177"}]], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439747085*^9}], Cell[BoxData["0.4879500364742665896771923181200500948`19.85387196432177"], \ "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439770501*^9}], Cell[BoxData["0.4879500364742665896771923181200500948`19.85387196432177"], \ "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439802829*^9}], Cell[BoxData[ RowBox[{"-", "0.4879500364742665896771923181200500948`19.85387196432177"}]], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.4152014398355417`*^9}], Cell[BoxData["0"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439869429*^9}], Cell[BoxData["0"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439903351*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.4152014399365*^9}], Cell[BoxData["\<\"Y\"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201439970228*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201440003234*^9}], Cell[BoxData[ RowBox[{"-", "0.4879500364742665896771923181200500948`19.85387196432177"}]], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.4152014400364943`*^9}], Cell[BoxData[ RowBox[{"-", "0.4879500364742665896771923181200500948`19.85387196432177"}]], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.41520144007041*^9}], Cell[BoxData["0.4879500364742665896771923181200500948`19.85387196432177"], \ "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201440104073*^9}], Cell[BoxData["0.4879500364742665896771923181200500948`19.85387196432177"], \ "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201440137632*^9}], Cell[BoxData["0"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201440189817*^9}], Cell[BoxData["0"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201440215687*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201440240959*^9}], Cell[BoxData["\<\"Z\"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.4152014402704353`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201440304159*^9}], Cell[BoxData["0.16666666666666666666666666666666666666`19.698970004336026"], \ "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201440337188*^9}], Cell[BoxData["0.16666666666666666666666666666666666666`19.698970004336026"], \ "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201440370392*^9}], Cell[BoxData["0.16666666666666666666666666666666666666`19.698970004336026"], \ "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201440403458*^9}], Cell[BoxData["0.16666666666666666666666666666666666666`19.698970004336026"], \ "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.415201440436821*^9}], Cell[BoxData["0.58333333333333333333333333333333333334`20.84509804001426"], \ "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.41520144047019*^9}], Cell[BoxData["0.75`20.477121254719663"], "Print", CellChangeTimes->{3.413655075568932*^9, 3.41365558859274*^9, 3.413813742476877*^9, 3.413819704095929*^9, 3.4146738728359632`*^9, 3.4152014405037746`*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"rule", " ", "=", " ", "8"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"numer", " ", "=", " ", "True"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"wsum", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "rule"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"w", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"3", "/", "8"}], " ", ")"}], " ", "*", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], ")"}], "^", "2"}], " ", "*", " ", "w"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "w", " ", "]"}], ";", "\[IndentingNewLine]", RowBox[{"wsum", " ", "=", " ", RowBox[{"wsum", " ", "+", "w"}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "wsum"}], " ", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"x", " ", "=", " ", RowBox[{"a", " ", "*", " ", RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], " ", ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "x", " ", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", " ", "=", " ", RowBox[{"b", " ", "*", " ", RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], " ", ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "y", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "z", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}]}], "Input", CellChangeTimes->{{3.413654633924918*^9, 3.413654799625798*^9}, { 3.41365484902993*^9, 3.413654878240939*^9}, {3.413654957362526*^9, 3.413654963399206*^9}, {3.413655009030931*^9, 3.413655015807214*^9}, { 3.413655100840743*^9, 3.413655101204001*^9}, {3.4136555727306747`*^9, 3.413655577706349*^9}, {3.4138136587778482`*^9, 3.413813668023497*^9}, { 3.413813750364406*^9, 3.413813756006405*^9}, {3.414673894829178*^9, 3.4146738952988033`*^9}, {3.415201447176774*^9, 3.415201449792779*^9}}], Cell[CellGroupData[{ Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451118979*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Pyramid rule number \"\>", "\[InvisibleSpace]", "8"}], SequenceForm["Pyramid rule number ", 8], Editable->False]], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014511507874`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451193362*^9}], Cell[BoxData["\<\"W\"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451223928*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451257547*^9}], Cell[BoxData["0.07558941155986907293751728836823841301`19.92310009707811"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451283646*^9}], Cell[BoxData["0.07558941155986907293751728836823841301`19.92310009707811"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014513166733`*^9}], Cell[BoxData["0.07558941155986907293751728836823841301`19.92310009707811"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014513501863`*^9}], Cell[BoxData["0.07558941155986907293751728836823841301`19.92310009707811"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451382707*^9}], Cell[BoxData["0.17441058844013092706248271163176122711`19.73039842905131"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451417008*^9}], Cell[BoxData["0.17441058844013092706248271163176122711`19.73039842905131"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451450328*^9}], Cell[BoxData["0.17441058844013092706248271163176122711`19.73039842905131"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451483902*^9}], Cell[BoxData["0.17441058844013092706248271163176122711`19.73039842905131"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014515169783`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"sum(w) = \"\>", "\[InvisibleSpace]", "0.99999999999999999999999999999999856048`19.780203919166894"}], SequenceForm["sum(w) = ", 1.`19.780203919166894], Editable->False]], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451550683*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451605604*^9}], Cell[BoxData["\<\"X\"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451629469*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014516531*^9}], Cell[BoxData[ RowBox[{"-", "0.26318405556971359557121701304556403629`19.959850197981535"}]], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014516836433`*^9}], Cell[BoxData["0.26318405556971359557121701304556403629`19.959850197981535"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014517178917`*^9}], Cell[BoxData["0.26318405556971359557121701304556403629`19.959850197981535"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451751396*^9}], Cell[BoxData[ RowBox[{"-", "0.26318405556971359557121701304556403629`19.959850197981535"}]], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014517849827`*^9}], Cell[BoxData[ RowBox[{"-", "0.50661630334978742377431469429035839296`19.844606322115272"}]], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451819026*^9}], Cell[BoxData["0.50661630334978742377431469429035839296`19.844606322115272"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014518524637`*^9}], Cell[BoxData["0.50661630334978742377431469429035839296`19.844606322115272"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451885275*^9}], Cell[BoxData[ RowBox[{"-", "0.50661630334978742377431469429035839296`19.844606322115272"}]], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014519183607`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451951606*^9}], Cell[BoxData["\<\"Y\"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201451985264*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201452019301*^9}], Cell[BoxData[ RowBox[{"-", "0.26318405556971359557121701304556403629`19.959850197981535"}]], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201452052656*^9}], Cell[BoxData[ RowBox[{"-", "0.26318405556971359557121701304556403629`19.959850197981535"}]], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014521055403`*^9}], Cell[BoxData["0.26318405556971359557121701304556403629`19.959850197981535"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201452131917*^9}], Cell[BoxData["0.26318405556971359557121701304556403629`19.959850197981535"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014521698627`*^9}], Cell[BoxData[ RowBox[{"-", "0.50661630334978742377431469429035839296`19.844606322115265"}]], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014522027893`*^9}], Cell[BoxData[ RowBox[{"-", "0.50661630334978742377431469429035839296`19.844606322115265"}]], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201452235969*^9}], Cell[BoxData["0.50661630334978742377431469429035839296`19.844606322115265"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201452269699*^9}], Cell[BoxData["0.50661630334978742377431469429035839296`19.844606322115265"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014523031073`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201452336748*^9}], Cell[BoxData["\<\"Z\"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201452369637*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201452403078*^9}], Cell[BoxData["0.54415184401122528879992623629551456891`21.09077125854498"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201452435375*^9}], Cell[BoxData["0.54415184401122528879992623629551456891`21.09077125854498"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014524697*^9}], Cell[BoxData["0.54415184401122528879992623629551456891`21.09077125854498"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014525028267`*^9}], Cell[BoxData["0.54415184401122528879992623629551456891`21.09077125854498"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014525360527`*^9}], Cell[BoxData["0.12251482265544137786674043037115209775`19.51128873278299"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201452587761*^9}], Cell[BoxData["0.12251482265544137786674043037115209775`19.51128873278299"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201452614842*^9}], Cell[BoxData["0.12251482265544137786674043037115209775`19.51128873278299"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.415201452652521*^9}], Cell[BoxData["0.12251482265544137786674043037115209775`19.51128873278299"], \ "Print", CellChangeTimes->{{3.413655102334674*^9, 3.413655123050221*^9}, 3.413655578412269*^9, 3.413813756875557*^9, 3.41467389649409*^9, 3.4152014526857862`*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"rule", " ", "=", " ", RowBox[{"-", "8"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"numer", " ", "=", " ", "True"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"wsum", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "rule"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"w", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"3", "/", "8"}], " ", ")"}], " ", "*", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], ")"}], "^", "2"}], " ", "*", " ", "w"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "w", " ", "]"}], ";", "\[IndentingNewLine]", RowBox[{"wsum", " ", "=", " ", RowBox[{"wsum", " ", "+", "w"}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "wsum"}], " ", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"x", " ", "=", " ", RowBox[{"a", " ", "*", " ", RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], " ", ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "x", " ", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", " ", "=", " ", RowBox[{"b", " ", "*", " ", RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], " ", ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "y", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "z", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}]}], "Input", CellChangeTimes->{{3.413654633924918*^9, 3.413654799625798*^9}, { 3.41365484902993*^9, 3.413654878240939*^9}, {3.413654957362526*^9, 3.413654963399206*^9}, {3.413655009030931*^9, 3.413655015807214*^9}, { 3.413655132528431*^9, 3.413655133586656*^9}, {3.4136555608281193`*^9, 3.413655564906551*^9}, {3.413813763053733*^9, 3.4138137700452633`*^9}, { 3.4146739027007504`*^9, 3.414673903114731*^9}, {3.415201460404704*^9, 3.415201463361573*^9}}], Cell[CellGroupData[{ Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201464611731*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Pyramid rule number \"\>", "\[InvisibleSpace]", RowBox[{"-", "8"}]}], SequenceForm["Pyramid rule number ", -8], Editable->False]], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201464649694*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201464690445*^9}], Cell[BoxData["\<\"W\"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201464721826*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201464756331*^9}], Cell[BoxData["0.16438287736328777572426139601856958058`19.726745121078157"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201464789424*^9}], Cell[BoxData["0.16438287736328777572426139601856958058`19.726745121078157"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201464820859*^9}], Cell[BoxData["0.16438287736328777572426139601856958058`19.726745121078157"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.4152014648481207`*^9}], Cell[BoxData["0.16438287736328777572426139601856958058`19.726745121078157"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201464881425*^9}], Cell[BoxData["0.08561712263671222427573860398143049219`19.968149294271626"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201464914887*^9}], Cell[BoxData["0.08561712263671222427573860398143049219`19.968149294271626"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201464947358*^9}], Cell[BoxData["0.08561712263671222427573860398143049219`19.968149294271626"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201464981051*^9}], Cell[BoxData["0.08561712263671222427573860398143049219`19.968149294271626"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465015274*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"sum(w) = \"\>", "\[InvisibleSpace]", "1.00000000000000000000000000000000029109`19.79530483982384"}], SequenceForm["sum(w) = ", 1.`19.79530483982384], Editable->False]], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465070298*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465094179*^9}], Cell[BoxData["\<\"X\"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465117901*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465149042*^9}], Cell[BoxData[ RowBox[{"-", "0.51197009372656270106787177114855164201`19.842226734065253"}]], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465181823*^9}], Cell[BoxData["0.51197009372656270106787177114855164201`19.842226734065253"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.41520146521516*^9}], Cell[BoxData["0.51197009372656270106787177114855164201`19.842226734065253"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465249669*^9}], Cell[BoxData[ RowBox[{"-", "0.51197009372656270106787177114855164201`19.842226734065253"}]], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465282353*^9}], Cell[BoxData[ RowBox[{"-", "0.28415447557052037455572583630241398546`19.983782725344142"}]], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465317027*^9}], Cell[BoxData["0.28415447557052037455572583630241398546`19.983782725344142"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465350422*^9}], Cell[BoxData["0.28415447557052037455572583630241398546`19.983782725344142"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465382992*^9}], Cell[BoxData[ RowBox[{"-", "0.28415447557052037455572583630241398546`19.983782725344142"}]], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465416317*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465449728*^9}], Cell[BoxData["\<\"Y\"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465483526*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465517473*^9}], Cell[BoxData[ RowBox[{"-", "0.51197009372656270106787177114855164201`19.842226734065243"}]], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465569152*^9}], Cell[BoxData[ RowBox[{"-", "0.51197009372656270106787177114855164201`19.842226734065243"}]], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.41520146559543*^9}], Cell[BoxData["0.51197009372656270106787177114855164201`19.842226734065243"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.4152014656213083`*^9}], Cell[BoxData["0.51197009372656270106787177114855164201`19.842226734065243"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465650537*^9}], Cell[BoxData[ RowBox[{"-", "0.28415447557052037455572583630241398546`19.983782725344142"}]], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465684061*^9}], Cell[BoxData[ RowBox[{"-", "0.28415447557052037455572583630241398546`19.983782725344142"}]], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.41520146571736*^9}], Cell[BoxData["0.28415447557052037455572583630241398546`19.983782725344142"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465750602*^9}], Cell[BoxData["0.28415447557052037455572583630241398546`19.983782725344142"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.4152014657840557`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465817657*^9}], Cell[BoxData["\<\"Z\"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.4152014658500223`*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465886146*^9}], Cell[BoxData["0.11024490204163285720001714824664956347`19.451566711477998"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465921301*^9}], Cell[BoxData["0.11024490204163285720001714824664956347`19.451566711477998"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465951128*^9}], Cell[BoxData["0.11024490204163285720001714824664956347`19.451566711477998"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201465983471*^9}], Cell[BoxData["0.11024490204163285720001714824664956347`19.451566711477998"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201466035845*^9}], Cell[BoxData["0.51832652652979571422855428032477900796`21.451523275513946"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201466063242*^9}], Cell[BoxData["0.51832652652979571422855428032477900796`21.451523275513946"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201466101202*^9}], Cell[BoxData["0.51832652652979571422855428032477900796`21.451523275513946"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.4152014661337433`*^9}], Cell[BoxData["0.51832652652979571422855428032477900796`21.451523275513946"], \ "Print", CellChangeTimes->{3.41365513440967*^9, 3.4136555657308187`*^9, 3.4138137716125402`*^9, 3.413819768218162*^9, 3.4146739044374323`*^9, 3.415201466167201*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"rule", " ", "=", " ", "9"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"numer", " ", "=", " ", "True"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"wsum", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "rule"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"w", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"3", "/", "8"}], " ", ")"}], " ", "*", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], ")"}], "^", "2"}], " ", "*", " ", "w"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "w", " ", "]"}], ";", "\[IndentingNewLine]", RowBox[{"wsum", " ", "=", " ", RowBox[{"wsum", " ", "+", "w"}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "wsum"}], " ", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"x", " ", "=", " ", RowBox[{"a", " ", "*", " ", RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], " ", ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "x", " ", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", " ", "=", " ", RowBox[{"b", " ", "*", " ", RowBox[{"(", " ", RowBox[{"1", " ", "-", " ", "z"}], " ", ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "y", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\< \>\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", " ", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", RowBox[{"Abs", "[", "rule", "]"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", "w"}], "}"}], " ", "=", " ", RowBox[{"PyraGaussRuleInfo", "[", RowBox[{ RowBox[{"{", RowBox[{"rule", ",", "numer"}], "}"}], ",", "i"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"z", " ", "=", " ", RowBox[{ RowBox[{"(", " ", RowBox[{"c", " ", "+", " ", "1"}], " ", ")"}], " ", "/", " ", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", " ", "[", " ", "z", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}]}], "Input", CellChangeTimes->{{3.413654633924918*^9, 3.413654799625798*^9}, { 3.41365484902993*^9, 3.413654878240939*^9}, {3.413654957362526*^9, 3.413654963399206*^9}, {3.413655009030931*^9, 3.413655015807214*^9}, { 3.413655147967842*^9, 3.413655148322613*^9}, {3.413655547762664*^9, 3.413655553355095*^9}, {3.413813783432456*^9, 3.41381378869302*^9}, { 3.4146739126612387`*^9, 3.414673913050136*^9}, {3.415201474170848*^9, 3.415201476734393*^9}}], Cell[CellGroupData[{ Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.4152014772080317`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Pyramid rule number \"\>", "\[InvisibleSpace]", "9"}], SequenceForm["Pyramid rule number ", 9], Editable->False]], "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.41520147724074*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.415201477271002*^9}], Cell[BoxData["\<\"W\"\>"], "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.415201477296763*^9}], Cell[BoxData["\<\" \"\>"], "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.41520147732861*^9}], Cell[BoxData["0.1307338967227594479125121507626381027`19.71875074396513"], \ "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.415201477362423*^9}], Cell[BoxData["0.1307338967227594479125121507626381027`19.71875074396513"], \ "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.415201477396188*^9}], Cell[BoxData["0.1307338967227594479125121507626381027`19.71875074396513"], \ "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.415201477429695*^9}], Cell[BoxData["0.1307338967227594479125121507626381027`19.71875074396513"], \ "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.415201477462536*^9}], Cell[BoxData["0.10989110327724055208748784923735608703`19.875400077308992"], \ "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.415201477495598*^9}], Cell[BoxData["0.10989110327724055208748784923735608703`19.875400077308992"], \ "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.415201477529195*^9}], Cell[BoxData["0.10989110327724055208748784923735608703`19.875400077308992"], \ "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.4152014775628757`*^9}], Cell[BoxData["0.10989110327724055208748784923735608703`19.875400077308992"], \ "Print", CellChangeTimes->{3.413655149293708*^9, 3.413655554481958*^9, 3.413813789705792*^9, 3.414673617993634*^9, 3.414673913657628*^9, 3.41520147759686*^9}], C