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

learn more… | top users | synonyms

2
votes
3answers
86 views

Symbolic integration conditional bug?

I just ran across this: Assuming[x > 0, Integrate[1/t^2, {t, 1, x}]] gives the output ...
2
votes
1answer
65 views

How to define the symbolic expression $V = \frac{1}{N} \sum_{i=1}^N ( y_i(g) - z_i)^2$ in Mathematica?

I want to define this expression symbolically: $V = \frac{1}{N} \sum_{i=1}^N ( y_i(g) - z_i)^2$ in which $\{y_i(g)\}_{i=1}^N$ is to be understood as a symbolic sequence that depends on a vector $g$. ...
1
vote
0answers
69 views

Summing the probability distribution to 1 to convince myself

I have worked out a probability distribution and want to check its sum which is necessarily 1. First we write $$ r \triangleq \frac{(2 \lambda + \mu)^2}{2(\mu + \lambda)^2}, \quad s \triangleq ...
0
votes
1answer
38 views

Taking function output as shown, and processing it regardless of its history

I am making some symbolical calculations for my quantum chemistry class. I am supposed to build a wave function composed of multiple other functions, each of which is also composite of some lower ...
1
vote
1answer
113 views

Why do expressions with D and Derivative giving so different results for me?

I tried to express with Derivative the same I expressed earlier by D, but failed: ...
4
votes
2answers
137 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 ...
2
votes
1answer
99 views

How to get a more compact form of this probability calculation?

Inspired by the probability calculation here, I am trying to solve a little general one: $$\mathbb{P} (\sum_{i=1}^{m-1} A_i + \sum_{i=1}^{m} S_i < L < \sum_{i=1}^{m} A_i + \sum_{i=1}^{m+1} ...
5
votes
3answers
255 views

Comparison of Dispatch objects ureliable in Mathematica 10

In 2008 I wrote a group theory package. I've recently started using it again, and I found that one (at least) of my functions is broken in Mathematica 10. The problem is complicated to describe, but ...
2
votes
0answers
59 views

How does Mathematica resolve symbolic systems of inequalities?

Before getting into my (probably annoying) question, let me provide some context: I recently started using Mathematica to automatically perform stability analyses that I previously did by hand. A ...
8
votes
2answers
184 views

How to do algebra on summations of variable expressions?

Is there any easy, general way to handle algebra on variable Sums? For examples, Gaussian Mixtures or Fourier Series: ...
0
votes
2answers
55 views

DSolve bvlim error due to Boundary condition at infinity

I am solving an Ode using Mathematica that has boundary condition at infinity which is why Mathematica is showing error message. Please help to tackle this problem. ...
5
votes
1answer
267 views

How to solve this probability symbolically or numerically?

I am trying to calculate the following probability $$\mathbb{P} \big(\sum_{i=1}^{m} (A_i + S_i) \le L < \sum_{i=1}^{m+1} (A_i + S_i) \big)$$ where, $$A_i \sim \exp(\lambda), \quad S_i \sim ...
-1
votes
0answers
27 views

Solving Symbolically in Mathmatica 10.1- [duplicate]

I'm trying to solve this equation symbolically for x :sqrt[29991348*x^2 - 13958305552*x + 992000000 + ((a - b)/x)/(0.02 + c - a) = 0. However, when I try ...
1
vote
2answers
44 views

Symbolic integration of a indexed product term?

I am trying to perform the symbolic integration of the following expression $$\int \prod_{i=1}^n p(x_i | z) p(z) dz$$ where $$p(x_i|z) = \frac{exp(\frac{-(x-z)^2}{2\sigma_i^2})}{\sqrt{2 \pi ...
0
votes
2answers
74 views

Symbolic matrice or vectors in Mathematica?

Suppose I would like to define an expression, containing and array or matrix of symbolic elements. Like this ...
1
vote
0answers
37 views

Symbolic Nullspace computation in parallel

I am trying to determine the nullspace of a large symbolic matrix. The single core evolution seems to take too long. I tried Parallelize[] on the commands ...
1
vote
0answers
44 views

Symbolic and numeric calculations (and plots) simultanuosly

I use Mathematica to do a bunch of symbolic calculations (integrals, ...). This is good because I found that sometimes, if I plug in numeric values, Mathematica takes much longer. However, sometimes ...
9
votes
0answers
85 views

Backslide of Limit

A friend of mine showed me this example: Limit[Sum[Sin[Pi*k/n]/(n + 1/k), {k, 1, n}], n -> Infinity] This sample calculates well in v8.0.4: but not in ...
3
votes
1answer
77 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 ...
0
votes
1answer
66 views

Wrong Definite Integral result [closed]

When I calculate the definite integrate of (a+bt)*cos(nt). Integrate[(a + b t) Cos[n t], {t, t1, t2}] MMA gives a correct answer: ...
0
votes
2answers
44 views

How to use Expected values in symbolic solving

How can I create a function that brackets off terms of an expectation value? For example, I have : f = a[x] + c1*b[x] + c2*c[x]; g = Expand[f*f]; ...
0
votes
1answer
44 views

Integral evalutation

I am trying to integrate the following expression over $L$. E^(-2 L n) (1 - L/(2 s))^(-1 + 4 n µ) (L/s)^(-1 + 4 n v) I did... ...
1
vote
0answers
38 views

Why does integration of a radical times HeavisideTheta give a conditional expression?

If I want to compute an integral over x^p HeavisideTheta[y-x] where p is some fraction, I do not get the desired result, but ...
14
votes
6answers
527 views

How can I evaluate only a single step of a recursive function?

Let's say have a simple recursive function for the Fibonacci sequence f[0] := 1 f[1] := 1 f[n_] := f[n - 1] + f[n - 2] but I want to see how it will expand in a ...
1
vote
1answer
65 views

Strange integration

I tried to evaluate this line Integrate[((-I/2) (E^((-I) x) - E^(I x)) + ((E^((-I) x) + E^(I x)) x)/ 2)/((-1 + E^x) x), {x, 0, Infinity}] Then I get $14$ ...
4
votes
2answers
102 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 ...
3
votes
3answers
681 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: ...
5
votes
0answers
145 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): ...
0
votes
1answer
73 views

About Symbolize and Variables

I do believe this sort of similar questions had been asked many times, so I have read a LOT. For example this and this as well as many other related posts. But a glance at what my main question is ...
3
votes
2answers
155 views

How to define even permutations correctly?

I define even permutations as following, but there may be some error. I use it in two different way and get different output. ...
0
votes
1answer
56 views

Algebraic manipulation of closed loop state space system [closed]

I have the following equations $$ \begin{align} x(k+1) &= A x(k) + B \tilde{u}(k) \\ y(k) &= C x(k) + D \tilde{u}(k) \\ w(k+1) &= E w(k) + F \tilde{e}(k) \\ u(k) &= G w(k) + H ...
1
vote
0answers
25 views

Derive logistic choice probabilities symbolically

More generally, I am interested in learning what the current limitations of Mathematica are when using it for doing pure mathematics. A recent blog post by Stephen Wolfram (see: ...
5
votes
0answers
145 views

Is it possible to implement symbolic sub matrices?

The accepted answer in the following link: http://stackoverflow.com/questions/5708208/symbolic-matrices-in-mathematica-with-unknown-dimensions provides a functionality to create symbolic matrices ...
0
votes
1answer
52 views

A problem with getting a coefficient

I try Coefficient[E^(I a (t - b)), E^(I a t)] and expect the output E^(-I a b) but in fact I get ...
2
votes
1answer
39 views

A question about expanding complex functions of real arguments

I have the following problem. There is such an expression as: P[x_,y_] := z[y] E^(I beta x) + Conjugate[z[y] E^(I beta x)]; The variables ...
1
vote
1answer
62 views

How to simplify symbolic matrix multiplication using the associativity of scalar multiplication?

My input is (2A).(3B). Is there a way to get 6A.B? All I could get back from Mathematica was ...
0
votes
0answers
61 views

No solution with DSOLVE

I have been trying to solve the following system of ODEs: \begin{eqnarray} ...
0
votes
1answer
192 views

Is this function convex?

How can I determine convexity of the function f = Log[ x, 1 + (x^a - 1) (x^b - 1)/(x - 1)] with the parameters $a,\,b$ belonging to the interval $(0,1)$ in ...
5
votes
2answers
289 views

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

Update I have tried like these. I think there is a bug. Plot[1/Sqrt[-1 + 2^2 Sech[x]^2], {x, 0, ArcCosh[2]}, Ticks -> {{ArcCosh[2]}, Automatic}] This ...
1
vote
0answers
31 views

Symbolic computation under certain assumptions [duplicate]

Assume that I have a matrix $M$ and I'm doing some operations on it. Is there a way to specify that two variables satisfy an equation say $x=y^2$ and then let Mathematica compute under that ...
8
votes
1answer
6k views

Extract real part of a complex expression better than Re does

I have a complex expression with real positive variables only. Mathematica Input Style: ...
6
votes
1answer
64 views

How to assign up-values for `Derivative`?

I have defined several custom analytic functions. Here is the simplest example: ln[x_, a_?NumericQ] := Piecewise[{{Log[x], Re[a] > 0}, {-Log[1/x], True}}] ...
18
votes
2answers
480 views

symbolic integration error

I'm running Mathematica 10.0.0 and encountered a disturbing error in the symbolic integration of a rather simple function ...
3
votes
1answer
80 views

Pulling powers to the outside of expressions (inverting PowerExpand)

Is there something like an inverse of PowerExpand that pulls powers to the outside of expressions: $$ x^2 y^2 \leadsto (x y)^2 \\ x^2/y^2 \leadsto (x/y)^2 $$ I ...
1
vote
0answers
96 views
7
votes
1answer
105 views

How to find a recurrence relation for a sequence?

I have a sequence given by an explicit formula for n-th term: ...
6
votes
1answer
159 views

Symbolic tensor simplifications and the identity matrix

How can I get Mathematica to simplify the following expressions, with $Assumptions including Element[n,Integers], n > 0, and ...
1
vote
0answers
58 views

Typing (and executing) expressions with multiple superscripts

I am trying to type (and evaluate) expressions of the following form: $$ G^{a,b} $$ into a mathematica. I've tried the obvious G^(a, b) or ...
18
votes
0answers
478 views

Proving the hairy ball theorem using xAct

I would like to formally prove the hairy ball theorem in Mathematica, initially just for $S^2$, and then see about generalizing. An approach I thought about to use the xAct package to define $S^2$ ...
3
votes
2answers
95 views

Expectation of a composite Markov-Gamma distribution

In a model I have a discrete two-state first order Markov process, defined by a (2x2) transition matrix with two free parameters. If the first state occurs then the process outputs zero for that ...