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

learn more… | top users | synonyms (5)

7
votes
0answers
115 views

How does Mathematica resolve symbolic systems of inequalities?

Before getting into my question, let me provide some context: I recently started using Mathematica to automatically perform stability analyses that I previously did by hand. A major part of the ...
6
votes
0answers
64 views

Choosing appropriate WorkingPrecision when solving numerical system of equations

Consider a linear system of equations, which is conveniently written as A.x=y. The matrix A has dimensions ...
5
votes
0answers
256 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
98 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
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 ...
5
votes
0answers
97 views

Faster way to perform SameQ[Reduce[…], Reduce[…]]

I'm have some expressions that I need to confirm the equality of. SameQ[Reduce[...], Reduce[...]] works like a charm for the more simple expressions, however when ...
4
votes
0answers
54 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 ...
4
votes
0answers
101 views

Synchrotron Radiation and ListDensityPlot

I'm trying to reproduce this plots: For this task i'm following this blog entry: Synchrotron radiation As you can see there, the task is to make a DensityPlot of a vector field essentially. This ...
4
votes
0answers
170 views

Solving a nasty equation

Here's an equation that I need to solve for $t_k$ $$ \theta_0 = \tan^{-1}\left(\frac{2(w_x - m)}{p - g {t_k}^2 - 2 t_k w_y}\right) - {\frac{1} {r}} \left( \left(\frac{m}{t_k}\right)^2 + \left( ...
4
votes
0answers
120 views

Symbolic solution to ODE, pure InverseFunction not evaluated

I'm new to Mathematica and I don't understand why in the solution of the following ODE, the #1 in the pure function is not immediately replaced by the corresponding ...
3
votes
0answers
64 views

NSolve doesn't work on an equation containing a numerical integral and constraints

I'm having trouble getting Mathematica to solve equations numerically. I know that its important to specify the type of variables for pattern testing (see e.g. here) but this doesn't seem to work. ...
3
votes
0answers
53 views

Method of characteristics, recursion formulas

I'm helping a friend out with a problem. He needs to plot out the solution to a water surge in a pipeline. I have all the parameters set and I've put the equations in Mathematica, but I can't get the ...
3
votes
0answers
71 views

Solving for the roots of a trilinear system of polynomials

I have been trying to solve for the roots of the following system of trilinear polynomials: ...
3
votes
0answers
91 views

Stability analysis of transcendental equation (stability crossing curves)

I am working with a non-linear delay system with three scalar delays. After taking the Laplace transform of the linearized system, the characteristic function is a transcendental equation with three ...
3
votes
0answers
125 views

Finding all zeroes - mysterious missing roots

Finding all the roots of a function of one real variable has been discussed extensively in e.g. Find all roots of an interpolating function (solution to a differential equation) (answer no. 2 by Jens ...
3
votes
0answers
196 views

Nonlinear FEM and FindRoot

I'm trying to develop a kind of nonlinear FEM application using mathematica to solve a bvp like the following: $$ \gamma(u') ~u^{iv} + 2 \gamma'(u') u''' u''+ u''^3 = f(x) $$ where $u = ...
3
votes
0answers
148 views

Solving Thue equations

A Thue equation is a 2-variable homogeneous integer polynomial of degree at least 3. It's well-known that such equations have only finitely many solutions over the integers. I'm trying to solve some ...
3
votes
0answers
285 views

Function with a sharp resonance which Mathematica fails to integrate. Why?

I'm a new user of Mathematica and I'm trying to use it to calculate the collisional cross-section as a function of energy for a given potential and decay rate. I know that the resulting function ...
3
votes
0answers
123 views

Solving huge equations

I have this expression: $d=-\frac{\sqrt{\frac{780 (f+130) \left(5 d^2 (13 c+f+130)-5 d (c (f+1300)-325 (f+160))+125 \left(2 (2 c+725) f+5 f^2+104000\right)+13 d^3\right) \left(5 d^2 (52 c+17 ...
3
votes
0answers
110 views

Deciding whether a semi-algebraic set is empty. Which command is best?

[Update] I was able to prove that the "larger" set $ 1\ge x \ge y,w \ge 0$ and $p(x,y,w)>0 $ is empty with another computer algebra system in my slow laptop in few minutes but in Mathematica even ...
3
votes
0answers
620 views

FindRoot gives a wrong solution which obviously should not be there

I got stuck on FindRoot and I didn't see any similar problem posted, so let me explain what I am trying to do and what problem I meet here. I try to find roots of a particular function, which in the ...
3
votes
0answers
153 views

RSolve not returning known solution

I have the following recurrence relationship for which I know the solution is given by a binomial distribution: ...
3
votes
0answers
1k views

Solving a system of nonlinear equations self-consistently

I am trying to solve a set of three non-linear equations in Mathematica. I need help with them. The Mathematica code (in plain text format) is attached below. If you copy & paste the code below ...
2
votes
0answers
74 views

Suggestions for solving a large linear system

I have a system of $\approx 200,000$ linear equations in $\approx 40,000$ variables (with rational coefficients) and I would like to determine the dimension of the solution space, which I know to be ...
2
votes
0answers
58 views

Is FindRoot wrong about its WorkingPrecision?

I noticed an odd mismatch between the behaviors of NDSolve and FindRoot. If you give NDSolve an equation of a certain precision, say 20, and tell it to use a higher WorkingPrecision, say 30, it will ...
2
votes
0answers
54 views

Problem with plotting a function that calls FindRoot

I am having difficulties in plotting the solution of FindRoot to solve 4 variables. I defined a function and then tried tp plot it, but got an error message ...
2
votes
0answers
67 views

How to handle solution returned by ParametricNDSolveValue in FindRoot

i was trying to solve a problem relative to root finding in a system of differential equations; here's a simpler case. With two distinct set of parametric differential equations, i'm able to find the ...
2
votes
0answers
76 views

Simplify large fractions

Hi i'm quite new with Mathematica and i'm experiencing some problems with Simplify and FullSimplify. I'm dealing with a sistem ...
2
votes
0answers
105 views

NDSolve is running an extremely long time: how can I save the existing data?

I am trying to solve a PDE by the following code. It takes 1 hour or so to reach t=25.72 but about 20 hours to reach t=25.72404031638060174049337306126853310997. Actually, the time step is extremely ...
2
votes
0answers
130 views

Proving large inequality in 6 variables

I have arrived at a very large inequality in 6 variables that I would like to prove is true. The command I have tried is the following: ...
2
votes
0answers
159 views

Apart for complex roots?

Let p[x] be a polynomial in x and consider the partial fraction decomposition of 1/p[x]. ...
2
votes
0answers
145 views

What are the steps using boolean algebra for transforming one equation into another?

I have a special case where I have two equations, which I have solved in a better and quicker way, but I need to solve with Boolean algebra. Using the laws of Boolean algebra, I have to prove that ...
2
votes
0answers
290 views

Solving a system of two nonlinear equations with many symbolic parameters

I have two nonlinear equations with two unknowns, $x$ and $y$, and many symbolic parameters, $a, b, c, d$. They look someting like: $\frac{1}{x}(1+x)^a (1-b)(cx+dy)=0$ ...
2
votes
0answers
127 views

Solving a system of difference equations in an exponential steady-state growth setting

I have a system of difference equations in an exponential steady-state balanced growth setting; all variables grow at the same steady-state rate, say, $g$. I use the notations as follows: Variable $X$ ...
2
votes
0answers
139 views

Error when using EquationTrekker

I just got to know EquationTrakker because I need to do phase graphs for my classes. When I load it ...
2
votes
0answers
71 views

How to solve for exponents in an infinite product?

Say I have a power series, $f(x,y) \in 1 + x\mathbb{Z}[[x,y]]$ Then it admits a unique expansion $$f(x,y) = \prod_{i=1}^\infty \prod_{j=0}^\infty (1 - x^i y^j)^{n_{i,j}}$$ for some integers, $i, ...
2
votes
0answers
103 views

How to handle infeasible points in FindRoot?

I am calling FindRoot[f[x,y],{{x,xInit,xMin,xMax},{y,yInit,yMin,yMax}}] where for some points {x,y}, ...
2
votes
0answers
641 views

FindRoot in range

I am trying to find all roots in a range, using y /. FindRoot[Sin[y] == 0, {y, #}] & /@ Range[7] which gives ...
2
votes
0answers
44 views

How can I select two lists from a list of lists so that I receive the system of equations has only rational solutions?

I want to make a sytem of equations has the form $$\begin{cases} a_1x^2+b_1y^2 + c_1xy+d_1x + e_1y+3=0,\\ a_2x^2+b_2y^2 + c_2xy+d_2x + e_2y+3=0. \end{cases}$$ so that the given system of equations has ...
2
votes
0answers
120 views

Question about Solve

I don't understand what Mathematica is doing to get these solutions. I ask it to solve two equations, the left hand side of both expressions are the same, the right hand sides differ by a ± sign. I ...
2
votes
0answers
131 views

Equation Solving: Solution in terms of predefined functions

Is there a way that I can have Mathematica give me the solution of a problem it has solved in terms of functions that I have predefined? The specific problem I'm trying to solve is: ...
2
votes
0answers
207 views

Problems with Simultaneous Equation Solving

Let me first explain what I am doing and then I will explain the problem that I can't figure out. I am basically trying to solve 2 equations for 2 unknowns (a1 and ...
2
votes
0answers
569 views

I can't solve this set of equations

I am having trouble solving a set of equations in mathematica. I do not get any errors, but i don't get any answers as well. I have some restrictions on my parameters and variables in the equations ...
2
votes
0answers
138 views

Nest for large value of n

I am working on a problem to find the Feigenbaum constant. In this problem I have to get $f^{2^n}(1/2)$ with $f(x) = kx(1-x)$. I use Nest to get $f^{2^n}(1/2)$ but ...
2
votes
0answers
349 views

Conditional statements in intial conditions?

This is potentially a daft question, but I thought I'd ask it; I have some material free to diffuse in a boundary between rn and ro; I've been able to get it working nicely for neumann type boundary ...
2
votes
0answers
395 views

differential equations with implicit functions

I have an ordinary differential equation that contains coefficient functions that depend implicitly on the independent variable via an algebraic equation. I am trying to go ahead and use NSolve to ...
2
votes
0answers
171 views

Why is FindRoot initial value far from the specified one?

I am trying to numerically find the root of a function that looks a bit like: 1/x - (SchurDecomposition[A[x]][[2]])[[1]], where ...
2
votes
0answers
195 views

Solving $k\log(n) - k\log(k) \leq u$

How can I persuade Mathematica to help me find solutions to the inequality $k\ln(n) - k\ln(k) \leq u$ over the reals, assuming $u \geq n \geq k \geq 1$. When I evaluate ...
2
votes
0answers
241 views

Am I entering this 2-variable recurrence correctly?

I'm trying to solve this recurrence with two variables. ...
2
votes
0answers
331 views

FindRoot - Speed and errors

I am using FindRoot[] to solve a complicated equation. It seems I get the correct answer even though I get errors about an ...