(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.0' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 30645, 795]*) (*NotebookOutlinePosition[ 31291, 817]*) (* CellTagsIndexPosition[ 31247, 813]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ \(Solve[ x\^m == \(\(\((\(-1\))\)\^\(n + m\)\) \@\(2 m\)\)\/\(2 m\^2 - m \((2 \ n - 1)\) - m \((2 k + 1)\) + \(1\/2\) \((2 n - 1)\) \((2 k + 1)\)\), x]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{\(Solve::"ifun"\), \(\(:\)\(\ \)\), "\<\"Inverse functions are \ being used by \\!\\(TraditionalForm\\`Solve\\), so some solutions may not be \ found; use Reduce for complete solution information. \ \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\ \\\", ButtonFrame->None, ButtonData:>\\\"Solve::ifun\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ \(TraditionalForm\`{{x \[Rule] 2\^\(1\/\(2\ m\)\)\ \((\(-\(\(\((\(-1\))\)\^\(m + n\)\ \ \@m\)\/\(\(-2\)\ m\^2 + 2\ k\ m + 2\ n\ m + k - 2\ k\ n - n + 1\/2\)\)\))\)\^\ \(1\/m\)}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Solve[ x\^m == \(\(\((\(-1\))\)\^\(n + m\)\) \@\(2 m\)\)\/\(2 m\^2 + m \((2 \ n - 1)\) + m \((2 k - 1)\) + \(1\/2\) \((2 n - 1)\) \((2 k - 1)\)\), x]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{\(Solve::"ifun"\), \(\(:\)\(\ \)\), "\<\"Inverse functions are \ being used by \\!\\(TraditionalForm\\`Solve\\), so some solutions may not be \ found; use Reduce for complete solution information. \ \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\ \\\", ButtonFrame->None, ButtonData:>\\\"Solve::ifun\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ \(TraditionalForm\`{{x \[Rule] 2\^\(1\/\(2\ m\)\)\ \((\(\((\(-1\))\)\^\(m + n\)\ \@m\)\/\(2\ m\^2 \ + 2\ k\ m + 2\ n\ m - 2\ m - k + 2\ k\ n - n + 1\/2\))\)\^\(1\/m\)}}\)], \ "Output"] }, Open ]], Cell["\<\ Trying to verify or otherwise generate a value for these sums.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(\(Clear[n, \ k, p];\)\[IndentingNewLine] \(k := 3;\)\[IndentingNewLine] \(p := \ 3;\)\[IndentingNewLine] \(n := 2;\)\[IndentingNewLine] \(F[x_]\ := \ \((\(1 + x\)\/\(1 - x\))\)\^\(1\/p\);\)\[IndentingNewLine] \(\((\(-1\))\)\^n\) \((Re[\((\(z\^\(2 k - 1\)\) \(\[Integral]\(z\^\(\(-2\) k + 2 n + 1\)\) \((\[Integral]\(z\^\(\(-2\) n - 2\)\) F[z] \[DifferentialD]z)\) \[DifferentialD]z\))\) /. z \[Rule] \[ImaginaryI]] - \((\((\(z\^\(2 k - 1\)\) \(\[Integral]\(z\^\(\(-2\) k + 2 n + 1\)\) \((\[Integral]\(z\^\(\(-2\) n - 2\)\) F[z] \[DifferentialD]z)\) \[DifferentialD]z\))\) /. z \[Rule] 0)\))\)\)\)\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{\(Integrate::"ilim"\), \(\(:\)\(\ \)\), "\<\"Invalid \ integration variable or limit(s) in \\!\\(TraditionalForm\\`\[ImaginaryI]\\). \ \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\ \\\", ButtonFrame->None, ButtonData:>\\\"Integrate::ilim\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ FormBox[ RowBox[{\(Power::"infy"\), \(\(:\)\(\ \)\), "\<\"Infinite expression \ \\!\\(TraditionalForm\\`1\\/0\\^5\\) encountered. \\!\\(\\*ButtonBox[\\\"More\ \[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"Power::infy\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ FormBox[ RowBox[{\(\[Infinity]::"indet"\), \(\(:\)\(\ \)\), "\<\"Indeterminate \ expression \\!\\(TraditionalForm\\`\\(0\\\\ \\@2\\%3\\\\ \[Infinity]\\)\\) \ encountered. \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::indet\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ FormBox[ RowBox[{\(\[Infinity]::"indet"\), \(\(:\)\(\ \)\), "\<\"Indeterminate \ expression \\!\\(TraditionalForm\\`\\(0\\\\ \\@2\\%3\\\\ \[Infinity]\\)\\) \ encountered. \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::indet\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ FormBox[ RowBox[{\(Integrate::"ilim"\), \(\(:\)\(\ \)\), "\<\"Invalid \ integration variable or limit(s) in \\!\\(TraditionalForm\\`0\\). \ \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\ \\\", ButtonFrame->None, ButtonData:>\\\"Integrate::ilim\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ FormBox[ RowBox[{"-", RowBox[{"Im", "(", RowBox[{"\[Integral]", RowBox[{ RowBox[{\(1\/2430\), RowBox[{"(", RowBox[{\(\@\(-1\)\%6\), " ", \(\((1 + \[ImaginaryI])\)\^\(2/3\)\), " ", RowBox[{"(", RowBox[{\(\((820 - 186\ \[ImaginaryI])\)\ \@\(1 + \ \[ImaginaryI]\)\%3\), "+", RowBox[{ "335", " ", "\[ImaginaryI]", " ", \(\@2\%3\), " ", RowBox[{ InterpretationBox[\(F\_1\), AppellF1, Editable->False, Selectable->False], "(", \(2\/3; \(-\(1\/3\)\), 1; 5\/3; 1\/2 - \[ImaginaryI]\/2, 1 - \[ImaginaryI]\), ")"}]}], "+", RowBox[{\((67 + 67\ \[ImaginaryI])\), " ", \(\@2\%3\), " ", RowBox[{ InterpretationBox[\(F\_1\), AppellF1, Editable->False, Selectable->False], "(", \(5\/3; 2\/3, 1; 8\/3; 1\/2 - \[ImaginaryI]\/2, 1 - \[ImaginaryI]\), ")"}]}]}], ")"}]}], ")"}]}], \(\[DifferentialD]\[ImaginaryI]\)}]}], ")"}]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(FullSimplify[\((\(z\^\(2 k - 1\)\) \(\[Integral]\(z\^\(\(-2\) k + 2 n + 1\)\) \((\[Integral]\(z\^\(\(-2\) n - 2\)\) F[z] \[DifferentialD]z)\) \[DifferentialD]z\))\)]\)], \ "Input"], Cell[BoxData[ FormBox[ RowBox[{\(z\^5\), " ", RowBox[{"\[Integral]", RowBox[{ RowBox[{"-", RowBox[{\(1\/\(2430\ z\^6\)\), RowBox[{"(", RowBox[{\(\@\(1\/\(1 - z\)\)\%3\), " ", \((z - 1)\), " ", RowBox[{"(", RowBox[{ RowBox[{"67", " ", \(\@2\%3\), " ", \(z\^5\), " ", RowBox[{"(", RowBox[{ RowBox[{\((z - 1)\), " ", RowBox[{ InterpretationBox[\(F\_1\), AppellF1, Editable->False, Selectable->False], "(", \(5\/3; 2\/3, 1; 8\/3; \(1 - z\)\/2, 1 - z\), ")"}]}], "-", RowBox[{"5", " ", RowBox[{ InterpretationBox[\(F\_1\), AppellF1, Editable->False, Selectable->False], "(", \(2\/3; \(-\(1\/3\)\), 1; 5\/3; \(1 - z\)\/2, 1 - z\), ")"}]}]}], ")"}]}], "-", \(\@\(z + 1\)\%3\ \((z\ \((z\ \((z\ \((919\ z + 753)\) + 585)\) + 567)\) + 486)\)\)}], ")"}]}], ")"}]}]}], \(\[DifferentialD]z\)}]}]}], TraditionalForm]], "Output"] }, Open ]], Cell["\<\ Let me to try this one step at a time. Verify (W-7).\ \>", "Text"], Cell[BoxData[ \(\(V[n_, \ p_] := \ \(1\/Factorial[n]\) Limit[D[\((\(1\ + \ x\)\/\(1\ - \ x\))\)\^\(1\/p\), \ {x, \ n}], \ x \[Rule] 0];\)\)], "Input"], Cell["\<\ Alternating series--some end up below the result below, some above. So far, \ it looks good!\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(N[\[Sum]\+\(m = 1\)\%14\(\((\(-1\))\)\^m\) V[2 m, \ 3]]\)\)\)], "Input"], Cell[BoxData[ \(TraditionalForm\`\(-0.11017180748936595`\)\)], "Output"] }, Open ]], Cell["\<\ Very good.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(N[ Re[\((\(1 + \[ImaginaryI]\)\/\(1 - \[ImaginaryI]\))\)\^\(1\/3\)] - \ \((\(1 + 0\)\/\(1 - 0\))\)\^\(1\/3\)]\)\)\)], "Input"], Cell[BoxData[ \(TraditionalForm\`\(-0.1339745962155614`\)\)], "Output"] }, Open ]], Cell["\<\ Now what if I considered just one of my two 1/x factors form (W-8). It looks \ like this converges really quickly! \ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Clear[n];\)\), "\[IndentingNewLine]", \(\(n\ := \ 3;\)\), "\[IndentingNewLine]", \(N[\[Sum]\+\(m = 1\)\%16\(\((\(-1\))\)\^m\/\(2 m - \((2 n + 1)\)\)\) V[2 m, \ 3]]\)}], "Input"], Cell[BoxData[ \(TraditionalForm\`0.2027028364602475`\)], "Output"] }, Open ]], Cell["\<\ Now, the issue with creating this through integration is that I need to kill \ off the odd Taylor components before integrating, not after. Note that I \ need the power on each particular x to be equal to (2m - 2n + 1) - 1. In \ other words, I need to adjust each one by -2n. To be fair, however, the \ primary issue at hand is only the Taylor term that will correspond to 1/x \ before integration. If I subtracted out that term then integrated, divided \ and evaluated, it would give the same result. Let me give it a shot. Start then with my series V(1/p,x). Multiply it by x^(-2n) Subtract x^-1 times the ( 2n - 1) th Taylor coefficient. Integrate Divide it by x^(-2n+1) to adjust it again Then take the real part evaluated at x=i.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Clear[n, p];\)\), "\[IndentingNewLine]", \(\(p := 3;\)\), "\[IndentingNewLine]", \(\(n\ := \ 3;\)\), "\[IndentingNewLine]", \(N[Re[\(\[Integral]\((\(x\^\(\(-2\) n\)\) \((\(1\ + \ x\)\/\(1\ - \ \ x\))\)\^\(1\/p\) - \(x\^\(-1\)\) V[2 n - 1, \ p])\) \[DifferentialD]x\)\/x\^\ \(\(-2\) n + 1\) /. x \[Rule] \[ImaginaryI]] - \((\(\[Integral]\((\(x\^\(\(-2\) n\)\) \ \((\(1\ + \ x\)\/\(1\ - \ x\))\)\^\(1\/p\) - \(x\^\(-1\)\) V[2 n - 1, \ \ p])\) \[DifferentialD]x\)\/x\^\(\(-2\) n + 1\) /. x \[Rule] 0)\)]\), "\[IndentingNewLine]", \(\)}], "Input"], Cell[BoxData[ FormBox[ RowBox[{\(Power::"infy"\), \(\(:\)\(\ \)\), "\<\"Infinite expression \ \\!\\(TraditionalForm\\`1\\/0\\^5\\) encountered. \\!\\(\\*ButtonBox[\\\"More\ \[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"Power::infy\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ FormBox[ RowBox[{\(\[Infinity]::"indet"\), \(\(:\)\(\ \)\), "\<\"Indeterminate \ expression \\!\\(TraditionalForm\\`\\(ComplexInfinity + \\(\\(-\[Infinity]\\)\ \\) + \[Infinity]\\)\\) encountered. \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\ \", ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::indet\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ \(TraditionalForm\`Indeterminate\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(\(Clear[n, p];\)\[IndentingNewLine] \(p := 3;\)\[IndentingNewLine] \(n\ := \ 3;\)\[IndentingNewLine] Series[\((\(1\ + \ x\)\/\(1\ - \ x\))\)\^\(1\/p\), \ {x, \ 0, \ 5}]\[IndentingNewLine] V[2 n - 1, p]\[IndentingNewLine] Series[\(x\^\(\(-2\) n - 2\)\) \((\(1\ + \ x\)\/\(1\ - \ x\))\)\^\(1\/p\), \ {x, \ 0, \ 5}]\[IndentingNewLine] Series[\(x\^\(\(-2\) n - 2\)\) \((\(1\ + \ x\)\/\(1\ - \ x\))\)\^\(1\/p\) - \ \(x\^\(-1\)\) V[2 n + 1, \ p], \ {x, \ 0, \ 5}]\)\)\)], "Input"], Cell[BoxData[ FormBox[ InterpretationBox[ RowBox[{ "1", "+", \(\(2\ x\)\/3\), "+", \(\(2\ x\^2\)\/9\), "+", \(\(22\ x\^3\)\/81\), "+", \(\(38\ x\^4\)\/243\), "+", \(\(134\ x\^5\)\/729\), "+", InterpretationBox[\(O(x\^6)\), SeriesData[ x, 0, {}, 0, 6, 1], Editable->False]}], SeriesData[ x, 0, {1, Rational[ 2, 3], Rational[ 2, 9], Rational[ 22, 81], Rational[ 38, 243], Rational[ 134, 729]}, 0, 6, 1], Editable->False], TraditionalForm]], "Output"], Cell[BoxData[ \(TraditionalForm\`134\/729\)], "Output"], Cell[BoxData[ FormBox[ InterpretationBox[ RowBox[{\(1\/x\^8\), "+", \(2\/\(3\ x\^7\)\), "+", \(2\/\(9\ x\^6\)\), "+", \(22\/\(81\ x\^5\)\), "+", \(38\/\(243\ x\^4\)\), "+", \(134\/\(729\ x\^3\)\), "+", \(818\/\(6561\ x\^2\)\), "+", \(2818\/\(19683\ x\)\), "+", \(6226\/59049\), "+", \(\(189986\ x\)\/1594323\), "+", \(\(441442\ x\^2\)\/4782969\), "+", \(\(1479250\ x\^3\)\/14348907\), "+", \(\(10672070\ x\^4\)\/129140163\), "+", \(\(35437030\ x\^5\)\/387420489\), "+", InterpretationBox[\(O(x\^6)\), SeriesData[ x, 0, {}, -8, 6, 1], Editable->False]}], SeriesData[ x, 0, {1, Rational[ 2, 3], Rational[ 2, 9], Rational[ 22, 81], Rational[ 38, 243], Rational[ 134, 729], Rational[ 818, 6561], Rational[ 2818, 19683], Rational[ 6226, 59049], Rational[ 189986, 1594323], Rational[ 441442, 4782969], Rational[ 1479250, 14348907], Rational[ 10672070, 129140163], Rational[ 35437030, 387420489]}, -8, 6, 1], Editable->False], TraditionalForm]], "Output"], Cell[BoxData[ FormBox[ InterpretationBox[ RowBox[{\(1\/x\^8\), "+", \(2\/\(3\ x\^7\)\), "+", \(2\/\(9\ x\^6\)\), "+", \(22\/\(81\ x\^5\)\), "+", \(38\/\(243\ x\^4\)\), "+", \(134\/\(729\ x\^3\)\), "+", \(818\/\(6561\ x\^2\)\), "+", \(6226\/59049\), "+", \(\(189986\ x\)\/1594323\), "+", \(\(441442\ x\^2\)\/4782969\), "+", \(\(1479250\ x\^3\)\/14348907\), "+", \(\(10672070\ x\^4\)\/129140163\), "+", \(\(35437030\ x\^5\)\/387420489\), "+", InterpretationBox[\(O(x\^6)\), SeriesData[ x, 0, {}, -8, 6, 1], Editable->False]}], SeriesData[ x, 0, {1, Rational[ 2, 3], Rational[ 2, 9], Rational[ 22, 81], Rational[ 38, 243], Rational[ 134, 729], Rational[ 818, 6561], 0, Rational[ 6226, 59049], Rational[ 189986, 1594323], Rational[ 441442, 4782969], Rational[ 1479250, 14348907], Rational[ 10672070, 129140163], Rational[ 35437030, 387420489]}, -8, 6, 1], Editable->False], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(N[ Re[\((\(x\^\(2 n + 1\)\) \(\[Integral]\((\(x\^\(\(-2\) n - 2\)\) \((\(1\ + \ x\)\/\(1\ - \ \ x\))\)\^\(1\/p\) - \(x\^\(-1\)\) V[2 n + 1, \ p])\) \[DifferentialD]x\))\) /. x \[Rule] \[ImaginaryI]]]\)\)\)], "Input"], Cell[BoxData[ \(TraditionalForm\`0.059004324210751476`\)], "Output"] }, Open ]], Cell["\<\ So this looks really, really good.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(N[\[Sum]\+\(m = \ 0\)\%15\(\((\(-1\))\)\^m\/\(2 m - \((2 n + 1)\)\)\) V[2 m, \ 3]]\)\)\)], "Input"], Cell[BoxData[ \(TraditionalForm\`0.058056145861710026`\)], "Output"] }, Open ]], Cell["\<\ This MUST give me the difference, when evaluated at x=0.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(FullSimplify[\(x\^\(2 n + 1\)\) \(\[Integral]\((\(x\^\(\(-2\) n - 2\)\) \((\(1\ + \ x\)\/\(1\ - \ x\))\)\^\(1\/p\) - \ \(x\^\(-1\)\) V[2 n + 1, \ p])\) \[DifferentialD]x\)]\)\)\)], "Input"], Cell[BoxData[ FormBox[ FractionBox[ RowBox[{\(x\^7\), " ", RowBox[{"(", RowBox[{\(\(3\ \((x - 1)\)\ \@\(\(x + 1\)\/\(1 - x\)\)\%3\ \((x\ \((x\ \((x\ \ \((x\ \((x\ \((140573\ x + 112611)\) + 89055)\) + 84969)\) + 73872)\) + 72900)\) + 65610)\)\)\/x\^7\), "-", RowBox[{ "29589", " ", \(\@2\%3\), " ", \(\@\(1\/\(1 - x\)\)\%3\), " ", \((x - 1)\), " ", RowBox[{"(", RowBox[{ RowBox[{\((x - 1)\), " ", RowBox[{ InterpretationBox[\(F\_1\), AppellF1, Editable->False, Selectable->False], "(", \(5\/3; 2\/3, 1; 8\/3; \(1 - x\)\/2, 1 - x\), ")"}]}], "-", RowBox[{"5", " ", RowBox[{ InterpretationBox[\(F\_1\), AppellF1, Editable->False, Selectable->False], "(", \(2\/3; \(-\(1\/3\)\), 1; 5\/3; \(1 - x\)\/2, 1 - x\), ")"}]}]}], ")"}]}], "-", \(197260\ \(log(x)\)\)}], ")"}]}], "1377810"], TraditionalForm]], "Output"] }, Open ]], Cell["\<\ Fine. It doesn't want to simplify this to a state where I can evaluate it at \ x=0. I can, however, examine the series above and see that for this \ particular choice of series, I always end up with an extra -1/7.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Clear[n, p];\)\), "\[IndentingNewLine]", \(\(p := 3;\)\), "\[IndentingNewLine]", \(\(n\ := \ 3;\)\), "\[IndentingNewLine]", \(N[Re[\((\(x\^\(2 n + 1\)\) \(\[Integral]\((\(x\^\(\(-2\) n - 2\)\) \((\(1\ + \ x\)\/\(1\ - \ x\))\)\^\(1\ \/p\) - \(x\^\(-1\)\) V[2 n + 1, \ p])\) \[DifferentialD]x\))\) /. x \[Rule] \[ImaginaryI]]\ + \ 1\/\(2\ n\ + \ 1\)]\)}], "Input"], Cell[BoxData[ \(TraditionalForm\`0.20186146706789432`\)], "Output"] }, Open ]], Cell["\<\ Very good. I see then that this works for the addition of one of the series. \ Now I should be able to repeat the procedure above for the other series. \ First, though, let me see if this is simplifiable.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(Clear[n, p]\[IndentingNewLine] FullSimplify[ Re[Limit[\((\(x\^\(2 n + 1\)\) \(\[Integral]\((\(x\^\(\(-2\) n - 2\)\) \((\(1\ + \ x\)\/\(1\ - \ x\))\)\^\(1\ \/p\) - \(x\^\(-1\)\) V[2 n + 1, \ p])\) \[DifferentialD]x\))\), x \[Rule] \[ImaginaryI]]]\ + \ 1\/\(2\ n\ + \ 1\), \ {{n, p} \[Epsilon]Integers, \ n \[GreaterEqual] 0, \ p > 0}]\)\)\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{\(General::"ivar"\), \(\(:\)\(\ \)\), \ "\<\"\\!\\(TraditionalForm\\`\\(\[ImaginaryI] + \ \[ExponentialE]\\^\\(Limit`t(40525)\\)\\)\\) is not a valid variable. \ \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\ \\\", ButtonFrame->None, ButtonData:>\\\"General::ivar\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ FormBox[ RowBox[{\(General::"ivar"\), \(\(:\)\(\ \)\), \ "\<\"\\!\\(TraditionalForm\\`\\(\[ImaginaryI] + \ \[ExponentialE]\\^\\(Limit`t(40525)\\)\\)\\) is not a valid variable. \ \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\ \\\", ButtonFrame->None, ButtonData:>\\\"General::ivar\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ FormBox[ RowBox[{\(Integrate::"ilim"\), \(\(:\)\(\ \)\), "\<\"Invalid \ integration variable or limit(s) in \\!\\(TraditionalForm\\`\\(\[ImaginaryI] \ + \[ExponentialE]\\^\\(Limit`t(40525)\\)\\)\\). \\!\\(\\*ButtonBox[\\\"More\ \[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"Integrate::ilim\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{"Re", "(", RowBox[{\(lim\+\(x \[Rule] \[ImaginaryI]\)\), "\[ThinSpace]", RowBox[{\(x\^\(2\ n + 1\)\), " ", RowBox[{"\[Integral]", RowBox[{ FractionBox[ RowBox[{\(x\^\(\(-2\)\ n\)\ \((\(x + 1\)\/\(1 - \ x\))\)\^\(1\/p\)\), "-", FractionBox[ RowBox[{"x", " ", RowBox[{"(", RowBox[{\(lim\+\(x \[Rule] 0\)\), "\[ThinSpace]", FractionBox[\(\[PartialD]\^\(2\ n + 1\)\((\(x + \ 1\)\/\(1 - x\))\)\^\(1\/p\)\), \(\[PartialD]x\^\(2\ n + 1\)\), MultilineFunction->None]}], ")"}]}], \(\[CapitalGamma]( 2\ \((n + 1)\))\)]}], \(x\^2\)], \ \(\[DifferentialD]x\)}]}]}]}], ")"}], "+", \(1\/\(2\ n + 1\)\)}], TraditionalForm]], "Output"] }, Open ]], Cell["\<\ Fine. It doesn't want to simplify this. At least I know it works--maybe Dr. \ A knows how to simplify this.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Clear[n];\)\), "\[IndentingNewLine]", \(\(n\ := \ 3;\)\), "\[IndentingNewLine]", \(\(k := 2;\)\), "\[IndentingNewLine]", \(\(p := 3;\)\), "\[IndentingNewLine]", \(N[\[Sum]\+\(m = 0\)\%16\( 1\/\(2 m - 2 k + 1\)\) \(\((\(-1\))\)\^m\/\(2 m - \((2 n + 1)\)\)\) V[2 m, \ 3]]\), "\[IndentingNewLine]", \(N[\[Sum]\+\(m = 1\)\%16\( 1\/\(2 m - 2 k + 1\)\) \(\((\(-1\))\)\^m\/\(2 m - \((2 n + 1)\)\)\) V[2 m, \ 3]]\)}], "Input"], Cell[BoxData[ \(TraditionalForm\`0.010547958346988526`\)], "Output"], Cell[BoxData[ \(TraditionalForm\`\(-0.0370710892720591`\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(N[Re[\((\(x\^\(2 k - 1\)\) \(\[Integral]\((\(x\^\(\(-2\) k\)\) \(x\^\(2 n + 1\)\) \(\[Integral]\((\(x\^\(\(-2\) n - 2\)\) \((\(1\ + \ x\)\/\(1\ - \ \ x\))\)\^\(1\/p\) - \(x\^\(-1\)\) V[2 n + 1, \ p])\) \[DifferentialD]x\))\) \ \[DifferentialD]x\))\)]\ /. x \[Rule] \[ImaginaryI]]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{\(Integrate::"ilim"\), \(\(:\)\(\ \)\), "\<\"Invalid \ integration variable or limit(s) in \\!\\(TraditionalForm\\`\[ImaginaryI]\\). \ \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\ \\\", ButtonFrame->None, ButtonData:>\\\"Integrate::ilim\\\"]\\)\"\>"}], TraditionalForm]], "Message"], Cell[BoxData[ FormBox[ RowBox[{"7.257894775041551`*^-7", " ", RowBox[{"Im", "(", RowBox[{"\[Integral]", RowBox[{ RowBox[{\(-\[ImaginaryI]\), " ", RowBox[{"(", RowBox[{\(\((480966 - 122286\ \[ImaginaryI])\)\ \@\(-1\)\%6\), "-", \(98630\ \[ImaginaryI]\ \[Pi]\), "+", FractionBox[ RowBox[{ "59178", " ", \(\@\(-1\)\%6\), " ", \(\@2\%3\), " ", RowBox[{"(", RowBox[{ RowBox[{\(-5\), " ", RowBox[{ InterpretationBox[\(F\_1\), AppellF1, Editable->False, Selectable->False], "(", \(2\/3; \(-\(1\/3\)\), 1; 5\/3; 1\/2 - \[ImaginaryI]\/2, 1 - \[ImaginaryI]\), ")"}]}], "-", RowBox[{\((1 - \[ImaginaryI])\), " ", RowBox[{ InterpretationBox[\(F\_1\), AppellF1, Editable->False, Selectable->False], "(", \(5\/3; 2\/3, 1; 8\/3; 1\/2 - \[ImaginaryI]\/2, 1 - \[ImaginaryI]\), ")"}]}]}], ")"}]}], \(\((1 + \[ImaginaryI])\)\^\(4/3\)\)]}], ")"}]}], \(\[DifferentialD]\[ImaginaryI]\)}]}], ")"}]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(FullSimplify[ Re[\((\(x\^\(2 k - 1\)\) \(\[Integral]\((\(x\^\(\(-2\) k\)\) \(x\^\(2 n + 1\)\) \(\[Integral]\((\(x\^\(\(-2\) n - 2\)\) \((\(1\ + \ x\)\/\(1\ - \ \ x\))\)\^\(1\/p\) - \(x\^\(-1\)\) V[2 n + 1, \ p])\) \[DifferentialD]x\))\) \ \[DifferentialD]x\))\)]]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{\(1\/1377810\), RowBox[{"Re", "(", RowBox[{\(x\^3\), " ", RowBox[{"\[Integral]", RowBox[{ RowBox[{\(x\^3\), " ", RowBox[{"(", RowBox[{\(\(3\ \((x - 1)\)\ \@\(\(x + 1\)\/\(1 - x\)\)\%3\ \((x\ \ \((x\ \((x\ \((x\ \((x\ \((140573\ x + 112611)\) + 89055)\) + 84969)\) + 73872)\) + 72900)\) + 65610)\)\)\/x\^7\), "-", RowBox[{ "29589", " ", \(\@2\%3\), " ", \(\@\(1\/\(1 - x\)\)\%3\), " ", \((x - 1)\), " ", RowBox[{"(", RowBox[{ RowBox[{\((x - 1)\), " ", RowBox[{ InterpretationBox[\(F\_1\), AppellF1, Editable->False, Selectable->False], "(", \(5\/3; 2\/3, 1; 8\/3; \(1 - x\)\/2, 1 - x\), ")"}]}], "-", RowBox[{"5", " ", RowBox[{ InterpretationBox[\(F\_1\), AppellF1, Editable->False, Selectable->False], "(", \(2\/3; \(-\(1\/3\)\), 1; 5\/3; \(1 - x\)\/2, 1 - x\), ")"}]}]}], ")"}]}], "-", \(197260\ \(log(x)\)\)}], ")"}]}], \(\[DifferentialD]x\)}]}]}], ")"}]}], TraditionalForm]], "Output"] }, Open ]], Cell["\<\ I'm pretty sure the above integral is right. I just need to figure out how \ to evaluate them (the ones above have to have the zero-mode term subtracted \ out. In this particular expansion, it will have subtracted from it the value \ of the coefficient on the m=0 item! \ \>", "Text"], Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(Re[\((\(x\^\(2 k - 1\)\) \(\[Integral]\((\(x\^\(\(-2\) k\)\) \(x\^\(2 n + 1\)\) \(\[Integral]\((\(x\^\(\(-2\) n - 2\)\) \((\(1\ + \ x\)\/\(1\ - \ \ x\))\)\^\(1\/p\) - \(x\^\(-1\)\) V[2 n + 1, \ p])\) \[DifferentialD]x\))\) \ \[DifferentialD]x\))\) /. x \[Rule] \[ImaginaryI]] - \(1\/\(-\((2 n + 1)\)\)\) 1\/\(\(-2\) k + 1\)\)\)\)], "Input"] }, FrontEndVersion->"5.0 for Microsoft Windows", ScreenRectangle->{{0, 1280}, {0, 951}}, WindowSize->{1109, 879}, WindowMargins->{{14, Automatic}, {Automatic, 15}} ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1776, 53, 201, 4, 57, "Input"], Cell[1980, 59, 424, 7, 23, "Message"], Cell[2407, 68, 220, 4, 67, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2664, 77, 201, 4, 57, "Input"], Cell[2868, 83, 424, 7, 23, "Message"], Cell[3295, 92, 213, 4, 67, "Output"] }, Open ]], Cell[3523, 99, 87, 3, 52, "Text"], Cell[CellGroupData[{ Cell[3635, 106, 877, 14, 182, "Input"], Cell[4515, 122, 364, 6, 23, "Message"], Cell[4882, 130, 333, 5, 37, "Message"], Cell[5218, 137, 377, 6, 29, "Message"], Cell[5598, 145, 377, 6, 29, "Message"], Cell[5978, 153, 352, 6, 23, "Message"], Cell[6333, 161, 1639, 36, 45, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8009, 202, 261, 5, 41, "Input"], Cell[8273, 209, 1672, 33, 55, "Output"] }, Open ]], Cell[9960, 245, 78, 3, 52, "Text"], Cell[10041, 250, 184, 3, 49, "Input"], Cell[10228, 255, 118, 4, 52, "Text"], Cell[CellGroupData[{ Cell[10371, 263, 131, 2, 73, "Input"], Cell[10505, 267, 76, 1, 29, "Output"] }, Open ]], Cell[10596, 271, 35, 3, 52, "Text"], Cell[CellGroupData[{ Cell[10656, 278, 179, 3, 70, "Input"], Cell[10838, 283, 75, 1, 29, "Output"] }, Open ]], Cell[10928, 287, 144, 4, 52, "Text"], Cell[CellGroupData[{ Cell[11097, 295, 225, 4, 93, "Input"], Cell[11325, 301, 70, 1, 29, "Output"] }, Open ]], Cell[11410, 305, 770, 21, 318, "Text"], Cell[CellGroupData[{ Cell[12205, 330, 619, 11, 173, "Input"], Cell[12827, 343, 333, 5, 37, "Message"], Cell[13163, 350, 403, 6, 23, "Message"], Cell[13569, 358, 64, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13670, 364, 604, 12, 233, "Input"], Cell[14277, 378, 602, 16, 47, "Output"], Cell[14882, 396, 59, 1, 45, "Output"], Cell[14944, 399, 1237, 28, 48, "Output"], Cell[16184, 429, 1176, 27, 48, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17397, 461, 336, 6, 74, "Input"], Cell[17736, 469, 72, 1, 29, "Output"] }, Open ]], Cell[17823, 473, 59, 3, 52, "Text"], Cell[CellGroupData[{ Cell[17907, 480, 160, 3, 73, "Input"], Cell[18070, 485, 72, 1, 29, "Output"] }, Open ]], Cell[18157, 489, 81, 3, 52, "Text"], Cell[CellGroupData[{ Cell[18263, 496, 280, 4, 74, "Input"], Cell[18546, 502, 1364, 30, 87, "Output"] }, Open ]], Cell[19925, 535, 239, 5, 52, "Text"], Cell[CellGroupData[{ Cell[20189, 544, 476, 8, 114, "Input"], Cell[20668, 554, 71, 1, 29, "Output"] }, Open ]], Cell[20754, 558, 232, 5, 52, "Text"], Cell[CellGroupData[{ Cell[21011, 567, 500, 9, 94, "Input"], Cell[21514, 578, 389, 7, 25, "Message"], Cell[21906, 587, 389, 7, 25, "Message"], Cell[22298, 596, 406, 6, 25, "Message"], Cell[22707, 604, 1144, 26, 162, "Output"] }, Open ]], Cell[23866, 633, 135, 4, 52, "Text"], Cell[CellGroupData[{ Cell[24026, 641, 563, 12, 178, "Input"], Cell[24592, 655, 72, 1, 29, "Output"], Cell[24667, 658, 75, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[24779, 664, 475, 8, 52, "Input"], Cell[25257, 674, 364, 6, 23, "Message"], Cell[25624, 682, 1788, 36, 68, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[27449, 723, 454, 9, 52, "Input"], Cell[27906, 734, 1849, 38, 185, "Output"] }, Open ]], Cell[29770, 775, 298, 6, 71, "Text"], Cell[30071, 783, 570, 10, 74, "Input"] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)