(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 3.0, MathReader 3.0, or any compatible application. The data for the notebook starts with the line of 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[ 36000, 1139]*) (*NotebookOutlinePosition[ 36711, 1164]*) (* CellTagsIndexPosition[ 36667, 1160]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ \(CrN = Series[BesselJ[1 - \[Nu], x*\[Beta]], {\[Beta], 0, 4}]\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(1 - \[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(2\^\(\(-1\) + \[Nu]\)\ x\^\(1 - \[Nu]\)\)\/Gamma[2 - \[Nu]]\), "-", \(\(2\^\(\(-3\) + \[Nu]\)\ x\^\(3 - \[Nu]\)\ \[Beta]\^2\)\/Gamma[3 - \[Nu]]\), "+", InterpretationBox[\(O[\[Beta]]\^4\), SeriesData[ \[Beta], 0, {}, 0, 4, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 2, Plus[ -1, \[Nu]]], Power[ x, Plus[ 1, Times[ -1, \[Nu]]]], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1]], 0, Times[ -1, Power[ 2, Plus[ -3, \[Nu]]], Power[ x, Plus[ 3, Times[ -1, \[Nu]]]], Power[ Gamma[ Plus[ 3, Times[ -1, \[Nu]]]], -1]]}, 0, 4, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(CrD = Series[BesselJ[\[Nu] - 1, x*\[Beta]], {\[Beta], 0, 4}]\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(\(-1\) + \[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(2\^\(1 - \[Nu]\)\ x\^\(\(-1\) + \[Nu]\)\)\/Gamma[\[Nu]]\), "-", \(\(2\^\(\(-1\) - \[Nu]\)\ x\^\(1 + \[Nu]\)\ \[Beta]\^2\)\/Gamma[ 1 + \[Nu]]\), "+", InterpretationBox[\(O[\[Beta]]\^4\), SeriesData[ \[Beta], 0, {}, 0, 4, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 2, Plus[ 1, Times[ -1, \[Nu]]]], Power[ x, Plus[ -1, \[Nu]]], Power[ Gamma[ \[Nu]], -1]], 0, Times[ -1, Power[ 2, Plus[ -1, Times[ -1, \[Nu]]]], Power[ x, Plus[ 1, \[Nu]]], Power[ Gamma[ Plus[ 1, \[Nu]]], -1]]}, 0, 4, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Cra = FullSimplify[CrN/CrD] /. x -> a\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(\(-2\)\ \[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(4\^\(\(-1\) + \[Nu]\)\ a\^\(2 - 2\ \[Nu]\)\ Gamma[\[Nu]]\ \[Beta]\^2\)\/Gamma[2 - \[Nu]]\), "+", \(4\^\(\(-2\) + \[Nu]\)\ a\^\(4 - 2\ \[Nu]\)\ \((1\/\(\[Nu]\ Gamma[2 - \[Nu]]\) - 1\/Gamma[3 - \[Nu]])\)\ Gamma[\[Nu]]\ \[Beta]\^4\), "+", InterpretationBox[\(O[\[Beta]]\^6\), SeriesData[ \[Beta], 0, {}, 2, 6, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 4, Plus[ -1, \[Nu]]], Power[ a, Plus[ 2, Times[ -2, \[Nu]]]], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1], Gamma[ \[Nu]]], 0, Times[ Power[ 4, Plus[ -2, \[Nu]]], Power[ a, Plus[ 4, Times[ -2, \[Nu]]]], Plus[ Times[ Power[ \[Nu], -1], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1]], Times[ -1, Power[ Gamma[ Plus[ 3, Times[ -1, \[Nu]]]], -1]]], Gamma[ \[Nu]]]}, 2, 6, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Crb = FullSimplify[CrN/CrD] /. x -> b\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(\(-2\)\ \[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(4\^\(\(-1\) + \[Nu]\)\ b\^\(2 - 2\ \[Nu]\)\ Gamma[\[Nu]]\ \[Beta]\^2\)\/Gamma[2 - \[Nu]]\), "+", \(4\^\(\(-2\) + \[Nu]\)\ b\^\(4 - 2\ \[Nu]\)\ \((1\/\(\[Nu]\ Gamma[2 - \[Nu]]\) - 1\/Gamma[3 - \[Nu]])\)\ Gamma[\[Nu]]\ \[Beta]\^4\), "+", InterpretationBox[\(O[\[Beta]]\^6\), SeriesData[ \[Beta], 0, {}, 2, 6, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 4, Plus[ -1, \[Nu]]], Power[ b, Plus[ 2, Times[ -2, \[Nu]]]], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1], Gamma[ \[Nu]]], 0, Times[ Power[ 4, Plus[ -2, \[Nu]]], Power[ b, Plus[ 4, Times[ -2, \[Nu]]]], Plus[ Times[ Power[ \[Nu], -1], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1]], Times[ -1, Power[ Gamma[ Plus[ 3, Times[ -1, \[Nu]]]], -1]]], Gamma[ \[Nu]]]}, 2, 6, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Jn\[Beta]x = Series[BesselJ[\(-\[Nu]\), x*\[Beta]], {\[Beta], 0, 4}]\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(-\[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(2\^\[Nu]\ x\^\(-\[Nu]\)\)\/Gamma[1 - \[Nu]]\), "-", \(\(2\^\(\(-2\) + \[Nu]\)\ x\^\(2 - \[Nu]\)\ \[Beta]\^2\)\/Gamma[ 2 - \[Nu]]\), "+", InterpretationBox[\(O[\[Beta]]\^4\), SeriesData[ \[Beta], 0, {}, 0, 4, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 2, \[Nu]], Power[ x, Times[ -1, \[Nu]]], Power[ Gamma[ Plus[ 1, Times[ -1, \[Nu]]]], -1]], 0, Times[ -1, Power[ 2, Plus[ -2, \[Nu]]], Power[ x, Plus[ 2, Times[ -1, \[Nu]]]], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1]]}, 0, 4, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(J\[Beta]x = Series[BesselJ[\[Nu], x*\[Beta]], {\[Beta], 0, 4}]\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\[Nu]\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(2\^\(-\[Nu]\)\ x\^\[Nu]\)\/Gamma[1 + \[Nu]]\), "-", \(\(2\^\(\(-2\) - \[Nu]\)\ x\^\(2 + \[Nu]\)\ \[Beta]\^2\)\/Gamma[ 2 + \[Nu]]\), "+", InterpretationBox[\(O[\[Beta]]\^4\), SeriesData[ \[Beta], 0, {}, 0, 4, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 2, Times[ -1, \[Nu]]], Power[ x, \[Nu]], Power[ Gamma[ Plus[ 1, \[Nu]]], -1]], 0, Times[ -1, Power[ 2, Plus[ -2, Times[ -1, \[Nu]]]], Power[ x, Plus[ 2, \[Nu]]], Power[ Gamma[ Plus[ 2, \[Nu]]], -1]]}, 0, 4, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(DJ\[Beta]x = D[J\[Beta]x, x]/\[Beta]\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\[Nu]\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(2\^\(-\[Nu]\)\ x\^\(\(-1\) + \[Nu]\)\ \[Nu]\)\/\(Gamma[1 + \[Nu]]\ \[Beta]\)\), "-", \(\(2\^\(\(-2\) - \[Nu]\)\ x\^\(1 + \[Nu]\)\ \((2 + \[Nu])\)\ \[Beta]\)\/Gamma[2 + \[Nu]]\), "+", InterpretationBox[\(O[\[Beta]]\^3\), SeriesData[ \[Beta], 0, {}, -1, 3, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 2, Times[ -1, \[Nu]]], Power[ x, Plus[ -1, \[Nu]]], \[Nu], Power[ Gamma[ Plus[ 1, \[Nu]]], -1]], 0, Times[ -1, Power[ 2, Plus[ -2, Times[ -1, \[Nu]]]], Power[ x, Plus[ 1, \[Nu]]], Plus[ 2, \[Nu]], Power[ Gamma[ Plus[ 2, \[Nu]]], -1]]}, -1, 3, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(DJn\[Beta]x = D[Jn\[Beta]x, x]/\[Beta]\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(-\[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(-\(\(2\^\[Nu]\ x\^\(\(-1\) - \[Nu]\)\ \[Nu]\)\/\(Gamma[1 - \[Nu]]\ \[Beta]\)\)\), "-", \(\(2\^\(\(-2\) + \[Nu]\)\ x\^\(1 - \[Nu]\)\ \((2 - \[Nu])\)\ \[Beta]\)\/Gamma[2 - \[Nu]]\), "+", InterpretationBox[\(O[\[Beta]]\^3\), SeriesData[ \[Beta], 0, {}, -1, 3, 1]]}], SeriesData[ \[Beta], 0, { Times[ -1, Power[ 2, \[Nu]], Power[ x, Plus[ -1, Times[ -1, \[Nu]]]], \[Nu], Power[ Gamma[ Plus[ 1, Times[ -1, \[Nu]]]], -1]], 0, Times[ -1, Power[ 2, Plus[ -2, \[Nu]]], Power[ x, Plus[ 1, Times[ -1, \[Nu]]]], Plus[ 2, Times[ -1, \[Nu]]], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1]]}, -1, 3, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Jsum\[Beta]xSq = FullSimplify[\((Jn\[Beta]x + Cra*J\[Beta]x)\)^2]\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(\(-2\)\ \[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(4\^\[Nu]\ x\^\(\(-2\)\ \[Nu]\)\)\/Gamma[1 - \[Nu]]\^2\), "+", \(\(2\^\(\(-1\) + 2\ \[Nu]\)\ \((a\^\(2 - 2\ \[Nu]\) - x\^\(2 - 2\ \[Nu]\)\ \[Nu])\)\ Gamma[\[Nu]]\ Sin[\[Pi]\ \[Nu]]\ \[Beta]\^2\)\/\(\[Pi]\ \[Nu]\ Gamma[2 - \[Nu]]\)\), "+", InterpretationBox[\(O[\[Beta]]\^4\), SeriesData[ \[Beta], 0, {}, 0, 4, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 4, \[Nu]], Power[ x, Times[ -2, \[Nu]]], Power[ Gamma[ Plus[ 1, Times[ -1, \[Nu]]]], -2]], 0, Times[ Power[ 2, Plus[ -1, Times[ 2, \[Nu]]]], Power[ Pi, -1], Power[ \[Nu], -1], Plus[ Power[ a, Plus[ 2, Times[ -2, \[Nu]]]], Times[ -1, Power[ x, Plus[ 2, Times[ -2, \[Nu]]]], \[Nu]]], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1], Gamma[ \[Nu]], Sin[ Times[ Pi, \[Nu]]]]}, 0, 4, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(DJsum\[Beta]xSq = FullSimplify[\((DJn\[Beta]x + Cra*DJ\[Beta]x)\)^2]\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(\(-2\)\ \[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(4\^\[Nu]\ x\^\(\(-2\)\ \((1 + \[Nu])\)\)\)\/\(Gamma[\(-\[Nu]\)]\^2\ \[Beta]\^2\)\), "+", \(\(2\^\(\(-1\) + 2\ \[Nu]\)\ \((\(-\(a\^\(2 - 2\ \[Nu]\)\/x\^2\)\) - x\^\(\(-2\)\ \[Nu]\)\ \((\(-2\) + \[Nu])\))\)\ Gamma[1 + \[Nu]]\ Sin[\[Pi]\ \[Nu]]\)\/\(\[Pi]\ Gamma[2 - \[Nu]]\)\), "+", InterpretationBox[\(O[\[Beta]]\^2\), SeriesData[ \[Beta], 0, {}, -2, 2, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 4, \[Nu]], Power[ x, Times[ -2, Plus[ 1, \[Nu]]]], Power[ Gamma[ Times[ -1, \[Nu]]], -2]], 0, Times[ Power[ 2, Plus[ -1, Times[ 2, \[Nu]]]], Power[ Pi, -1], Plus[ Times[ -1, Power[ a, Plus[ 2, Times[ -2, \[Nu]]]], Power[ x, -2]], Times[ -1, Power[ x, Times[ -2, \[Nu]]], Plus[ -2, \[Nu]]]], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1], Gamma[ Plus[ 1, \[Nu]]], Sin[ Times[ Pi, \[Nu]]]]}, -2, 2, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Jsum\[Beta]bSq = Jsum\[Beta]xSq /. x \[Rule] b\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(\(-2\)\ \[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(4\^\[Nu]\ b\^\(\(-2\)\ \[Nu]\)\)\/Gamma[1 - \[Nu]]\^2\), "+", \(\(2\^\(\(-1\) + 2\ \[Nu]\)\ \((a\^\(2 - 2\ \[Nu]\) - b\^\(2 - 2\ \[Nu]\)\ \[Nu])\)\ Gamma[\[Nu]]\ Sin[\[Pi]\ \[Nu]]\ \[Beta]\^2\)\/\(\[Pi]\ \[Nu]\ Gamma[2 - \[Nu]]\)\), "+", InterpretationBox[\(O[\[Beta]]\^4\), SeriesData[ \[Beta], 0, {}, 0, 4, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 4, \[Nu]], Power[ b, Times[ -2, \[Nu]]], Power[ Gamma[ Plus[ 1, Times[ -1, \[Nu]]]], -2]], 0, Times[ Power[ 2, Plus[ -1, Times[ 2, \[Nu]]]], Power[ Pi, -1], Power[ \[Nu], -1], Plus[ Power[ a, Plus[ 2, Times[ -2, \[Nu]]]], Times[ -1, Power[ b, Plus[ 2, Times[ -2, \[Nu]]]], \[Nu]]], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1], Gamma[ \[Nu]], Sin[ Times[ Pi, \[Nu]]]]}, 0, 4, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Jsum\[Beta]aSq = Jsum\[Beta]xSq /. x \[Rule] a\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(\(-2\)\ \[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(4\^\[Nu]\ a\^\(\(-2\)\ \[Nu]\)\)\/Gamma[1 - \[Nu]]\^2\), "+", \(\(2\^\(\(-1\) + 2\ \[Nu]\)\ \((a\^\(2 - 2\ \[Nu]\) - a\^\(2 - 2\ \[Nu]\)\ \[Nu])\)\ Gamma[\[Nu]]\ Sin[\[Pi]\ \[Nu]]\ \[Beta]\^2\)\/\(\[Pi]\ \[Nu]\ Gamma[2 - \[Nu]]\)\), "+", InterpretationBox[\(O[\[Beta]]\^4\), SeriesData[ \[Beta], 0, {}, 0, 4, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 4, \[Nu]], Power[ a, Times[ -2, \[Nu]]], Power[ Gamma[ Plus[ 1, Times[ -1, \[Nu]]]], -2]], 0, Times[ Power[ 2, Plus[ -1, Times[ 2, \[Nu]]]], Power[ Pi, -1], Power[ \[Nu], -1], Plus[ Power[ a, Plus[ 2, Times[ -2, \[Nu]]]], Times[ -1, Power[ a, Plus[ 2, Times[ -2, \[Nu]]]], \[Nu]]], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1], Gamma[ \[Nu]], Sin[ Times[ Pi, \[Nu]]]]}, 0, 4, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(DJsum\[Beta]bSq = DJsum\[Beta]xSq /. x \[Rule] b\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(\(-2\)\ \[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(4\^\[Nu]\ b\^\(\(-2\)\ \((1 + \[Nu])\)\)\)\/\(Gamma[\(-\[Nu]\)]\^2\ \[Beta]\^2\)\), "+", \(\(2\^\(\(-1\) + 2\ \[Nu]\)\ \((\(-\(a\^\(2 - 2\ \[Nu]\)\/b\^2\)\) - b\^\(\(-2\)\ \[Nu]\)\ \((\(-2\) + \[Nu])\))\)\ Gamma[1 + \[Nu]]\ Sin[\[Pi]\ \[Nu]]\)\/\(\[Pi]\ Gamma[2 - \[Nu]]\)\), "+", InterpretationBox[\(O[\[Beta]]\^2\), SeriesData[ \[Beta], 0, {}, -2, 2, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 4, \[Nu]], Power[ b, Times[ -2, Plus[ 1, \[Nu]]]], Power[ Gamma[ Times[ -1, \[Nu]]], -2]], 0, Times[ Power[ 2, Plus[ -1, Times[ 2, \[Nu]]]], Power[ Pi, -1], Plus[ Times[ -1, Power[ a, Plus[ 2, Times[ -2, \[Nu]]]], Power[ b, -2]], Times[ -1, Power[ b, Times[ -2, \[Nu]]], Plus[ -2, \[Nu]]]], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1], Gamma[ Plus[ 1, \[Nu]]], Sin[ Times[ Pi, \[Nu]]]]}, -2, 2, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(DJsum\[Beta]aSq = DJsum\[Beta]xSq /. x \[Rule] a\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(\(-2\)\ \[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(4\^\[Nu]\ a\^\(\(-2\)\ \((1 + \[Nu])\)\)\)\/\(Gamma[\(-\[Nu]\)]\^2\ \[Beta]\^2\)\), "+", \(\(2\^\(\(-1\) + 2\ \[Nu]\)\ \((\(-a\^\(\(-2\)\ \[Nu]\)\) - a\^\(\(-2\)\ \[Nu]\)\ \((\(-2\) + \[Nu])\))\)\ Gamma[1 + \[Nu]]\ Sin[\[Pi]\ \[Nu]]\)\/\(\[Pi]\ Gamma[2 - \[Nu]]\)\), "+", InterpretationBox[\(O[\[Beta]]\^2\), SeriesData[ \[Beta], 0, {}, -2, 2, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 4, \[Nu]], Power[ a, Times[ -2, Plus[ 1, \[Nu]]]], Power[ Gamma[ Times[ -1, \[Nu]]], -2]], 0, Times[ Power[ 2, Plus[ -1, Times[ 2, \[Nu]]]], Power[ Pi, -1], Plus[ Times[ -1, Power[ a, Times[ -2, \[Nu]]]], Times[ -1, Power[ a, Times[ -2, \[Nu]]], Plus[ -2, \[Nu]]]], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1], Gamma[ Plus[ 1, \[Nu]]], Sin[ Times[ Pi, \[Nu]]]]}, -2, 2, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Norm\[Beta]b = FullSimplify[ \((b^2)\)/2* \((DJsum\[Beta]bSq + \((1 - \((\[Nu]/\((b*\[Beta])\))\)^2)\)*Jsum\[Beta]bSq)\)]\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(\(-2\)\ \[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(2\^\(\(-1\) + 2\ \[Nu]\)\ \((\((\(-a\^\(2 - 2\ \[Nu]\)\) + b\^\(2 - 2\ \[Nu]\))\)\ \[Nu]\ Gamma[1 - \[Nu]] + b\^\(2 - 2\ \[Nu]\)\ Gamma[2 - \[Nu]]) \)\)\/\(Gamma[1 - \[Nu]]\^2\ Gamma[2 - \[Nu]]\)\), "+", InterpretationBox[\(O[\[Beta]]\^2\), SeriesData[ \[Beta], 0, {}, 0, 2, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 2, Plus[ -1, Times[ 2, \[Nu]]]], Power[ Gamma[ Plus[ 1, Times[ -1, \[Nu]]]], -2], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1], Plus[ Times[ Plus[ Times[ -1, Power[ a, Plus[ 2, Times[ -2, \[Nu]]]]], Power[ b, Plus[ 2, Times[ -2, \[Nu]]]]], \[Nu], Gamma[ Plus[ 1, Times[ -1, \[Nu]]]]], Times[ Power[ b, Plus[ 2, Times[ -2, \[Nu]]]], Gamma[ Plus[ 2, Times[ -1, \[Nu]]]]]]]}, 0, 2, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Norm\[Beta]a = FullSimplify[ \((a^2)\)/2* \((DJsum\[Beta]aSq + \((1 - \((\[Nu]/\((a*\[Beta])\))\)^2)\)*Jsum\[Beta]aSq)\)]\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(\(-2\)\ \[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(2\^\(\(-1\) + 2\ \[Nu]\)\ a\^\(2 - 2\ \[Nu]\)\)\/Gamma[1 - \[Nu]]\^2\), "+", InterpretationBox[\(O[\[Beta]]\^2\), SeriesData[ \[Beta], 0, {}, 0, 2, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 2, Plus[ -1, Times[ 2, \[Nu]]]], Power[ a, Plus[ 2, Times[ -2, \[Nu]]]], Power[ Gamma[ Plus[ 1, Times[ -1, \[Nu]]]], -2]]}, 0, 2, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Cintegrand1 = FullSimplify[\((x^\[Nu])\)*\((Jn\[Beta]x + Crb*J\[Beta]x)\)]\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\(-\[Nu]\)\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(2\^\[Nu]\/Gamma[1 - \[Nu]]\), "+", \(\(2\^\(\(-2\) + \[Nu]\)\ \((b\^\(2 - 2\ \[Nu]\)\ x\^\(2\ \[Nu]\)\ Gamma[\[Nu]] - x\^2\ Gamma[1 + \[Nu]])\)\ \[Beta]\^2\)\/\(Gamma[2 - \[Nu]]\ Gamma[1 + \[Nu]]\)\), "+", InterpretationBox[\(O[\[Beta]]\^4\), SeriesData[ \[Beta], 0, {}, 0, 4, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 2, \[Nu]], Power[ Gamma[ Plus[ 1, Times[ -1, \[Nu]]]], -1]], 0, Times[ Power[ 2, Plus[ -2, \[Nu]]], Power[ Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], -1], Power[ Gamma[ Plus[ 1, \[Nu]]], -1], Plus[ Times[ Power[ b, Plus[ 2, Times[ -2, \[Nu]]]], Power[ x, Times[ 2, \[Nu]]], Gamma[ \[Nu]]], Times[ -1, Power[ x, 2], Gamma[ Plus[ 1, \[Nu]]]]]]}, 0, 4, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Cintgnd1norm = FullSimplify[Cintegrand1/\((Norm\[Beta]b - Norm\[Beta]a)\)]\)], "Input"], Cell[BoxData[ RowBox[{\(\[Beta]\^\[Nu]\), " ", RowBox[{"(", InterpretationBox[ RowBox[{ \(\(2\^\(1 - \[Nu]\)\ a\^\(2\ \[Nu]\)\ b\^\(2\ \[Nu]\)\ Gamma[1 - \[Nu]]\ Gamma[2 - \[Nu]]\)\/\(\(( a\^\(2\ \[Nu]\)\ b\^2 - a\^2\ b\^\(2\ \[Nu]\))\)\ \((\[Nu]\ Gamma[1 - \[Nu]] + Gamma[2 - \[Nu]])\)\)\), "+", InterpretationBox[\(O[\[Beta]]\^2\), SeriesData[ \[Beta], 0, {}, 0, 2, 1]]}], SeriesData[ \[Beta], 0, { Times[ Power[ 2, Plus[ 1, Times[ -1, \[Nu]]]], Power[ a, Times[ 2, \[Nu]]], Power[ b, Times[ 2, \[Nu]]], Power[ Plus[ Times[ Power[ a, Times[ 2, \[Nu]]], Power[ b, 2]], Times[ -1, Power[ a, 2], Power[ b, Times[ 2, \[Nu]]]]], -1], Gamma[ Plus[ 1, Times[ -1, \[Nu]]]], Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], Power[ Plus[ Times[ \[Nu], Gamma[ Plus[ 1, Times[ -1, \[Nu]]]]], Gamma[ Plus[ 2, Times[ -1, \[Nu]]]]], -1]]}, 0, 2, 1]], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Cintgnd1norm /. \[Beta] \[Rule] 0\)], "Input"], Cell[BoxData[ \(\(0\^\[Nu]\ 2\^\(1 - \[Nu]\)\ a\^\(2\ \[Nu]\)\ b\^\(2\ \[Nu]\)\ Gamma[1 - \[Nu]]\ Gamma[2 - \[Nu]]\)\/\(\(( a\^\(2\ \[Nu]\)\ b\^2 - a\^2\ b\^\(2\ \[Nu]\))\)\ \((\[Nu]\ Gamma[1 - \[Nu]] + Gamma[2 - \[Nu]])\)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TermN0 = FullSimplify[Cintgnd1norm*\((x^\[Alpha])\)*Jn\[Beta]x]\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{ \(\(2\ a\^\(2\ \[Nu]\)\ b\^\(2\ \[Nu]\)\ x\^\(\[Alpha] - \[Nu]\)\ Gamma[2 - \[Nu]]\)\/\(\(( a\^\(2\ \[Nu]\)\ b\^2 - a\^2\ b\^\(2\ \[Nu]\))\)\ \((\[Nu]\ Gamma[1 - \[Nu]] + Gamma[2 - \[Nu]])\)\)\), "+", InterpretationBox[\(O[\[Beta]]\^2\), SeriesData[ \[Beta], 0, {}, 0, 2, 1]]}], SeriesData[ \[Beta], 0, { Times[ 2, Power[ a, Times[ 2, \[Nu]]], Power[ b, Times[ 2, \[Nu]]], Power[ Plus[ Times[ Power[ a, Times[ 2, \[Nu]]], Power[ b, 2]], Times[ -1, Power[ a, 2], Power[ b, Times[ 2, \[Nu]]]]], -1], Power[ x, Plus[ \[Alpha], Times[ -1, \[Nu]]]], Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], Power[ Plus[ Times[ \[Nu], Gamma[ Plus[ 1, Times[ -1, \[Nu]]]]], Gamma[ Plus[ 2, Times[ -1, \[Nu]]]]], -1]]}, 0, 2, 1]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Term00 = Cintgnd1norm*Crb*\((x^\[Alpha])\)*J\[Beta]x\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{ \(\(a\^\(2\ \[Nu]\)\ b\^2\ x\^\(\[Alpha] + \[Nu]\)\ Gamma[1 - \[Nu]]\ Gamma[\[Nu]]\ \[Beta]\^2\)\/\(2\ \((a\^\(2\ \[Nu]\)\ b\^2 - a\^2\ b\^\(2\ \[Nu]\))\)\ \((\[Nu]\ Gamma[1 - \[Nu]] + Gamma[2 - \[Nu]])\)\ Gamma[1 + \[Nu]]\)\), "+", InterpretationBox[\(O[\[Beta]]\^4\), SeriesData[ \[Beta], 0, {}, 2, 4, 1]]}], SeriesData[ \[Beta], 0, { Times[ Rational[ 1, 2], Power[ a, Times[ 2, \[Nu]]], Power[ b, 2], Power[ Plus[ Times[ Power[ a, Times[ 2, \[Nu]]], Power[ b, 2]], Times[ -1, Power[ a, 2], Power[ b, Times[ 2, \[Nu]]]]], -1], Power[ x, Plus[ \[Alpha], \[Nu]]], Gamma[ Plus[ 1, Times[ -1, \[Nu]]]], Power[ Plus[ Times[ \[Nu], Gamma[ Plus[ 1, Times[ -1, \[Nu]]]]], Gamma[ Plus[ 2, Times[ -1, \[Nu]]]]], -1], Gamma[ \[Nu]], Power[ Gamma[ Plus[ 1, \[Nu]]], -1]]}, 2, 4, 1]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Term0 = FullSimplify[Term00]\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{ \(\(a\^\(2\ \[Nu]\)\ b\^2\ x\^\(\[Alpha] + \[Nu]\)\ Gamma[\(-\[Nu]\)]\ \[Beta]\^2\)\/\(2\ \((\(-a\^\(2\ \[Nu]\)\)\ b\^2 + a\^2\ b\^\(2\ \[Nu]\))\)\ \((\[Nu]\ Gamma[1 - \[Nu]] + Gamma[2 - \[Nu]])\)\)\), "+", InterpretationBox[\(O[\[Beta]]\^4\), SeriesData[ \[Beta], 0, {}, 2, 4, 1]]}], SeriesData[ \[Beta], 0, { Times[ Rational[ 1, 2], Power[ a, Times[ 2, \[Nu]]], Power[ b, 2], Power[ Plus[ Times[ -1, Power[ a, Times[ 2, \[Nu]]], Power[ b, 2]], Times[ Power[ a, 2], Power[ b, Times[ 2, \[Nu]]]]], -1], Power[ x, Plus[ \[Alpha], \[Nu]]], Power[ Plus[ Times[ \[Nu], Gamma[ Plus[ 1, Times[ -1, \[Nu]]]]], Gamma[ Plus[ 2, Times[ -1, \[Nu]]]]], -1], Gamma[ Times[ -1, \[Nu]]]]}, 2, 4, 1]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TermT = FullSimplify[Term0 + TermN0]\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{ \(\(2\ a\^\(2\ \[Nu]\)\ b\^\(2\ \[Nu]\)\ x\^\(\[Alpha] - \[Nu]\)\ Gamma[2 - \[Nu]]\)\/\(\(( a\^\(2\ \[Nu]\)\ b\^2 - a\^2\ b\^\(2\ \[Nu]\))\)\ \((\[Nu]\ Gamma[1 - \[Nu]] + Gamma[2 - \[Nu]])\)\)\), "+", InterpretationBox[\(O[\[Beta]]\^2\), SeriesData[ \[Beta], 0, {}, 0, 2, 1]]}], SeriesData[ \[Beta], 0, { Times[ 2, Power[ a, Times[ 2, \[Nu]]], Power[ b, Times[ 2, \[Nu]]], Power[ Plus[ Times[ Power[ a, Times[ 2, \[Nu]]], Power[ b, 2]], Times[ -1, Power[ a, 2], Power[ b, Times[ 2, \[Nu]]]]], -1], Power[ x, Plus[ \[Alpha], Times[ -1, \[Nu]]]], Gamma[ Plus[ 2, Times[ -1, \[Nu]]]], Power[ Plus[ Times[ \[Nu], Gamma[ Plus[ 1, Times[ -1, \[Nu]]]]], Gamma[ Plus[ 2, Times[ -1, \[Nu]]]]], -1]]}, 0, 2, 1]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(FullSimplify[TermT /. \[Beta] \[Rule] 0]\)], "Input"], Cell[BoxData[ \(\(2\ a\^\(2\ \[Nu]\)\ b\^\(2\ \[Nu]\)\ x\^\(\[Alpha] - \[Nu]\)\ Gamma[2 - \[Nu]]\)\/\(\(( a\^\(2\ \[Nu]\)\ b\^2 - a\^2\ b\^\(2\ \[Nu]\))\)\ \((\[Nu]\ Gamma[1 - \[Nu]] + Gamma[2 - \[Nu]])\)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TermT1 = FullSimplify[\ Gamma[2 - \[Nu]]/\((\ Gamma[2 - \[Nu]]\ + \[Nu]*Gamma[1 - \[Nu]])\)] \)], "Input"], Cell[BoxData[ \(Gamma[2 - \[Nu]]\/\(\[Nu]\ Gamma[1 - \[Nu]] + Gamma[2 - \[Nu]]\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TermT1 = FullSimplify[\ 1/\((\ 1 + \((\[Nu]*Gamma[1 - \[Nu]])\)/Gamma[2 - \[Nu]])\)]\)], "Input"], Cell[BoxData[ \(1 - \[Nu]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TermT2 = FullSimplify[ \((a\^\(2\ \[Nu]\)\ b\^\(2\ \[Nu]\))\)/ \((a\^\(2\ \[Nu]\)\ b\^2 - a\^2\ b\^\(2\ \[Nu]\))\)]\)], "Input"], Cell[BoxData[ \(1\/\(\(-a\^\(2 - 2\ \[Nu]\)\) + b\^\(2 - 2\ \[Nu]\)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TermT = TermT1*TermT2*2*x^\((\[Alpha] - \[Nu])\)\)], "Input"], Cell[BoxData[ \(\(2\ x\^\(\[Alpha] - \[Nu]\)\ \((1 - \[Nu])\)\)\/\(\(-a\^\(2 - 2\ \[Nu]\)\) + b\^\(2 - 2\ \[Nu]\)\)\)], "Output"] }, Open ]] }, FrontEndVersion->"Microsoft Windows 3.0", ScreenRectangle->{{0, 1024}, {0, 712}}, WindowSize->{747, 584}, WindowMargins->{{2, Automatic}, {Automatic, 5}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic} ] (*********************************************************************** 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[1731, 51, 93, 1, 30, "Input"], Cell[1827, 54, 1090, 30, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2954, 89, 93, 1, 30, "Input"], Cell[3050, 92, 996, 27, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4083, 124, 70, 1, 30, "Input"], Cell[4156, 127, 1531, 42, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5724, 174, 70, 1, 30, "Input"], Cell[5797, 177, 1531, 42, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7365, 224, 108, 2, 30, "Input"], Cell[7476, 228, 1014, 28, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8527, 261, 98, 2, 30, "Input"], Cell[8628, 265, 943, 26, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9608, 296, 69, 1, 30, "Input"], Cell[9680, 299, 1078, 29, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10795, 333, 71, 1, 30, "Input"], Cell[10869, 336, 1186, 32, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12092, 373, 101, 2, 30, "Input"], Cell[12196, 377, 1506, 41, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13739, 423, 108, 2, 30, "Input"], Cell[13850, 427, 1669, 46, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15556, 478, 79, 1, 30, "Input"], Cell[15638, 481, 1506, 41, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17181, 527, 79, 1, 30, "Input"], Cell[17263, 530, 1506, 41, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18806, 576, 81, 1, 30, "Input"], Cell[18890, 579, 1669, 46, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[20596, 630, 81, 1, 30, "Input"], Cell[20680, 633, 1593, 44, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[22310, 682, 197, 6, 30, "Input"], Cell[22510, 690, 1625, 44, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[24172, 739, 197, 6, 30, "Input"], Cell[24372, 747, 735, 20, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[25144, 772, 114, 2, 30, "Input"], Cell[25261, 776, 1448, 41, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[26746, 822, 114, 2, 30, "Input"], Cell[26863, 826, 1614, 46, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[28514, 877, 66, 1, 30, "Input"], Cell[28583, 880, 279, 5, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[28899, 890, 99, 2, 30, "Input"], Cell[29001, 894, 1283, 39, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[30321, 938, 85, 1, 30, "Input"], Cell[30409, 941, 1401, 43, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[31847, 989, 61, 1, 30, "Input"], Cell[31911, 992, 1233, 37, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[33181, 1034, 69, 1, 30, "Input"], Cell[33253, 1037, 1283, 39, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[34573, 1081, 73, 1, 30, "Input"], Cell[34649, 1084, 252, 4, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[34938, 1093, 151, 4, 30, "Input"], Cell[35092, 1099, 101, 2, 44, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[35230, 1106, 136, 4, 30, "Input"], Cell[35369, 1112, 43, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[35449, 1118, 173, 4, 31, "Input"], Cell[35625, 1124, 87, 1, 42, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[35749, 1130, 81, 1, 30, "Input"], Cell[35833, 1133, 151, 3, 43, "Output"] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)