Questions tagged [equation-solving]

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

626 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
13
votes
0answers
119 views

Reduce[Sin[x]==0,x] splitting odd and even multiples of Pi

as from Title, I'm wondering why Mathematica does not automatically give a single, simple solution for the problem Reduce[Sin[x] == 0, x] for which I obtain the ...
7
votes
0answers
921 views

Does NSolve find all solutions?

Is the solution set returned by NSolve usually complete? Can I assume that there are no more solutions than what it returns? Consider systems of equations e.g. ...
7
votes
0answers
255 views

Strange behaviour of NotElement in Reduce

Bug introduced in 8.0 or earlier and fixed in 11.0.0. Confirmed, as bug, by WRI: [CASE:3527549] I have observed curious behaviour on domains supplied to NotElement...
6
votes
0answers
52 views

Limitations in using Reduce for Integer systems

The System For example, when solving the following system using Reduce in Integers: ...
6
votes
0answers
94 views

When does Root have a third argument

In Mathematica "11.0.1 for Microsoft Windows (64-bit) (September 20, 2016)", Root[#^4 + 1 &, 2]; actually has three arguments, as can be seen from ...
6
votes
0answers
177 views

Spurious DSolve Solution

Bug introduced in 8.0.4 or earlier, persisting through 12.0. DSolve quickly returns solutions to the following PDE (which is the homogeneous portion of the PDE in ...
6
votes
0answers
482 views

What is the algorithm used by Reduce for this equation?

Mathematica can use Reduce to arrive at an answer in the following equation: ...
6
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 ...
6
votes
0answers
344 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
364 views

How does NSolve find the roots to an equation?

I'm using NSolve to find the roots to a transcendental equation. For this, I need to specify a bounded region where roots are searched. What I want to know is how <...
5
votes
0answers
767 views

Solving a system of differential algebraic equations (DAE)

I am trying to solve a system of 8 differential algebraic equations, where equations 3 and 5 are differential equations and the rest are constraints which need to be satisfied. Also I only know the ...
5
votes
0answers
164 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
292 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
278 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 ...
5
votes
0answers
124 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
516 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 ...
5
votes
0answers
111 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]}] ...
4
votes
0answers
48 views

Force Solve to respect unit system

I give Solve quantities in FPS and it returns resulting quantities in SI, necessitating that I add an entire section to my workflow where I pick apart the solution ...
4
votes
0answers
62 views

Indicate Unit Vectors in Solve

I'm solving a system defined by two rays. The rays originate from two points {x,y,z} and {x,y,z-h}. They go in directions (unit ...
4
votes
0answers
128 views

DSolve, NDSolve with WhenEvent Give Incorrect Solution for Simple ODE

NDSolve Results On the course of addressing question 181974, I encountered the following problem. ...
4
votes
0answers
174 views

Discordance in plots

Let us solve the system with Floors ...
4
votes
0answers
53 views

A simple expression on which `Reduce` fails

The conditions in x ∉ Integers && y ∈ Integers && y == -(1/2) + x/2 are inconsistent but Reduce does not solve ...
4
votes
0answers
228 views

Why does Reduce[x + 1/x == 1/x, x] give the result x == 0, instead of the result False?

Why does Reduce[x + 1/x == 1/x, x] give the result x == 0 instead of the result False? There is no value of x that causes ...
4
votes
0answers
108 views

Order of evaluation of FunctionDomain, Reduce, Solve, etc

For the input (mine is much more complicated, this is a toy example): x/x /. x -> 0 the output of 1 (as opposed to ...
4
votes
0answers
114 views

How to extend a root finding method from 1D to multi dimensions

The number (0.25 Pi) is a third order root of f[x] below. Clear[f]; f[x_] := Sin[x] - Cos[x] - Sqrt[2] (x - 0.25 Pi); Chop@Normal@Series[f[x], {x, 0.25 Pi, 4}] -...
4
votes
0answers
396 views

Homotopy Continuation solution of system of polynomials

I have very large systems (>20) of polynomial (max degree 3) equations that I would like to find a solution to. I'm not interested in all solutions as presumably there are too many (a huge number ...
4
votes
0answers
121 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
355 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
372 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 ...
4
votes
0answers
324 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 ...
4
votes
0answers
270 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 ...
4
votes
0answers
408 views

Am I entering this 2-variable recurrence correctly?

I'm trying to solve this recurrence with two variables. RSolve[g[d, n] == 1 + g[d - 1, n] + g[d - 1, n - 1] && g[d, 0] == g[0, n] == 0, g[d, n], {d, n}] ...
4
votes
0answers
161 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
60 views

Different output to Solve command in Mathematica changing a variable for a subscripted variable

I have a system of equations that is solved with the Solve command in Mathematica 10.2, when one of the variables is changed for a subscripted variable the answer is not the same, the difference is ...
3
votes
0answers
122 views

Errors Using FindFit and FindRoot

I'm trying to fit the following equation to data: where and x and y are the dimensionless variables x=V/Vr and y=Ir/(2*I0). n=1 and I'm trying to find Vr and R. The blue is an example of the curve ...
3
votes
0answers
138 views

Unexpected Behavior of Parametric Sensitivity in ParametricNDSolveValue

Bug introduced in 10.4 or earlier and continuing through 11.3 Submitted as CASE:3916971 While exploring alternative methods of solving 33538, I encountered difficulties with the parametric ...
3
votes
0answers
133 views

Solving Monster Equations with Mathematica

I would like to solve the following equations in Mathematica. Here are the formulas: $u=0.5$ $\alpha=0.5$ $\epsilon_0=0.1$ $\epsilon_1=0.1$ $f_0(y)=\mbox{Gaussian density with mean -1 and ...
3
votes
0answers
150 views

Replace all (/.) parameters with equation or condition rather than single values

Suppose that one wants to calculate the surface charge density of a polarized dielectric ellipsoid, prolate ellipsoid with long and short axes as $a$ and $c$ where $a$ is parallel to the $z$ axis. $a$ ...
3
votes
0answers
436 views

Solve PDE over time-dependent region

I have some experience with C++ programming, but I am quite new to Mathematica and I think the latter requires a different mindset which I am still not used to. I would like to simulate the effect ...
3
votes
0answers
97 views

Does Reduce have problems with greek symbols that lead to wrong results?

I call Reduce on the following simple system of inequalities ...
3
votes
0answers
100 views

Spheroidal Eigenfrequency

I am trying to create a Mathematica function which will calculate the eigenfequencies of a spheroidal cavity. The expressions I use are based both on the book "Spheroidal Wave functions in ...
3
votes
0answers
81 views

Generating All Regular Multigraphs — Issue with Solve and/or FindInstance

I'm looking to generate all regular multi-graphs over $n$ vertexes and of degree $d$ in mathematica (side note: $n,d$ are fairly small, this problem gets unmanagable very fast). So here $n,d$ are my ...
3
votes
0answers
105 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
320 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
191 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$ ...
3
votes
0answers
294 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
451 views

Numerical solution to integral equation

I have an integral equation where one factor $f(t)$ in the integrand is defined in terms of an integral equation where it is again a factor. $\quad \quad y[t]=\int _0^tf[\tau ]g[t-\tau ]d\tau + h[t]$...
3
votes
0answers
272 views

How can I improve the solution of my PDE?

I want to improve the last code to find a better solution for umi. For now I'm using only MaxStepFraction. Also, when I put no ...
3
votes
0answers
1k 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+...
3
votes
0answers
524 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)...

1
2 3 4 5
13