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
77 views

Finding linearly independent combination of vectors with a certain symmetry property

I have two linearly independent vectors $v_1$ and $v_2$. Is there a way for me using Mathematica to find two linear combinations $w_1$ and $w_2$ such that: $w_1$ and $w_2$ are linearly independent. ...
3
votes
1answer
102 views
3
votes
1answer
293 views

Minimizing a potential using Solve

I am having problems minimizing a potential: $\text{V}(h,\eta)=\gamma \left(-h^2\right) \left(\eta ^2 \cos ^2(\theta )+\eta \cos (\delta ) \sqrt{-\eta ^2-h^2+1} \sin (2 \theta )+\left(-\eta ^2-h^2+...
3
votes
1answer
362 views

ParametricNDSolve and Constraint satisfaction problem

I am trying to find the values of two parameters that allow for a specific result of my differential equations, under different initial conditions. Given a system of equations: ...
3
votes
1answer
2k views

NDSolve does not respond

For some sets of constants, NDSolve gives me true solutions, but when I try for example, T = 1/(2*2200), Mathematica does not respond. What can I do? The code below ...
3
votes
1answer
310 views

Result of Reduce contains duplicate terms

On my system, the result of this: f[x_] := 2 Sin[x] + Sin[x]^2; Reduce[f'[x] == 0, x] contains the following expression twice: ...
3
votes
1answer
167 views

Solve problem: Inconsistent or redundant transcendental equation

I would like to mechanically Solve[a y == (1 - a) x && u0 == x^a y^(1 - a),{x,y}], given 1>a>0 and u0>0. The emphasis is on mechanically: I know it is ...
3
votes
1answer
76 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
176 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
225 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
130 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 quick/...
3
votes
1answer
541 views

Mathematica takes too long to solve this

I have the following system of equations: ...
3
votes
1answer
218 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
143 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
162 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
91 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
252 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
126 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
145 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\...
3
votes
0answers
80 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 for:...
3
votes
0answers
172 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
86 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
126 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
215 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
148 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
123 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
272 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
318 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
154 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
755 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
3answers
421 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
450 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
339 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
597 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
394 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
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
141 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
3k 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
3answers
133 views

Solving my graduation exam geometry problems with Mathematica [on hold]

How capable is Mathematica when it comes to typical school geometry problems? Is there a way to throw constraints at it and get it to spit out a solution? Here are two geometry problems from my ...