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

learn more… | top users | synonyms

5
votes
1answer
123 views

Checking if a symbolic matrix is positive semi-definite

As far as I can tell, there is no way to tell the builtin PositiveSemidefiniteQ about assumptions on symbols. For example, the matrix ...
5
votes
1answer
135 views

MarginalDistribution with Symbolic range in ProbabilityDistribution

Given the following joint density: for 10 < x < 20, x/2 < y < x. To find the marginal density for X we do: ...
5
votes
1answer
436 views
5
votes
3answers
309 views

Running out of variables

I am running symbolic calculations with Mathematica. Simplification of algebraic expressions is an important part of them as very often final results can be simplified considerably. Typically my ...
5
votes
1answer
84 views

Why is NHoldFirst not propagated to symbolic derivatives?

I encountered a nasty problem that N cannot evaluate expressions containing a symbolic Derivative of a multi-parameter function ...
5
votes
2answers
169 views

How to check an algebraic number for membership in a list

I need to check an algebraic number for membership in a list of algebraic numbers. The numbers can be expressed in different forms (combinations of radicals, Root ...
5
votes
1answer
69 views

Rearrange an algebraic expression so that each variable appears only once

Mathematica gave me a solution to an equation in this form: ans1 = (b c)/(a - a c + b c) $\text{ans1} = \frac{b c}{a - a c + b c}$ From solving the same equation ...
5
votes
3answers
369 views

Different results of a definite integral $\int_0^{\cosh ^{-1}(a)} \frac{1}{\sqrt{a^2 \text{sech}^2(x)-1}}\, dx$

fixed in 10.0.2 Update I have tried like these. I think there is a bug. ...
5
votes
1answer
53 views

Simplify complex answer given by DSolve[]

I tried the technique in this question, but that did not work for the solution given by the MWE: DSolve[{f'[t] + f[t]*g[t] == h[t], f[0] == f0}, {f[t]}, t] How ...
5
votes
1answer
74 views

How to bail out from a non-applicable transformation function?

I provide a custom function in TransformationFunctions options of FullSimplify that can do a transformation only for some values ...
5
votes
1answer
247 views

VariationalD giving the wrong solution?

EDIT: As pointed out in the comments, VariationalD gives a variational derivative (which I don't want), not a derivative with respect to a function (i.e. ...
5
votes
1answer
372 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
136 views

Slow-down encountered when using NMinimize with a compiled function

I'm minimizing a large function. Evaluating the function is very slow symbolically but using Compile I can get ~100x speed-up on this step. However, despite this, the NMinimize process is only ...
5
votes
1answer
375 views
5
votes
2answers
559 views

Problems with Symbolic summation over unknown values

I'm having some real trouble with Mathematica wrongly evaluating various symbolic sums at the moment. I have this function: $$h_{ij}(x) = ...
5
votes
2answers
443 views

Reducing exponential inequalities fails

I am quite stumped by this problem : Reduce[ N^(x-y) <1 && N > 0 , x, Reals] gives the expected result ...
5
votes
0answers
60 views

Operator which can be interpreted as binary and unary

I'm a bit lost with the way how e.g. the + operator is implemented in Mathematica as binary (infix) and unary (prefix) operator depending on the context, since I would like to define a similiar ...
5
votes
0answers
99 views

Integral formula for the inverse Laplace transform doesn't work?

The direct implementation of the definition of the inverse Laplace transform using Integrate fails in the following case: ...
5
votes
0answers
220 views

Integrate wrong for absolute value of trig function

I was trying to get $\int_0^1 \lvert \cos(2 \pi k x) \rvert \,\mathrm{d}x$ for $k \in \mathbb{Z}$, and was surprised by the result (using Mathematica 10.0.1.0): ...
5
votes
0answers
1k views

Solve huge symbolic system of linear equations

I have a system of 76 symbolic linear equations (i.e. some coefficients are symbolic) with a sparse coefficient matrix. However, neither Solve[] nor ...
4
votes
7answers
443 views

Given a symbolic expression how to find if starts with a minus or not?

I am using a Mathematica function which returns some error term in symbolic form. I needed a way to determine if this term starts with a minus sign or not. There will be only one term. This is to ...
4
votes
3answers
483 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 $$ ...
4
votes
3answers
406 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 ...
4
votes
2answers
127 views

Well-defined symbolic integral leading to ConditionalExpression

I would like to determine a closed-form expression for the following symbolic integral $$ \int_{-1/2}^{1/2} \!\!\!\! \mathrm{d} x \int_{-1/2}^{1/2} \!\!\!\! \mathrm{d} y \, \frac{1 + b x + c y}{1 + e ...
4
votes
1answer
252 views

Symbolic bit vectors

I'd like to see how addition and xoring bitvectors mix together. To do this, I implemented (a primitive) vec_add and vec_xor: ...
4
votes
1answer
683 views

Does the Im function work with symbolic arguments?

Does the Im function work with symbolic arguments? ...
4
votes
2answers
216 views

Why doesn't FullSimplify simplify expressions with DiracDelta?

I want to simplify a complicated expression with some Dirac delta distributions, but FullSimplify does not do what I want. Specifically, I want ...
4
votes
1answer
727 views

Definite and Indefinite integral give different results for piecewise function

I have the following function: $$ f(q,y)= \begin{cases} \tfrac{11720+p}{37791360} & -11720<p<-7720 \\ 0 & \text{True} \end{cases} $$ where $p = 443\ y-777600\ \sin^{-1}\left(\frac{q ...
4
votes
3answers
217 views

Why the inequality does not take into account the domain?

I have this inequality: Reduce[(4000-1000k)/(k-4) < 0] and the answer is k ∈ Reals I would expect ...
4
votes
1answer
173 views

Volume within parameter space

Imagine a parameter space with variable 0<p<1, 0<e1<1/2 and 0<e2<1/2. ...
4
votes
1answer
746 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?
4
votes
1answer
64 views

Factoring out parts of an expression

I'm looking to factor out the parts of a multiplication type of expression. For example for the expression: $$20 e^{-t^2}\frac{t^2}{\sqrt{t^2+a}} erf(b+t)$$ I'd like to get a list in the form of: ...
4
votes
1answer
134 views

Factoring a differential equation into an operator product

Following the fundamental theorem of algebra, I can (as stated here in Sec. 1.1.4) factor an $n$th-order linear ordinary differential equation $$ a_n \frac{d^n x}{dt^n} + a_{n-1} \frac{d^{n-1} ...
4
votes
2answers
322 views

How do I evaluate a symbolic integral involving Hermite polynomials?

I want to test a difficult integral : Integral on all reals of some complicated function involving the Hermitian polynomials, exponentials, squares, factorials, and being general considering any ...
4
votes
1answer
148 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
1answer
901 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 ...
4
votes
1answer
448 views

How to extract a possible closed form from WolframAlpha[] output

To find a possible closed form of a number, I can use the function WolframAlpha["6.38905609893065", IncludePods -> "PossibleClosedForm"] It returns a result ...
4
votes
2answers
248 views

Problem with adding vector to symbolic function (for NDSolve)

I'm trying to set up a system of differential equations for passing to NDSolve. Note that my initial conditions are vector valued so Mathematica should know that ...
4
votes
1answer
388 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 ...
4
votes
2answers
66 views

Using LeviCivitaTensor

I have 2 questions about using LeviCivitaTensor, based on the following session: The signs for the cross product seems to come out "opposite" of what I expected. Why is that? Did I miss anything? ...
4
votes
2answers
270 views

Coordinate free differential forms package

I am looking for a package in Mathematica which can handle differential forms in a coordinate free manner. I am aware of several packages which do differential forms, but it seems that for all of them ...
4
votes
1answer
315 views

How to check if a given expression is an “explicit algebraic number”?

The documentation for PossibleZeroQ says: With the setting Method...
4
votes
2answers
132 views

Remove redundant parameters from equation

I have random expressions like (b1 x + b2 x)/(b3 + b4/b5 x) + Sin[b6 x] where $x$ is a variable and {b1, ..., b6} are parameters to be fitted. This expression is ...
4
votes
1answer
101 views

Avoid writing explicit form of operator

I would like to evaluate the following (simplified) expression with Mathematica : $\frac{\delta}{\delta J} \exp[(J(x)+K(x)) \Delta (J(x)+K(x))]$ where $\Delta$ is a differential operator independent ...
4
votes
0answers
95 views

Symbolic integration of elliptic functions

Is there some clever way to integrate products of elliptic functions $\wp(z;g_2,g_3)$ or $\zeta(z;g_2,g_3)$ in Mathematica? Mathematica seems to be able to integrate functions like ...
4
votes
0answers
68 views

Does it help to do any preprocessing before `FindSequenceFunction`?

I am always amazed how smart FindSequenceFunction can be in discovering general term formula for a sequence of rational numbers. The documentation says it can work ...
4
votes
0answers
600 views

Analytically solve the eigenvalue problem with infinite dimensions by Mathematica?

If I am given a symbolic expression of all the matrix elements in an infinite-dimensional space, e.g., the Hamiltonian of a quantum mechanical system, is it possible to get the symbolic expression for ...
4
votes
0answers
95 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] ...
3
votes
2answers
481 views

Write a function that returns the logarithmic derivative

How can we write a function that if we input an expression f, it returns the log derivative $\frac{1}{f} \frac{df}{dx}$. We have to use a conditional or pattern test so that the function accepts any ...