Questions on the analytic and numeric equation solving functions of Mathematica (Solve, Reduce, NSolve, FindRoot, DSolve, RSolve, etc.).

learn more… | top users | synonyms (5)

5
votes
1answer
2k views

Error/warning when using NSolve for simple equation

I am using NSolve to solve an equation, as shown here: ...
5
votes
1answer
661 views

Find asymptotics of a function

Can Mathematica find the asymptotics of a function in the following sense? I have Log[1/n^2]/Log[4^(-Sqrt[Log[n]]) + (2^(-Sqrt[Log[n]]) - 1)^2] and I would like ...
5
votes
2answers
2k views

Solving a trig equation

I am trying to find an explicit solution to a trigonometric equation. It is running for a long time and I didn't get any result, but Maple solves it quickly. Why ? ...
5
votes
1answer
122 views

SystemException[“MemoryAllocationFailure”] Uncaught SystemException returned to top level from NSolve

Bug introduced in 10.0.0 When trying to solve 8 equations in 8 unknowns, I get SystemException["MemoryAllocationFailure"]. I am not sure if I am doing something ...
5
votes
2answers
200 views

Series solution of the Lane-Emden equation

I have the following problem: I'd like to show how the Lane-Emden equation would look like if someone would solve it as a series. Here's the code: ...
5
votes
1answer
225 views

Solution to $y = x^x (1-x)^{(1-x)}$ for $x$ [closed]

Does anyone have any insight to a symbolic Mathematica solution to this equation for $x$ in terms of y? Symbolic of any sort, just not numeric. $$y = x^x (1-x)^{(1-x)}$$
5
votes
2answers
416 views

Solve for x given an interpolation function, y

Given a list of data, lst={{x1,y2},{x2,y2},...,{xn,yn}} and the interpolation function fun, y=Interpolation[lst] which ...
5
votes
1answer
413 views

How to solve this Integral equation

D[x[t] - x[t - 1]/(2 E), {t, 3}] + Integrate[E^(-δ)*x[t - δ]/5^t, {δ, 2, 2.5}] == 0 I found solve this problem is hard with Mathematica. I also find a article ...
5
votes
3answers
890 views

Equation is solved quickly using Solve but takes too long using NSolve

I have the following simple command to solve for h in an equation: Solve[(1 - 2.25577*^-5*h)^5.25588 == 0.9644952131579817, h] ...
5
votes
2answers
291 views

How to get $x$ in terms of $y$ and $w$ from $w = a p x^{p-1} \left(a x^p+b y^p\right)^{\frac{1}{p}-1}$

I'm having trouble to solve the following equation for $x$ both manually and in Mathematica: $\quad \quad w = a p x^{p-1} \left(a x^p+b y^p\right)^{\frac{1}{p}-1}$ In Mathematica notation: ...
5
votes
2answers
382 views

Why doesn't Roots work on a certain quartic polynomial equation?

In everything that follows I am using Roots to get the solutions to the equations. I start off with this equation: $$ \frac{A}{3}x^4 - x^2 + 2ax - b^2 = 0$$ Now ...
5
votes
2answers
222 views

Reduce with Integer domain, but receive a rational solution

Bug introduced in 6.0 and fixed in 10.0 I want to find the integer values $k$ so that D[f[x], x] has integer solutions. I tried ...
5
votes
2answers
544 views

Intersection of two splines

I have a collection of sets of points, each set defining the segment of a (nice) curve on the unit sphere $S^2 \subset \mathbb{R}^3$. The points are computed numerically. I now want to compute the ...
5
votes
4answers
106 views

Using roots to get original equation

I have a list of four rules, that I would like to use to find the original polynomial that produced them. The rules are of the form ...
5
votes
3answers
122 views

On Solutions of Diophantine Equation

By Mathematica, I try to obtain solutions of Diophantine equations such as: $$F_{n_1}F_{n_2}\cdots F_{n_k}+1=F_t$$ where the sequence $\{F_n\}$ is the Fibonacci sequence,$n_1 < n_2 < \cdots &...
5
votes
1answer
186 views

$x^6-2\varphi^5x^5+2\varphi x+\varphi^6=0$ solving in radicals

How to solve the equation $x^6-2\varphi^5x^5+2\varphi x+\varphi^6=0$ in radicals using Mathematica? where $\varphi$ is the golden ratio. Solution is Ramanujan's Class Invariant $G_{125}$, Class ...
5
votes
2answers
927 views
5
votes
1answer
289 views

A problem about Steffensen method

Two days ago,I asked a question about bisection algorithm. Thanks for Michael E2 giving me a good solution My algorithm 1, I use the bisection algorithm to generate a initial value that contain a ...
5
votes
1answer
538 views

Solving a nonlinear, second order d.e. in Mathematica

I am trying to solve the following differential equation $$ s''(x) = \frac{8}{9}s(x) -2s^2(x) + s^3(x), $$ with boundary conditions $s(-30) = 4/3 $ and $s(30)=0$, in Mathematica by using the following ...
5
votes
1answer
3k views

finding poles of a function

Is there a command to find the poles of a function $f=f(z)$? example: let $$f(z) = \frac{1}{z^2-1}$$ then we know that the poles are at $z=\pm 1$ but is there a special command in mathematica to do ...
5
votes
2answers
2k views

Find all roots in range

Is there any way to check for all roots in a range? Jens' findAllRoots function is pretty good, but runs at approx. 10% of roots missed when I ran a quick check on ...
5
votes
1answer
174 views

NDSolve solution violates initial conditions

I have the following code: ...
5
votes
2answers
267 views

How to “ReplaceAll” using a set of conditional expressions instead of a list of rules [closed]

Suppose that I have two linear functions f[x_] := f0 + f1 x g[x_] := g0 + g1 x and a (possibly rather complicated) set of conditional expressions, obtained ...
5
votes
2answers
315 views

Computing the sign of the real part of eigenvalues in a 3D linearized system with 3 parameters

So I have this dynamical system given by: $$ \left\{\begin{aligned} x' &= a(y-\phi(x))\\ y' &= x-y+z\\ z' &= -by \end{aligned}\right. $$ where $\phi(x) = \mu x^3 - \nu x$ and $a,b,\mu,\nu$ ...
5
votes
1answer
6k views

Can Mathematica help to rearrange an equation?

I have an equation which I want to rearrange so one variable is on the left side by itself. Normally, I would do this with basic algebra and identities. However, is it possible to use Mathematica to ...
5
votes
2answers
108 views

How to stop solving after one solution is found?

With[{P = 12}, Solve[{a^2 + b^2 == c^2, a + b + c == P, 0 < b <= a, c > 0}, {a, b, c}, Integers]] {{a -> 4, b -> 3, c -> 5}} ...
5
votes
2answers
162 views

How to find an estimated root efficiently with binary stochastic estimation of the function?

I am trying to construct certain nice geometric patterns, and one particular feature I'm trying to optimize is the balance of black and white corresponding to point being inside, or outside a ...
5
votes
1answer
155 views

Solving an identity

I am trying to solve this identity using Mathematica: $$3 \sqrt{5 \sqrt[3]{37}-16}=\sqrt[3]{a}-\sqrt[3]{b}-c $$ I know that a, ...
5
votes
1answer
93 views

Passing starting values to FindRoot when called by the Shooting Method of NDSolve

According to the documentation of the NDSolve Method -> "Shooting" option, it is possible to pass to ...
5
votes
1answer
315 views

Calculate max value

Update Using ...
5
votes
1answer
203 views

Find the domain in which a sequence converges

I seem to have a lot of problems in the format of: Given the sequence $a_n=\frac{n^\alpha+1}{2n^8+5}$, for $\alpha > 0$, find the values of $\alpha$ such that the sequence converges. However, ...
5
votes
1answer
732 views

Solving recurrence relation using Mathematica defined in a piecewise way

I have a recurrence relation defined as following: RSolve[ { p[0] == p0, p[1] == λ p[0]/μ, p[i + 1] == λ p[i]/(2 μ) }, p[i], i ] Note that the relation ...
5
votes
1answer
67 views

How to solve an equation system whose equations are based on Mod[x,y]

I'm trying to find the parameters of a specific linear congruential generator (pseudo random number generator) based on a list of the first five numbers generated by it. The list is this one: ...
5
votes
1answer
213 views

Finding and plotting roots with parameter

I have the following function: $$u(x,t)=\sum_{k=1}^\infty\frac{v_0 L}{\lambda_k c \sin\left(\lambda_k\right)}\sin\left(\frac{\lambda_kx}{L}\right)\sin\left(\frac{\lambda_kct}{L}\right)$$ and this ...
5
votes
1answer
393 views

Use Pattern with Vectors in Assumptions

I had a problem using Pattern with Vectors in Assumptions. Here's how I use ...
5
votes
1answer
227 views

Verify Harmonic Addition Theorem with Mathematica

I've calculated the harmonic addition theorem and would just like to check my result with Mathematica. However, ...
5
votes
2answers
656 views

Extracting Reduce results

I'm solving a Diophantine equation inside of a function using Reduce but I'm having trouble extracting the necessary parts of the answer. For example, if my input ...
5
votes
2answers
501 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
1answer
252 views

Incompleteness of FindInstance

I've been playing around with Resolve and FindInstance. I know that Resolve should always ...
5
votes
0answers
51 views

Is there a way to tell how much time FindRoot is spending calculating Jacobians? [closed]

I am solving a large number of trascendental equations by scanning FindRoot over a bunch of seeds and a bunch of different parameter values, and I would like to ...
5
votes
0answers
213 views

Strange behaviour of NotElement in Reduce

Bug introduced in 8.0 or earlier and persisting through 10.4 or later Confirmed, as bug, by WRI: [CASE:3527549] I have observed curious behaviour on domains supplied to ...
5
votes
0answers
165 views

VERY slow to solve a nonlinear equation involving matrix eigenvalues

I am trying to solve a problem similar to the following: ...
5
votes
0answers
97 views

Find regions in which the roots of a third degree polynomial are real

I have to find the roots of a third degree polynomial in $\phi$ that depends from 3 parameters, namely $t,s,w\in \mathbb R$. In order to do that I've used the command ...
5
votes
0answers
115 views

How to find an algebraic representation

I want to find matrix representation of some algebra (Clifford algebra Cl(3,1) in this case). Here is an example, which I would like to get and then extend to higher dimensional matrices. Suppose we ...
5
votes
0answers
74 views

Reduce doesn't solve this equations with condition

Reduce[{4==Abs[-1+2 Cos[2A]+2 Cos[2B]+2 Cos[2(A+B)]]},{A,B},Reals] worked well, but ...
5
votes
0answers
292 views

What is the algorithm behind Mathematicas Reduce in this equation?

Mathematica can use Reduce to arrive at an answer in this equation ...
5
votes
0answers
100 views

Using quantifiers with free indices as declarative conditions in equation solving

This is a toy example: Is there a way to rewrite Reduce[{Norm@{c[1], c[2], c[3]} == 1, 0 <= c[1] < c[2] < c[3]}, {c[1], c[2], c[3]}] ...
5
votes
0answers
2k 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
3answers
496 views

How to obtain a fixed number of solutions to $x+y=1$?

How to obtain a fixed number of solutions to $x+y=1$? For example, I want 4 solutions: $(1,0),(2,-1),(3,-2),(-2,3)$ I want to do this also on bigger equations such as $x+y+...=1$ and then I want to ...
4
votes
7answers
6k views

What is the equation to plot a vertical line?

Please help. Trying to find area between three curves, e^-x, x = 2, y = 1. Can't find out how to plot x = 2. Don't want to use <...