For questions about symbolic computation, as opposed to numerical computations.

learn more… | top users | synonyms

5
votes
3answers
435 views

Why the difference?

When I do the double sum using the sigma notation I get $$1 + \sum_{n=0}^{\infty}\sum_{k = n}^{\infty} \frac{1}{(k+2)k!}$$ $1 + e - \cosh[1]$ When I do the sums as below, I get the expected ...
1
vote
3answers
924 views

Sqrt — how to get negative branch?

One of the many attractions of Mathematica is that integers are generally treated as exact symbolic entities, and not just as numbers. I was just playing with Sqrt, and was puzzled by what appeared to ...
1
vote
1answer
195 views

Why is a non-convergent series simplified to zero?

Sum[Sin[n*x]*Sin[n] /. x -> 1, {n, 1, Infinity}] == Sum[Sin[n*x]*Sin[n], {n, 1, Infinity}] /. x -> 1 Sum::div: Sum does not converge. >> Sum::div:...
4
votes
1answer
1k views

Non commutative multiply- expand expression

I began to use Mathematica a few days. My problem is: How to expand expression like $(a+b)*(a+b)$, where this multiplication is non commutative? Mathematica can do this?
1
vote
1answer
554 views

Generating a list of symbolic coefficients [duplicate]

I am trying to do something which I assume should be very simple: I want to generate a list with symbolic coefficients that looks like this: ...
1
vote
2answers
926 views

Dropping Higher Order Terms in symbolic evaluation

I have a series of symbolic manipulations (along with some coefficients). As a sample, here is one of the intermediate forms ...
0
votes
0answers
41 views

Symbolic Calculations with formal matrices [duplicate]

I want to make calculations with matrices in the form (for lovers of statistical mechanics, is just the first step replica symmetry breaking) ...
4
votes
1answer
592 views

Is it possible to create a compiled function with some symbolic arguments?

I am trying to create a compiled function that takes in several arguments. However, some of the arguments contain symbolic entries and thus I get the following error message when executing the cell ...
2
votes
1answer
160 views

Differentiating with respect to vectors of unspecified length

Given a function defined by f[x_,y_] := Sum[ x[[i]] (g[y[[i]]] - g[Total[y]] - Log[x[[i]]]), {i, 1, Length[x]}] I can easily calculate the derivative for a fixed ...
1
vote
1answer
154 views

Only perform a symbolic differentiation once

I want to define a function that involves a differentiation step that Mathematica can do easily, which might be of the form ...
0
votes
1answer
234 views

Alleviating a bottleneck involving a function call with symbolic derivatives

I have a bottleneck problem involving a function call with symbolic derivatives. I don't have much experience with Mathematica (or any dynamic language) and I'm certain that the way I've set this up ...
10
votes
2answers
919 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 ...
11
votes
4answers
2k views

Extracting variables from an expression

I'm looking for a way to extract a list of variables from an expression, for example with an input like: ...
13
votes
1answer
408 views

What is the confidence limit on this convergence?

Bug introduced in 7.0 and fixed in 10.0.0 When I run this, Product[n^MoebiusMu[n],{n,1,Infinity}] I get $\frac{1}{4 \pi^{2}}$ Over on Math Overflow ...
3
votes
1answer
306 views

Sum of positive terms gives negative answer

Bug introduced in 7.0 and fixed in 9.0 Mathematica evaluates Sum[((n - y - 1)*(n - y)^2*n^y)/y!, {y, 0, n - 2}] as -2 e^n n. ...
7
votes
3answers
581 views

Efficient code for solve this equation

We have $a*b*c=-1$, $\frac{a^2}{c}+\frac{b}{c^2}=1$, $a^2 b+a c^2+b^2 c=t$ What's the value of $a^5 c+a b^5+b c^5$? I tried ...
4
votes
1answer
169 views

Identical code, different answers?

I'm having some trouble with identical code giving different answers. On a fresh kernel (MM 9.0.0.0, Windows 64-bit), running the same code, copy-paste, gives two different answers: ...
4
votes
0answers
180 views

Version 8.0 integrates but Version 9.0.1 doesn't

I am trying to run the following integral in version 9.0, but it fails: ...
2
votes
1answer
1k views

Variable name string to expression

I'm not sure how to name this problem appropriately. But anyway. I have, say, three variables (arrays in my case), which I call: A1, A2, A3 I would like to ...
6
votes
1answer
2k views

Non-commutative algebra

I'm constantly dealing with non-commutative algebras. ** is inbuilt, non-commutative and associative. That's good :-) But it is not distributive. Rats. ...
2
votes
3answers
2k views

How to calculate partial derivative of an unknown function

How can I force Mathematica to calculate symbolically the partial derivative of a function u[x,y] with respect to a variable ...
4
votes
3answers
571 views

Symbolic derivative of $n$-term product

I want to determine the relationship that must exist between the $x_i$ and $y_i$ such that $$ \frac{\partial}{\partial\theta} \prod_{i=1}^n \frac{f(x_i,\theta)}{f(y_i,\theta)} = 0, $$ where $$ f(x,\...
3
votes
1answer
290 views

What properties make this equation difficult to symbolically solve?

I tried to solve the following equation with Mathematica's solve: Solve[K*(2*Tan[L/2*Sqrt[P/(EI)]]-L*Sqrt[P/(EI)])+4*P*Sqrt[P/(EI)] == 0, P] It gave the ...
2
votes
1answer
283 views

Can I reduce a matrix inequality such as $\mathbf x^\prime\mathbf A\mathbf x > \mathbf x^\prime\mathbf x$?

I'm new to Mathematica. When I do linear algebra, I wonder if I can have an inequality such as $\mathbf x^\prime\mathbf A\mathbf x > \mathbf x^\prime\mathbf x$, where $\mathbf x$ is a column vector ...
4
votes
0answers
96 views

Using Resolve and ForAll to prove takes a really long time

I've been trying to prove a lemma for my paper using Mathematica... basically that $$\forall \{n, d_i, d_j\} \in \mathbb{Z},\ n \ge d_i > d_j \ge 2$$ it's true that $$V[1, n, d_i-1, d_j-1] >...
5
votes
1answer
882 views

symbolic integration of a definite integral, how to speed up

I am wondering what are the general tips & hints for speeding up symbolic integration. As far as I understand, it solely depends on the form a function is represented before it is processed by <...
6
votes
1answer
482 views

Simplification of double symbolic sums containing a DiscreteDelta without explicit summation range

I am trying to get Mathematica to automatically do simplifications like the following: $$\sum\limits_{q}^{q\in qV}\sum\limits_{q'}^{q'\in q'V}{f(q)g(q')\delta(q-q')}=\sum_{q}^{q\in qV}{f(q)g(q)}.$$ ...
5
votes
1answer
199 views

A triple sum related question

I'm trying to compute the triple sum Sum[ 1/(i! j! k! ), {i, 1, Infinity}, {j, i + 1, Infinity}, {k, j + 1, Infinity}] but Mathematica doesn't return any value. ...
3
votes
1answer
662 views

Nested Sums to multiple sum

I would like to automatically "move nested sums to the left". I mean, just take out of an expression all the summations and go from a nested Sum to a multiple sum. Something like starting with: ...
2
votes
1answer
681 views

How can I solve an exponential inequality in Mathematica?

I would like to solve the equation: A Exp[m t + s Sqrt[t] x - 1/2 s^2 t] >= K with respect to x (which is in the exponent) ...
11
votes
5answers
2k views

How to get exact roots of this polynomial?

The equation $$ 64x^7 -112x^5 -8x^4 +56x^3 +8x^2 -7x - 1 = 0 $$ has seven solutions $x = 1$, $x = -\dfrac{1}{2}$ and $x = \cos \dfrac{2n\pi}{11}$, where $n$ runs from $1$ to $5$. With ...
12
votes
1answer
268 views

Symbolic computations with already assigned variables

For example I have assigned some variables that could defined in terms of each other. Then I for some reason need to derive some formula using this variables. How to show the result of this ...
5
votes
4answers
2k views

How to find the sum all even numbers of this sequence?

I have a sequence $(u_{n})$ $$u_1= 1, \quad u_2 = 2, \quad u_3 = 3, \quad u_{n}= -u_{n-3} + 3u _{n-2} +2 u_{n-1}, \quad \forall n \geqslant 4.$$ I want to list the first $20$ terms of this sequence ...
0
votes
1answer
807 views

How to efficiently fill in matrix values?

I would like to efficiently fill in or assign values to matrix m: m={ {a[1,1],a[1,2],a[1,3],a[1,4]}, {a[2,1],a[2,2],a[2,3],a[2,4]}, {a[3,1],a[3,2],a[3,3],a[3,4]}, {a[4,1],a[4,2],a[4,3],a[4,...
3
votes
1answer
534 views

Mathematica 9 and later behavior with derivative of a sum

Bug introduced in 9.0 or earlier and persisting through 10.3.0 or later ...
1
vote
1answer
231 views

Symbolic calculation

Is there any software where this kind of symbolic calculation is possible? s = 0 for i = 1 to m j = 2 i + a - t if (j > 0) s = s + a I want s as a ...
2
votes
1answer
229 views

Functions that operate on symbolic matrices?

I'd like to write functions that operate on symbolic matrices, and do nothing when fed anything else. ...
9
votes
1answer
382 views

Verifying and deriving basic (block) matrix identities

How can I use the new symbolic matrix/tensor capabilities to verify matrix identities, such as (1) or (2) Even better, how can I ask Mathematica to derive expressions for X, Y, Z, and U like ...
17
votes
3answers
2k views

Is learning to use Mathematica useful for pure theoretical research in Mathematics and Computer Science?

I am looking for opinions from Mathematica users about Mathematica itself. After reading the faq, I thought that in some sense I "wish to solve a problem using Mathematica" ... although I understand ...
4
votes
1answer
1k views

Computation of parametric integral

I am trying to compute the integral Integrate[(g^(u^(g - 1)))/(1 + u^g), {u, 0, t}] but as an answer I get my input expression. There must be something wrong ...
1
vote
1answer
214 views
0
votes
1answer
148 views

Industrial Level Applications. Recipe for mixed notation of equations set

I am working with large (linear) equations set within Mathematica in numerical notation. For example, set from 4056 eq. is solved for a second, no more. There is no doubt, result is great. But even ...
16
votes
1answer
521 views

Speeding up trigonometric integral

Context On a possible non trivial toric topology for the Universe (nothing less!). Problem I would like to carry out the following integral for $\ell=2,4\cdots 20$. $$\int _0^{\pi }\int _0^{2 \pi }...
8
votes
6answers
2k views

Prevent Part[] from trying to extract parts of symbolic expressions

If you have a list, e.g. {1, 2, 3} then you can extract the $k$th part using Part (...
2
votes
1answer
253 views

Difference between Map[f[#] + g[#] &, {a, b, c}] and Map[f[x]+g[x]&, {a, b, c}]

Map[f[#] + g[#] &, {a, b, c}] vs Map[f[x]+g[x]&, {a, b, c}] My question is: why is the output different? ...
10
votes
1answer
484 views

Symbolic Integration along contour: branch cut problem?

Context Following this question on path integrals in the complex plane, having defined again a numerical and symbolic integrator along a path as ...
20
votes
2answers
5k views

How do I work with Root objects?

I want to solve the trigonometric equation : $$(3-\cos 4x )\cdot (\sin x - \cos x ) = 2.$$ I tried Solve[(3 - Cos[4*x])*(Sin[x] - Cos[x]) == 2, x] It returns the ...
5
votes
1answer
1k views

Symbolic integration in the complex plane

Context While answering this question, I defined (symbolic and numerical) path integrations as follows ...
3
votes
0answers
247 views

Numerical-Symbolical Integration (Calculus)

I created a simple numeric-symbolic integration. Here you can use symbolical and numerical techniques at the same time. You can also interpolate numerical integrals. The problem with my function is ...
9
votes
5answers
2k views

Superscript prime symbol

x\[Prime] looks like $x_{'}$, ugly right? Is there a way to make a symbol with prime to look like $x'$? That's what I'm trying right now: ...