Questions on dealing with series data and constructing power series expansions of functions.

learn more… | top users | synonyms (1)

12
votes
4answers
790 views

Evaluation of a triple sum does not finish in reasonable time

I'm trying to compute the following triple sum, but no result is produced within a reasonable amount of time. What to do? ...
10
votes
3answers
2k views

Multivariable Taylor expansion does not work as expected

The basic multivariable Taylor expansion formula around a point is as follows: $$ f(\mathbf r + \mathbf a) = f(\mathbf r) + (\mathbf a \cdot \nabla )f(\mathbf r) + \frac{1}{2!}(\mathbf a \cdot ...
10
votes
1answer
263 views

Why does Mathematica fail to series expand this simple expression?

I wanted to expand the function $(x+2)^{x+2}$ around $x = -1$, that is, using Series[(x + 2)^(x + 2), {x, -1, 2}] and Mathematica returns the same expression. ...
9
votes
2answers
509 views

Find closed form expression for series expansion coefficients [duplicate]

Is there a built-in function that will find a general expression for the coefficient of the series expansion of a function? Series will only give the explicit ...
9
votes
2answers
570 views

Series expansion of an inverse

I have to find the series expansion of the inverse function of : $\arctan\left(\frac{\ln(1+x)}{1+x}\right)$ How do I find out the series expansion of any inverse ? Note: The inverse of a function ...
9
votes
1answer
537 views

How does Mathematica understand branchcuts of the complex logarithm?

Say I have the function $f(x) = x \tanh(\pi x) \log (x^2 +a^2)$ where $a$ is some positive real number. Then it seems to be me that Mathematica when given such a ...
9
votes
2answers
365 views

InverseSeries of multiple variables and multiple equations

CONTEXT Let us consider a bit of the Universe in which we draw spheres (see a high resolution image here). Astronomers have shown that the density within these spheres could be predicted quite ...
9
votes
2answers
143 views

Unexpected behaviour of SeriesCoefficient?

Following this question/answer I discovered and played with SeriesCoefficient. In particular, I tried ...
9
votes
0answers
198 views

Why does $\frac{\partial}{\partial x}O\left(\left(\frac{1}{x}\right)^0\right)$ equal $O\left(\left(\frac{1}{x}\right)^0\right)$ in a series expansion?

When taking the derivative of a series expansion around a finite point, the $O(x^n)$ part is differentiated as expected. $O(x^n)$ becomes $O(x^{n-1})$ except $O(x^0)$ which stays $O(x^0)$. When ...
8
votes
3answers
223 views

Declaration of abstract matrices to perform series expansion on them

I would like to have abstract matrices M and S to get out the coefficients of matrix power series however it treats M and S as numbers even if i checked that M.S - S.M != 0. I attach my code below: ...
8
votes
2answers
349 views

Series expansion with irrational power

I need the series expansion of a fairly nasty function and its derivative: ...
8
votes
1answer
1k views

Solving an ODE in power series

How do I find a series solution to an ODE? I do not mean taking the Taylor series of an exact solution; I want to solve nasty nonlinear differential equations locally via plug and chug. Surely, that ...
8
votes
1answer
189 views

Why is this infinite series wrongly computed by Mathematica?

Could you let me know if Mathematica (newer versions) is able to correctly compute this one? Sum[(-1)^(n + 1) Cos[3^n x]^3/3^n, {n, 1, Infinity}]
8
votes
0answers
207 views

How does Mathematica find a series expansion of expressions containing logarithms when there is a singularity at the expansion point?

I am looking for a good approximation to a function containing logarithms, especially at values close to zero. When I used Mathematica's Series command I found ...
7
votes
2answers
1k views

How to study asymptotic behavior, built-in functions

My question is as follows. Suppose we have a function $f(r)$ and we want to study its asymptotic behavior at infinity ($r\rightarrow \infty$). For example, the function may reduce to $-\frac{a}{r}$ or ...
7
votes
1answer
224 views

Series expansion for small ratio of variables

I have a messy expression of variables that I would like to simplify under the assumption that certain ratios of the variables are small. For example, consider $\sqrt{x+y}\;$ expanded for small ...
6
votes
3answers
335 views

How to make all numbers equal to one in a Series?

I have many outputs of one-dimensional Series expansions for which I am only interested in the general tending with the variable. For example, I would like to be able to transform something like ...
6
votes
1answer
747 views

How to expand a function into a power series with negative powers?

Is there any way to expand this expression a+b(1-Exp[-T/(b c)]/(z-Exp[-T/ (b c)]) (where a, ...
6
votes
1answer
179 views

Sophistication of Series[…]

I'll give a concrete example and I hope that my general question will be clear. Say I have three variables, $f$, $g$, $h$, and I know that $f=\mathcal O(x)$, $g=\mathcal O(x^2)$, $h=\mathcal O(x^3)$ ...
6
votes
1answer
76 views

Is it possible to find a function from first few terms in the expansion

Is it possible to find a function if first few terms of the expansion is known. For example if I have this series $f(x)=\frac{k^3 x^2}{6}-\frac{k^5 x^4}{120}+\frac{k^7 x^6}{5040}-\frac{k^9 ...
6
votes
1answer
111 views

Is it possible to circumvent a bug inside SeriesCoefficient?

As far as I can tell, there seems to be a bug in SeriesCoefficient: ...
6
votes
2answers
455 views

Solution of equation with power series (perturbation)

So I need to use Mathematica to find the solution of $y=x- \epsilon \sin(2y)$ as a power series in terms of $\epsilon$. I'd assume I'd need to create an equation $f=x-y- \epsilon \sin(2y)$, then ...
6
votes
0answers
108 views

Design considerations behind `O` (a.k.a. BigOh, a.k.a. Landau Order)

This works without any warnings: O[Log[x]]. This raises a warning: O[x^2]. I have a few questions around this: Why is it a ...
5
votes
5answers
1k views

Series expansion in terms of Hermite polynomials

I am trying to expand a polynomial in terms of orthogonal polynomials (in my case, Hermite). Maple has a nice built-in function for this, ChangeBasis. Is there a ...
5
votes
2answers
154 views

How to get the series expansion of $e^{x^a}$ at $x=0$?

I want to have a series expansion of $e^{x^a}$ or $e^{c_1x^a+c_2x^b}$ at $x=0$, but Series cannot give any useful result even if the assumption $a>0$ is ...
5
votes
1answer
156 views

How to compute the residue of $e^{z-\frac{1}{z}}$ at z=0?

I've tried the following but it didn't work: Residue[Exp[z - 1/z], {z, 0}] not even this: Residue[Exp[1/z], {z, 0}] ...
5
votes
1answer
112 views

Apparently contradictory output from Series and D

I defined a function as: g[b_] := Integrate[f[n]Exp[-b DA (n.u)^2], n] Obviously D[g[b], b] /. b -> 0 gives ...
5
votes
0answers
99 views

Expansion of $E(i c \mid m)$ at $c\to\infty$?

Currently, I am using a Windows machine with Mathematica 8. I noticed a difference in a series expansion of the function EllipticE[] in comparison with a result ...
4
votes
2answers
473 views

How do I solve N simultaneous equations for N variables?

I have a function: f[x_] := x + 31 x^3 + 5 x^25 Which I want to find an expansion for: ...
4
votes
3answers
105 views

Get interval of series

It seems to be a stupid question but I wonder how to get an interval of a series expansion. The current series command Series[f, {x, x0, n}] only give series ...
4
votes
1answer
199 views

Asymptotic expansion

I wanted to expand a function of $x$ about $x=\infty$ and see its coefficients as a function of the parameters $m,n,q,y$ and I wrote this - but it didn't work! It gave me back a very complicated ...
4
votes
1answer
199 views

Computing a series in terms of exponential function

Is there any way to compute the following series in terms of exponential function ? $$\sum_{k=0}^\infty Y_1(k)\;x^k$$ where $$Y_1(k) = \frac{(k - 1)!}{k!}Y_3(k - 1)$$ $$Y_2(k) = \frac{(k - ...
4
votes
1answer
156 views

List manipulation - various

Update OK, Please forgive the messiness of this, but I am working with something like: ...
4
votes
1answer
244 views

Erfi[z] expansion in Mathematica: is this a bug?

Looking at the expansion: Series[Erfi[x], {x, Infinity, 1}] I obtain -I+E^x^2 (1/(Sqrt[\[Pi]] x)+O[1/x]^2) (note the ...
4
votes
1answer
86 views

Series expansion for $\frac{x}{1- \frac{1}{x}}$

I would like to expand $\frac{x}{1- \frac{1}{x}}$ as $$\frac{x}{1- \frac{1}{x}} = x \left( 1+ \frac{1}{x} + \frac{1}{x^2} +\frac{1}{x^3} + \cdots \right) = x + 1 + \frac{1}{x} + \frac{1}{x^2} + ...
4
votes
0answers
108 views

Proper treatment of roots and powers in Series?

I have the following problem in Mathematica 9 on Linux. I let Mathematica compute the Series expansion: ...
3
votes
3answers
436 views

Limiting form of a polynomial expression

When simplifying an expression by hand, a trick that is often used is to remove terms that are lower powers of the independent variable, for instance, as $x \rightarrow \infty$, $x^2 + x$ becomes ...
3
votes
3answers
183 views

What is the formula for this numerical series?

I'm developing a questions game. My goal is that the score for each correct answer will increase as the user answers more questions. Initially there are 15 points for each correct answer. Every 4 ...
3
votes
4answers
153 views

Exponential of a Differential Operator

In Mathematica, is it possible to exponentiate a differential operator such that the operator will act on a function, $f(x,p)$? Specifically, I wondering if I can get Mathematica to do this: ...
3
votes
2answers
105 views

Zero order term in an exprerssion

I have a very long expression like a + b x/y + f + c y/x + d (x z)/y + ... where a , b, c, f ... are coefficients and ...
3
votes
1answer
245 views

How to get rid of all complex numbers and functions?

For almost esthetical reasons, I always have been fascinated by infinite sums of series and I always wonder if theyhave a closed form. A couple of days ago, I found on MSE a question related to the ...
3
votes
2answers
119 views

Change all values of a series to positive values

I would like to change all values of a series to positive values. I refer to my previous question, answered in the most part, very elegantly by RunnyKine. However, there remains a slight glitch that ...
3
votes
3answers
166 views

Series coefficients for an infinite sum?

I'm working through Carl Bender's Mathematical Physics lectures on YouTube (which are great fun), and I'd like Mathematica's help solving terms in the perturbation series. It would be convenient if ...
3
votes
1answer
140 views

Problems with series of generalized hypergeometric functions

I have been trying to do the following series: Series[HypergeometricPFQ[{1, 4, 4}, {4 - Sqrt[3], 4 + Sqrt[3]}, z],{z,1,0}] Mathematica says that the result is ...
2
votes
1answer
131 views

Is there a function that, given a fraction, will return the general term of its infinite series expansion?

Is here some way to expand a fraction to an infinite sum in mathematica, i.e., a series? I want the general term of the series. For example, $\frac{2}{3(x-1)^3}$
2
votes
1answer
99 views

Remove singularity at zero from $-((i^{-n} (-1 + i^n)^2)/(n^2 π^2))$

The expression h[n_] = -((I^-n (-1 + I^n)^2)/(n^2 π^2)) is real for all integers n. Although indeterminate at ...
2
votes
1answer
935 views

Laurent series expansion

Can someone share how to find a Laurent series expansion of $$f(z)=\frac{1}{(z^2-1)(z^2-4)}$$ centered at zero on the annular disk $1<|z|<2$?
2
votes
1answer
114 views

Matrix exponential MatrixExp[] vs Sum[MatrixPower[]] doesn't match?

I might be an idiot, but I cannot get the manual expansion of e^At to match the MatrixExp[A t] result. For example, I have the following: ...
2
votes
2answers
102 views

Bessel derivatives

I tried to expand BesselJ[k,x] function into a Taylor series with Series command. Here both ...
2
votes
1answer
76 views

Taylor expansion of function of vectors—simplifying the output

I've got a real-valued function of several vectors $f(u,v,w)$ formed by taking scalar products of linear combinations of the vectors, I want to Taylor expand around small $v$ by writing $$f(u,\delta ...