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

learn more… | top users | synonyms (5)

6
votes
0answers
222 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( \frac{...
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
0answers
84 views

Improve Accuracy of FindRoot

This is a follow up question to an earlier questions I asked. The link is: Product of $N\ 2 \times 2$ matrices and subsequently solving an equation dependent on the product As I am new here, I am not ...
4
votes
0answers
107 views

Why does Solve work in v9.0.1 but not in v10.3.1

I want to solve the following DE $$ P'=k(M-P)P, P[0]=53.08, P[50]=23.192, P[100]=76.212.$$ The Mathematica command is ...
4
votes
0answers
149 views

Complex boundary conditions and NDEigensystem

I am having difficulties implementing a Neumann value when numerically solving the Navier equation using NDEigensystem. The Navier equation is given by $\nabla^2 \vec u + (p^2 - 1) \nabla(\nabla\...
4
votes
0answers
79 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 ...
4
votes
0answers
127 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
75 views

Low accuracy while solving a multi-variable equation numerically

My question is as following: I was trying to make a function Vs2[r] with two parameters fit a set of conditions. Vs2[r]: <...
3
votes
0answers
135 views

How to solve partial differential equation

How to solve in Mathematica this partial differential equation: $0.5\frac{\partial t(x,y)}{\partial x}+1.5\frac{\partial t(x,y)}{\partial y}+t(x,y)=y\cdot \sqrt{1+x^{3}}$ with condition $t(1,y)=y+2$? ...
3
votes
0answers
173 views

solve tricky nonlinear ODE on a grid

I've been working on the same question for a while now and would appreciate any advice on what I'm doing wrong or any direction on how to actually solve the problem. :) On xzczd's advice, I retried ...
3
votes
0answers
95 views

Extending Bresenham's algorithm for subpixel rendering

Original Bresenham's line drawing algorithm it limited to drawing lines between two integer pixel positions and does not allow subpixel rendering. However with my reformulation of this algorithm as a ...
3
votes
0answers
128 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
124 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 ...
3
votes
0answers
154 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 ...
3
votes
0answers
124 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
86 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, j$....
3
votes
0answers
285 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 = \tilde{u}(x)...
3
votes
0answers
172 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
321 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
155 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 f+...
3
votes
0answers
129 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
767 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
166 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
87 views

How to define amazing Goodstein's sequence?

There is an amazing and counterintuitive theorem: For all $n$, there exists a $k$ such that the $k$-th term of the Goodstein sequence $G_k(n)=0$. In other words, every Goodstein sequence ...
2
votes
0answers
67 views

Solving an equation entered by the user via the Input function

I am new to Mathematica. I am trying to create a program that takes in an equation from user input and then gives the answer. I started with a manipulate function as I created a similar program for ...
2
votes
0answers
56 views

Unusual error message from Solve when called by DSolve

Question 111150 presents an equation nonlinear in its leading derivative ...
2
votes
0answers
145 views

Can equation solving be done in parallel?

I'm solving a set of equations that look like: ...
2
votes
0answers
99 views

Solving this system of nonlinear equations

Edit Interestingly, I realized solving my problem with $2n$ equations appeared to be more efficient than using mathematical results on the solutions to reduce the problem to $n$ equations (the $n$ ...
2
votes
0answers
99 views

Solving a non-linear equation in Mathematica with many variables

I have a beautiful equation, where I am trying to compute for R2. I am using Mathematica. So far every single time that I have tried to use Solve[] or Reduce [] the computation takes forever. I hope ...
2
votes
0answers
64 views

NIntegrate over statistical variables with parameter

Essentially, what I am trying to do is the following: Integrate a function which contains the parameter k ...
2
votes
0answers
48 views

speeding up Solve(ing) the Diophantine equation

Trying to solve the Diophantine equation (or variation?) using Solve in order to generate some chemical formulas. I need to generate all possible solutions. The ...
2
votes
0answers
121 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
104 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
217 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
140 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
416 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
173 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
488 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$ $(\frac{1}{x^2}+b)(1+x)^d+(1+...
2
votes
0answers
195 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
210 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
1k 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
57 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 ...