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

learn more… | top users | synonyms

1
vote
1answer
52 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
66 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
35 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
111 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: ...
3
votes
2answers
133 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
248 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
57 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 ...
7
votes
2answers
169 views
+100

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
73 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
82 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
76 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
526 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
677 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
144 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
71 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
153 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 [on hold]

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
143 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
191 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
288 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
477 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
56 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
477 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 ...
0
votes
0answers
82 views

Calculate the power spectral density of a Markov chain

I would like to calculate the symbolic power spectral density of a two state Markov process with a symbolic transition matrix characterised by two parameters. I have tried the below, but it doesn't ...