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

Tracking FindRoot errors in NestWhileList and Table

I'm solving many systems of two equations by varying the initial parameters. I'm running FindRoot in a NestWhileList where the ...
3
votes
1answer
787 views

wrong use of _?NumberQ?

I am having trouble with a numerical exercise. Here's a simplified version of my problem. I define two functions to find the solution to two constrained maximization problems: ...
3
votes
1answer
81 views

Solve equation for a function as solution

I need to solve the following function: $$ h(z) = \sum^{40}_{k=0} \frac{k \, p(k)}{z- k \, h(z)} $$ where $p(k)$ is a probability distribution (in this case binomial with n = 40). Note that $h(z)$ is ...
3
votes
1answer
57 views

How can I choose the coeffcients so that the givent system of equations has four integer solutions?

Consider the system of equations with unknown are $(x,y)$. $$\begin{cases} a_1x^2+b_1y^2 + c_1xy+d_1x + e_1y+f_1=0,\\ a_2x^2+b_2y^2 + c_2xy+d_2x + e_2y+f_2=0. \end{cases}$$ I want to find the ...
3
votes
1answer
94 views

How to efficiently find only the rational roots of a rational complex polynomial?

I need to find all the rational roots of a bunch of high-order polynomials with rational complex coefficients. Is there an efficient way to do this? My best effort on an example polynomial: ...
3
votes
1answer
180 views

How to tell Solve[] to satisfy as many constraints as it can?

Is there any way to tell Solve (or something similar) to return values for variables that solved as many constraints as it happened to satisfy? Edit: I had ...
3
votes
1answer
67 views

Reduce is not giving answers, only verifies it's true

I'm trying to find a range for {q1,q2} with the conditions listed such that U1R >= U1S ...
3
votes
1answer
68 views

Finding inconsistencies in a system of equations

I am a rookie in Mathematica and apologize in advance. Here is my question: Are there any tools, built in the program, of detecting inconsistency in a nonlinear system of equations? Whilst finding ...
3
votes
1answer
98 views

Unable to solve equation for a variable by any method

I am looking to solve for the value of r when the derivative of the expression ...
3
votes
1answer
72 views

SolveAlways 'never' completes

After allowing Mathematica to run on ...
3
votes
1answer
105 views

Proof or find counterexample to inequality (in three vars)

I'm trying to show \begin{eqnarray*} \left(\frac{1-q^r}{1-q^n}\right)^{\frac{1}{n-r}} \left(\frac{q^n \log (q)}{\left(1-q^n\right) (n-r)}-\frac{\log \left(\frac{1-q^r}{1-q^n}\right)}{(n-r)^2}\right)&...
3
votes
1answer
118 views

Solve Differential Equations with Boundary Condition at Variable Point

I have a set of differential equations: ...
3
votes
2answers
191 views

How do I solve for a variable inside a sum but independent of the summation index?

How would I go about solving the following for c? Solve[0 == Sum[(t[i]*m[i] - c*t[i]^2)/s[i]^2, {i, 1, n}], c, Reals] I get ...
3
votes
1answer
108 views

How can I solve my easy system of equations?

I am new to Mathematica and am afraid this is a silly question. I have this very easy substitution, it has three equations on three unknowns {w1, l1, t1}. I want to ...
3
votes
2answers
171 views

Numerical root-finding: Find leftmost root in a domain

This is surely a duplicate, but I cannot seem to find a similar question. I use NDSolve to find an InterpolatingFunction y[t] ...
3
votes
1answer
315 views

Solving a large system of nonlinear equations

I have a system of 80 nonlinear equations with 64 variables to solve in Mathematica but it takes a lot of time (about two hours) and then it return "no more memory available" knowing that I have 8GB ...
3
votes
1answer
137 views

Solving a system of generated equations?

I would like to generate a function in the following form, where the number of terms can be specified arbitrarily: ...
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
366 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
311 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
178 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
179 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
105 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
234 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
131 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
549 views

Mathematica takes too long to solve this

I have the following system of equations: ...
3
votes
1answer
225 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
144 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
167 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
253 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
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
83 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
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
1answer
231 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
125 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 ...