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

learn more… | top users | synonyms (5)

3
votes
1answer
70 views

Using result of Solve in further calculations

I'm trying to make a code for Gauss-Legendre quadrature. First I find the roots of the polyinomial with roots[n_,x_]:=Solve[LegendreP[n, x] == 0] Then I need to ...
3
votes
3answers
171 views

Reduce gives wrong solutions

Given the follwing system of equations eqs = { x/(x - 2) - 4/y == 1, x/(2 + 1/y) == y/(3 - 1/x) } Although the domain of this system does not contain ...
3
votes
1answer
104 views

Tree-like nesting of function

Is there an elegant way to nest a function in a similar fashion to tree graphs; i.e., create and follow branches until some criterion is met (until we reach a leaf)? An example: We want to ...
3
votes
1answer
205 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
1answer
129 views

Solving simultaneous systems of linear equations

I have several systems of linear equations in matrix form, and I would like to find the $\vec{x}$ that simultaneously solves them: $$A_1\vec{x} = b_1$$ $$\vdots$$ $$A_n\vec{x} = b_n$$ Is there a ...
3
votes
1answer
514 views

Mathematica takes too long to solve this

I have the following system of equations: ...
3
votes
1answer
212 views

Aborts and beeps without error warning

I have a piece of code that generates a figure, performing a calculation to generate each pixel. At a certain resolution it runs fine, but after a high enough resolution it fails. Mathematica beeps ...
3
votes
1answer
138 views

How to transform an equation so that the coefficient of a certain term becomes 1

Say I have an equation like A + 5 B + 10 C == 0 and I want to have Mathematica write this as (A/5) + B + 2 C == 0 i.e., I ...
3
votes
1answer
156 views

FindRoot: Minimum Execution Time

Is there a minimum execution time for the command FindRoot? If so, can I get rid of it? For example, if I compute the time to solve the fundtion x=0, it takes around 0.27ms per evaluation: ...
3
votes
1answer
90 views

Have transformation rules to be explicitly set?

I think this problem is not the hardest one, but I can't hit on it. I want to NSolve an equation. For this equation I use a transformation rule "ampTimes". First I used it like this: ...
3
votes
1answer
250 views

Using WhenEvent for derivative of discontinuous function

I have a discontinuous function ($u(t)$, a square wave) and I would like WhenEvent to trigger when the signal goes high/low, i.e. when the value of $u(t)$ changes. ...
3
votes
1answer
1k views

Solve pair of recurrence relations

[Corrected equations and added simple example] Can you solve a system of (loosely) coupled recurrence relations like this in Mathematica somehow? ...
3
votes
0answers
74 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
117 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
140 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) ...
3
votes
0answers
75 views

Equilibrium solution to system of 4 ODEs with 6 parameters [closed]

I have the following system of 4 ODEs of 4 state variables $(M_{U},M_{R},F_{N},F_{M})$ and 6 parameters $(\alpha,b,\gamma,\mu_{M},\delta,\mu_{F})$ that I am trying to find the equilibrium solution ...
3
votes
0answers
170 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
80 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
124 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
1answer
208 views

How to change parameters to equations of physical properties they belong to in an equation?

First off, I want to apologize for my poor English as I am not a native speaker. In my physics research, I had to solve a system of equations, from which I would get a solution for the speed of an ...
3
votes
0answers
121 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
138 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
120 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, ...
3
votes
0answers
265 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
171 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
317 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
151 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
128 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
744 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
165 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
3answers
413 views

Solve an equation that include Gamma

I want to solve the following equation Solve[Gamma[1 + x]/Gamma[x - 1/2] + 1 == 0, x] I have answer x=-0.25, but I can not obtain this answer with Mathematica. ...
2
votes
3answers
445 views

Solving an equation involving an integral

How do I get the following into Mathematica, solving for $a$: $$ 0.7 = 1 - \frac{2}{a} \times \left[ \frac{1}{a} \int_0^a \frac{x}{\exp(x)-1}\mathrm dx + \frac{a}{6} - 1\right] $$
2
votes
4answers
310 views

solve equation with loop

I want to solve equation with different variable and list answers in a table how can Solve it in a loop (Without repeating the equation) ...
2
votes
3answers
280 views

Minimization problem

How could you find the smallest constant $c$ that satisfies $$\frac{3^{3k}e\sqrt{3}}{\pi\sqrt{k}\;2^{3/2+2k}} \leq 2^{\;c\;k}$$ assuming $k\geq 1$ is an integer. You can of course take logs of ...
2
votes
2answers
235 views

Finding integral points on a surface

Suppose I have a dimension formula (for a Lie algebra representation) that is as follows: $$ d(a,b) = {(a+1)(b+1)(a+b+2) \over 2} $$ Now consider the surface $F(a,b,n) = 0 = d(a,b) -n$ where $n \in ...
2
votes
4answers
587 views

find the number of integral solutions a+b+c+d+e+f = 18 [duplicate]

Find the number of integral solutions of a + b + c + d + e + f = 18 where a, b, c, d, e, f are elements of the range ...
2
votes
2answers
2k views

Why doesn't Mathematica solve $x=\cos\,x$ properly? [duplicate]

Why doesn't Mathematica solve x == Cos[x] properly? Both Solve and NSolve fail with the ...
2
votes
2answers
139 views

Solving an equation and finding a couple of results given that they are integer

My problem is to find all solutions to an equation with 2 variables between two values that define the range of my search. e.g., $(x+y-y²)/(1-y-x²)=0$ with solutions ${x,y} \in \mathbb{Z}$.
2
votes
3answers
2k views

Solve and plot an equation that is dependent on a transcendental equation

I am trying to plot a function that uses the values of a transcendental equation as its input values: y = Tanh[(2 y + x)/2.5] z = 2.5 * Log[2*Cosh[(2 y + x)/2.5]] ...
2
votes
2answers
382 views

Get the approximate analytical solution to an algebra equation

Is there a way to get the analytical form to the root of this equation: 2 A e (-ArcTanh[smax/s0] + ArcTanh[(smax Cos[w ζ])/s0]) == F0 ...
2
votes
2answers
287 views

Check if integer has even number of digits

I want to check if a number has an even or odd number of digits. 10 would be even because it has 2 digits 300 would be odd ...
2
votes
1answer
2k views

Solving functional equations in Mathematica

I am trying to solve functional equations in Mathematica, but got nowhere. Just a simple example: $$f(x+y)^2=f(x)^2+f(y)^2$$ for all $x,y \in R$, assuming f is a real-valued function. How should I ...
2
votes
3answers
1k views

Finding roots of nontrivial equation

I am trying to solve the following nontrivial equation: ...
2
votes
2answers
234 views

Simultaneous equations—classify the solution set as finite, infinite, or none [duplicate]

How do I find the values of $b$ and $c$ for the equations $x + 5 y = 4$ and $2 x + by = c$ with a unique solution, an infinite set of solutions and no solutions with Mathematica?
2
votes
1answer
4k views

Forcing FindRoot to return only real solutions

FindRoot documentation reports that if the equation and the initial point are reals, the solutions are searched in the real domain. However, in the following case I ...
2
votes
3answers
227 views

Find roots of a peak to calculate its width at 10% peak height

I have numeric data points (x,y) values from a peak. I constructed a interpolation to get an interpolating function. Interpolation[data] What I want to do is the ...
2
votes
1answer
2k views

Mathematica Can't Evaluate This?

I am trying to find the inverse function of $\frac{\ln{x}}{x^2}$ in Mathematica. I am using the Solve[] function to find the inverse (yes, I know about ...
2
votes
2answers
201 views

Study the possible coincidences of the three clock hands hour, minute, and second

Inspired by the very nice problem solutions of the problem of coinciding clock hands, I'd like to extend the question including a "second" hand. While it is obvious that all three hand can coincide ...