Series expansion from the red book on special functions by Richard Ask

Click For Summary
SUMMARY

The discussion focuses on verifying a series expansion from Richard Ask's book on special functions using Mathematica. The correct expansion for the expression $$\frac{(1+\frac{1}{j})^x}{1+x/j}$$ is confirmed to be $$1+\frac{x(x-1)}{2j^2}+\mathcal{O}(1/j^3)$$. A participant suggests using the Mathematica command Series to obtain the series expansion about $$k=0$$, where $$k \equiv 1/j$$. The incorrect term $$-\frac{x(x+1)}{2j^2}$$ arises from omitting the $$x^2/j^2$$ term during multiplication of the expansions.

PREREQUISITES
  • Understanding of series expansions in calculus
  • Familiarity with Mathematica software
  • Knowledge of asymptotic notation, specifically $$\mathcal{O}$$ notation
  • Basic algebraic manipulation of series
NEXT STEPS
  • Learn how to use the Series command in Mathematica for Taylor expansions
  • Study the properties of asymptotic expansions in mathematical analysis
  • Explore the concept of higher-order terms in series expansions
  • Investigate the implications of neglecting terms in series multiplication
USEFUL FOR

Mathematicians, physicists, and students involved in advanced calculus or mathematical analysis, particularly those working with series expansions and Mathematica.

MathematicalPhysicist
Science Advisor
Gold Member
Messages
4,662
Reaction score
372
I want to check my calculations via mathematica.

In the book I am reading there's this expansion:
$$\frac{(1+\frac{1}{j})^x}{1+x/j}=1+\frac{x(x-1)}{2j^2}+\mathcal{O}(1/j^3)$$

though I get instead of the term ##\frac{x(x-1)}{2j^2}## in the rhs the term: ##-\frac{x(x+1)}{2j^2}##.

So I want to check by mathematica if my calculations are correct, how do you suggest me to implement it in mathematica?

Thanks!
 
Last edited by a moderator:
Physics news on Phys.org
I don't have experience with Mathematica, but I suggest you try the example x=2 to see if your expansion is right.
 
  • Like
Likes   Reactions: DrClaude and FactChecker
Define ##k \equiv 1/j## then use the command Series to get the series expansion about ##k=0##.
The book is correct.
 
  • Like
Likes   Reactions: docnet and FactChecker
DrClaude said:
Define ##k \equiv 1/j## then use the command Series to get the series expansion about ##k=0##.
The book is correct.
Yeah I know.
I thought to myself to expand ##(1+1/j)^x = 1+x/j+x(x-1)/(2j^2)+\mathcal{O}(1/j^3)## and ##(1+x/j)^{-1}=1-x/j+\mathcal{O}(1/j^2)## and then multiply both factors; but it seems not to be the correct approach.
 
MathematicalPhysicist said:
Yeah I know.
I thought to myself to expand ##(1+1/j)^x = 1+x/j+x(x-1)/(2j^2)+\mathcal{O}(1/j^3)## and ##(1+x/j)^{-1}=1-x/j+\mathcal{O}(1/j^2)## and then multiply both factors; but it seems not to be the correct approach.

You need to expand both to order j^{-2}: <br /> \begin{align*}<br /> \left(1 + \frac 1j\right)^x &amp;= 1 + \frac xj + \frac{x(x-1)}{2j^2} + O(j^{-3}) \\<br /> \left(1 + \frac xj \right)^{-1} &amp;= 1 - \frac xj + \frac{x^2}{j^2} + O(j^{-3}) \end{align*}. Now multiply these together, ignoring any term of order j^{-3} or higher: \begin{split}<br /> \left(1 - \frac xj + \frac{x^2}{j^2} + O(j^{-3})\right) + \frac xj \left( 1 - \frac xj + O(j^{-2}) \right) + <br /> \frac{x(x-1)}{2j^2} \left(1 + O(j^{-1})\right) \\= 1 + \frac{x(x-1)}{2j^2} + O(j^{-3}).\end{split} Omitting +x^2/j^2 from the first bracket will give the incorrect <br /> 1 - \frac{x(x+1)}{2j^2} + O(j^{-3}) from your original post.
 
  • Like
Likes   Reactions: docnet and DrClaude

Similar threads

  • · Replies 13 ·
Replies
13
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
4K
Replies
3
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 3 ·
Replies
3
Views
1K
Replies
1
Views
2K
  • · Replies 12 ·
Replies
12
Views
2K