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
2answers
290 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
380 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
889 views

Methods for NSolve

I've been using NSolve a moderate (50-100 equations) size system of linear equations and it has been working splendidly (Solve ...
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
531 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
3answers
120 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
185 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
906 views
5
votes
1answer
287 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
523 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
4answers
104 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
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
107 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
161 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
92 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
728 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
212 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
225 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
645 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
499 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? [on hold]

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
161 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
96 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
113 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
291 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
495 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 <...
4
votes
2answers
1k views

Can mathematica solve this equation?

I have the following Mathematica code to solve an equation coming from isoparameteric finite element for quadratic tetrahedron element. I need to get the iteroplation fast and accurate. I only find ...
4
votes
3answers
2k views

Using DSolve to solve for x[t,parameter]

I'm trying to solve an ODE with two independent variables (a cannon firing from a cliff incorporating wind resistance dependent on velocity). I've tried the following for the x-component: ...
4
votes
3answers
567 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,\...
4
votes
2answers
1k views

How can I solve the equation with integers as a solution?

I want to solve the equation $$(x-1)^2 + (y-1)^2 + (z-1)^2 = 49$$ where $x$, $y$, $z$ are integer and $x \neq 1$, $y \neq 1$, $x \neq 1$. How do I tell Mathematica to do that?
4
votes
3answers
4k views

How to solve this system of equations?

I want to solve the following system of equations: $$\begin{cases} 1+\sqrt{2 x+y+1}=4 (2 x+y)^2+\sqrt{6 x+3 y},\\ (x+1) \sqrt{2 x^2-x+4}+8 x^2+4 x y=4. \end{cases}$$ I tried ...
4
votes
2answers
726 views

How do I solve N simultaneous equations for N variables?

I have a function: f[x_] := x + 31 x^3 + 5 x^25 Which I want to find an expansion for: ...
4
votes
3answers
513 views

Finding parameters making real part of eigenvalues vanish

I have the following $\;3\times3$ matrix: $\left( \begin{array}{ccc} 0.04 -0.4 b & 0 & 0.04 -0.4 b \\ 0 & -0.08-1.2 b & -0.06-0.9 b \\ 1.04 -0.4 b & 2.08 -0.8 b & 0 \end{...
4
votes
5answers
378 views

Solve for exponents

I have this problem that I"m trying to work out: Here's my attemp: me = 0.51109991 mp = 938.2723 mpe = mp/me Solve[(2^n)*(3^mm)*(\[Pi]^p) == mpe, {n, mm, p}] ...
4
votes
5answers
358 views

On finding all the positive integral solutions of $x^2+y^2=z^2+1$

I am a new to Mathematica. My goal is to find many (if not all) positive integer solutions to the equation: $x^2+y^2=z^2+1$ using Mathematica. However the problem is that I can only find a ...
4
votes
1answer
2k views

Solve system of equations related to perspective projection

I'm trying to solve system of equations. I have 4 points before perspective projection and after and I whant to obtain matrix of perspective projection $m_{i,j}$ ($x,y,u,v$ are known , $m_{i,j}$ ...