Questions using the Sum command, especially for series and other algebraic objects, and related functions such as SumConvergence

learn more… | top users | synonyms

11
votes
2answers
153 views

How to correctly implement in a new function the scoping behavior of Table, Sum and other commands that use Block to localize iterators?

It is documented that "Block is automatically used to localize values of iterators in iteration constructs such as Do, Sum, and Table." Therefore the dummy index (iterator) in a Sum is shielded ...
0
votes
0answers
34 views

Creating a non-zeta-ified sum

The expression Sum[i^(-s), {i, 1, ∞}] evaluates to Zeta[s] But, of course, this is not strictly correct as ...
0
votes
1answer
87 views

Accessing the elements in a table is slower than calculating each time?

This is a follow-up of How to accelerate combinations and sum calculations here? I have tried all tricks but when N reaches 35 the time is still intimidating in a 16-core server, which is far better ...
4
votes
0answers
62 views

Can a single Sum with multiple iterators be different from nested Sums?

Multiple sums are documented with two or more iterators, for example: Sum[1/(j^2 (i + 1)^2), {i, 1, Infinity}, {j, 1, i}] however the same answer can be ...
2
votes
1answer
56 views

How to accelerate combinations and sum calculations here?

I have a 4d matrix H defined as below (embedded with combinations and sums) ...
0
votes
2answers
50 views

How do I add the nth element for every line in a text file?

Suppose I have a text file, in this format: ...
1
vote
0answers
26 views

Partial Sum of Binary Sequence not Working

I have the following code: ...
0
votes
0answers
33 views

Problem of Summation

I want to calculate something like this: $$w_t = \sum_{i+j=t}c_{5-i,4-j}.$$ for $t=0,1\cdots8$. How I can write the code? I assume that $i$, $j$ and $t$ are all positive integers. Thank you very ...
2
votes
1answer
67 views

Working with tensor algebra

My question is really easy for experienced users. In my tensor equations I have an unknown tensor Q (symmetric and traceless): $Q=\begin{pmatrix} n1(x,y) & n2(x,y) \\ n2(x,y) & -n1(x,y) ...
1
vote
2answers
78 views

Double sum in mathematica

It may look like a stupid question, but I was really confused by how double sums works in mathematica: After simplification, the function actually has nothing to do with n2? Also, when I manually ...
2
votes
0answers
54 views

Why can Mathematica compute numerical sums more efficiently when they are written as matrix operations?

Let $f(n)$ and $K(n,m)$ be functions such that the double sum, which we wish to evaluate numerically, $$ \sum_{n=1}^a \sum_{m=1}^a f(n) f(m) K(n,m) $$ exists when $a$ is some large positive number. I ...
0
votes
1answer
39 views

Sums with index expressions instead of variables

Often I replace $k$ with, say, $m-k^2$ in a sum, obtaining something like this: $$\sum_{0\leq m-k^2\leq n} f(k)$$ Is there a way to input these without manually solving the inequalities?
3
votes
1answer
97 views

Speed up computation of sum from large matrix

I have a 3x3 matrix (252^3) of data (densities) and I want to compute a correlation xi from nearest neighbours as well as a sum involving a check whether the density is in a certain bin. The ...
1
vote
0answers
56 views

Assumptions aren't working in this sum

I have a sum that should be real: FullSimplify@Sum[1/(α n^4 + β), {n, 1, ∞}, Assumptions -> α > 0 && β > 0] But the result have involved even ...
3
votes
1answer
54 views

Equivalent (?) definitions of function gives different answers

I have a function as follows: ...
4
votes
1answer
83 views

double sum with condition

I was trying to compute the following sum: $$\phi(x,y)=\frac{1}{4\pi^2}\sum_{(n,m)\neq(0,0)}\frac{1}{n^2+m^2}\exp(i(nx+my))$$ where the range of indices is, say, $-10\leq n, \,m\leq 10$. But I don't ...
7
votes
1answer
285 views

Differing answers when comparing Wolfram Alpha and Mathematica v.10.2

Out of curiousity, please consider following expression: Sum[(-1)^(n + 1)/n, {n, 1, 100000}] When evaluated using Wolfram Alpha: Result: ...
0
votes
0answers
53 views

Sum of Infinite convergent series

I have a series as Sum[(x^1 - (x - 1)^1) (b^Log[(1 + (x))]), {x, 1, Infinity}] which is convergent when b=0.3 and ...
2
votes
2answers
56 views

Why does ExpandAll not work within a symbolic sum?

If I evaluate Sum[(x + Subscript[y, n])^2, n] + (y + z)^2 // ExpandAll then the expression within the Sum is not expanded, yet ...
8
votes
3answers
179 views

Expand series unevaluated

I'm newbie in Mathematica. I'd like to obtain nice and verbose output for any series calculation. For example, given a simple sequence n*(-1)^(n-1) and ...
3
votes
2answers
136 views

Summing matrix products

I need to compute a double sum over a weighted matrix product: $L[M]=\sum_{i,j}^{N}\Lambda[[i]]\;\omega[[i]].M.\omega^{\dagger}[[j]]$. $\Lambda$ is a list of with N complex values(weights) and ...
9
votes
1answer
131 views

SumConvergence[((-1)^n)/(Sqrt[n] + (-1)^n), n] returns True in Version 10.2?

Bug persisting through 10.4.1 I claim that the series $\sum_{n=2}^{\infty}\frac{(-1)^n}{\sqrt{n}+(-1)^n}$ diverges. To see this, rewrite the $n^{th}$ term as follows: \begin{equation*} ...
3
votes
2answers
87 views

DifferenceRoot question

I was doing the following sum: $$\sum_{i=2}^k \frac{(-1)^i}{i-1} \binom{2k-i-1}{k-1}x^i$$ First, Mathematica simplifies it to some DifferenceRoot function: ...
0
votes
0answers
41 views

Wolfram Mathematica: Sum of a function over a list [duplicate]

I have the function: $I_q(n)=\frac{1}{n}\sum_{d|n}\mu(d)q^{n/d}$, where $\mu(d)$ is the Möbius function. How do I get it in Wolfram Mathematica? It should be something like that: ...
1
vote
1answer
118 views

Summing infinite series that converge only for some parameter values

The input Sum[d^t,{t,0,Infinity}] produces output 1/(1-d) which is correct for $|d|<1$. But for $|d|\geq1$ the sum does ...
0
votes
1answer
52 views

combining function and plotting

How do i merge these 3 function onto 1 set of function to plot? Construct the series representation of a function $f(x)$ using up to $N0$ terms: ...
0
votes
0answers
30 views

Simplify expression?

Suppose I have expressions that involve terms like: p1.p1+p2.p2+p3.p3+p4.p4 p1, p2, p3, and p4 are all elements of an array: ...
5
votes
2answers
244 views

Speed-up the computation of this sum of small matrices

Given a matrix mat with a large number of rows (a few thousands) and a few columns (between 2 and 10), I'd like to compute the sum of the "small" matrices obtained ...
5
votes
2answers
122 views

Problem summing an infinite series

Calculating this sum on Mathematica 10.3 Sum[(-1)^(r - 1)/((a^2 + r^2)r), {r, 1, Infinity}] gives the answer ...
12
votes
3answers
357 views

Double Sum Involving Condition

I would like to compute the dimensions of some small free nilpotent Lie algebras. However, I am totally new to this and I could not figure out how to write the double sum which gives the dimension of ...
5
votes
3answers
227 views

sum involving binomial coefficient

To evaluate $ \sum_{\substack{\text{even }k\\ 2\leq k\leq m-2}} \binom{m-2}{k-1}(k-1) $, where $m$ is an even integer. Running the following in Mathematica 10.3.0 (October 9, 2015) ...
3
votes
1answer
88 views

Euler-Maclaurin summation

I want to compute asymptotic approximations to partial sum of harmonic series in Mathematica, using Euler-Maclaurin summation formula. ...
0
votes
1answer
147 views

Plotting a Double Sum

I am attempting to plot a graph of a function that is the absolute value of the double sum of an exponential function, and I keep getting errors. Here is the code: ...
0
votes
1answer
71 views

Slow sum computation

I'm doing a simple OLS regression and I strikes my attention that this very simple computation takes several seconds to perform in Mathematica ...
5
votes
3answers
139 views

incorrect exact sum of a series

Sum[(-1)^m*(1/Binomial[2*m + 2, 2] + 1/Binomial[2*m + 3, 2]), {m, 0, Infinity}] gives the right value -2 + π but the same ...
1
vote
2answers
67 views

Engel expansion

How do I do an Engel expansion in Mathematica? And find the unique non-decreasing sequence of positive integers $ \{ a_1 ,a_2,a_3,\dots \}$ , such $$\frac{e}{\pi}=\frac{1}{a_1}+\frac{1}{a_1 ...
5
votes
3answers
240 views

Make a vector of sums of matrix rows

I have a matrix in Mathematica: ...
0
votes
0answers
21 views

Adding a list up to a starting value [duplicate]

The code that i need has a starting value of 5. From this starting value on a list of numbers is added consecutively like this list {-0.1,-0,3,-0.1,0.1,-0.1} ...
4
votes
2answers
62 views

Controlling evaluation in a custom summation function

I'm doing some computations that involve summations with many indices over the same range - consider the simplified example Sum[a[i,j],{i,1,2},{j,1,2}]. To prettify ...
4
votes
1answer
119 views

How to increase the evaluation speed of this somewhat complicated table / matrix operation?

I'd like to calculate the following one-dimensional array: OneDimArray = Table[(Norm[Sum[ MatrixA[[i,j]]VectorB[[j]],{j,N}]])^2,{i,N}] However this takes a very ...
2
votes
1answer
69 views

Constrained Sum evaluation [duplicate]

I need to sum up F[x1,x2,...] over x1,x2,... where I have a constraint: x1+x2+...=n (n=positive integer number). For example, I have F[x1,x2] and x1+x2=3, so for the sum I get F[3,0] + F[2,1] + F[1,2] ...
33
votes
3answers
605 views

$\sum_{k=1}^{\infty }\left\lfloor\frac{5}{5^k}\right\rfloor$ giving wrong answer?

Bug introduced in 7.0 or earlier and persisting through 10.4 or later When I try to evaluate the following: $$\sum_{k=1}^{\infty }\Bigg\lfloor\frac{5}{5^k}\Bigg\rfloor$$ using ...
0
votes
2answers
82 views

plotting summation to check for convergence

Our assignment was to enter in a series and check to see if it converges by graphing. However, when I attempt to graph it, I get a bunch of errors that I'm not sure what they mean. Best way to show ...
5
votes
0answers
71 views

What is the underlying algorithm to simplify sums of reciprocals of polynomials?

Flipping through Wolfram's blog entry on Leibniz, W noted Huygens' interview test for the young Leibniz, namely to determine: $$ \sum_{n\ge2} \frac{1}{{n \choose 2}} $$ It's one thing to do this by ...
0
votes
0answers
53 views

Extract coefficients from sum

I would like to extract coefficients of indexed variables from a large sum. The following is an example motivated by linear regression to illustrate the problem. Consider the (unnormalised) log ...
1
vote
2answers
80 views

Sum of a set of elements in a array?

For a function I do this (Integrate over circular range): NIntegrate[x*y Boole[Sqrt[x^2 + y^2] < 3], {x, 1, 5}, {y, 1, 5}] ...
5
votes
2answers
164 views

Why this upvalue doesn't escape from Sum?

Probably a hard question, but I decide to cry out loud :). This is actually another problem I encountered when answering this question. Consider the following transform rule stored as an upvalue: ...
2
votes
5answers
77 views

tabulating summation series with restrictions

Trying to solve this question for a good 5 hours but remained stuck. My intuition tells me I am missing something trivial that I've overlooked or forgotten. I am restricted only to the function ...
1
vote
0answers
46 views
1
vote
5answers
190 views

Getting $e$ closer than $0,001$?

Everyone knows that $\sum \frac{1}{n!} =e$. How should I make a program in Wolfram Mathematica, that tells me, how many members I need from the sum, to get a number, with mistake smaller, than ...