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

learn more… | top users | synonyms (4)

7
votes
0answers
437 views

Creating and using an explicit piecewise function in a convenient way

I have a set of data points that define a function in the form curvePts = {data1,data2,...} where ...
7
votes
0answers
184 views

What Method options does FindInstance accept?

While answering Trouble getting FindInstance to return multiple results for certain constraints it became clear that for this problem FindInstance would be better ...
6
votes
0answers
127 views

FindRoot kills kernel

Bug introduced in 7.0.1 or earlier. Persists through at least 10.0.2 I'm having problems with FindRoot killing the kernel. In my real application, I want to FindRoot on a black-box function ...
5
votes
0answers
80 views

Using quantifiers with free indices as declarative conditions in equation solving

This is a toy example: Is there a way to rewrite ...
5
votes
0answers
982 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 ...
5
votes
0answers
86 views

Faster way to perform SameQ[Reduce[…], Reduce[…]]

I'm have some expressions that I need to confirm the equality of. SameQ[Reduce[...], Reduce[...]] works like a charm for the more simple expressions, however when ...
4
votes
0answers
213 views

What is the algorithm behind Mathematicas Reduce in this equation?

Mathematica can use Reduce to arrive at an answer in this equation ...
4
votes
0answers
214 views

How to get vertices of a polytope given by equalities and inequalities?

I have a polytope given by some equalities and inequalities, how can I get all the vertices of the polytope? If there's a way to deal with the one given by inequalities only, how can I reduce the ...
4
votes
0answers
148 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( ...
4
votes
0answers
114 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
225 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
106 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
89 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
430 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
147 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
997 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
98 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
108 views

Volterra integral equation

I have to find an approximate numerical solution for the equation $$ F(x) - \lambda \int\limits_1^{x} \text{d}s \;s^2 F(s) Z(x-s) = G(x) $$ $$Z(s) = (\psi''(1-2\ h\ i\ s)- 0.5 \psi''(1-2\ h\ i\ s))$$ ...
2
votes
0answers
64 views

How does Mathematica resolve symbolic systems of inequalities?

Before getting into my (probably annoying) question, let me provide some context: I recently started using Mathematica to automatically perform stability analyses that I previously did by hand. A ...
2
votes
0answers
118 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$ ...
2
votes
0answers
83 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
74 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
61 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, ...
2
votes
0answers
83 views

Quartic equation

After using Solve I am choosing the physical root of my equation and dealing further with this root. However, when I am changing parameters using Manipulate, I clearly see not smooth behaviour. Is ...
2
votes
0answers
82 views

How to handle infeasible points in FindRoot?

I am calling FindRoot[f[x,y],{{x,xInit,xMin,xMax},{y,yInit,yMin,yMax}}] where for some points {x,y}, ...
2
votes
0answers
277 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
51 views

Equation involving hypergeometric functions

I want to solve this equation but nor NSolve nor Solve are able to do this. ...
2
votes
0answers
116 views

Question about Solve

I don't understand what Mathematica is doing to get these solutions. I ask it to solve two equations, the left hand side of both expressions are the same, the right hand sides differ by a ± sign. I ...
2
votes
0answers
109 views

Equation Solving: Solution in terms of predefined functions

Is there a way that I can have Mathematica give me the solution of a problem it has solved in terms of functions that I have predefined? The specific problem I'm trying to solve is: ...
2
votes
0answers
167 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. ...
2
votes
0answers
186 views

Problems with Simultaneous Equation Solving

Let me first explain what I am doing and then I will explain the problem that I can't figure out. I am basically trying to solve 2 equations for 2 unknowns (a1 and ...
2
votes
0answers
417 views

I can't solve this set of equations

I am having trouble solving a set of equations in mathematica. I do not get any errors, but i don't get any answers as well. I have some restrictions on my parameters and variables in the equations ...
2
votes
0answers
133 views

Nest for large value of n

I am working on a problem to find the Feigenbaum constant. In this problem I have to get $f^{2^n}(1/2)$ with $f(x) = kx(1-x)$. I use Nest to get $f^{2^n}(1/2)$ but ...
2
votes
0answers
285 views

Conditional statements in intial conditions?

This is potentially a daft question, but I thought I'd ask it; I have some material free to diffuse in a boundary between rn and ro; I've been able to get it working nicely for neumann type boundary ...
2
votes
0answers
310 views

differential equations with implicit functions

I have an ordinary differential equation that contains coefficient functions that depend implicitly on the independent variable via an algebraic equation. I am trying to go ahead and use NSolve to ...
2
votes
0answers
149 views

Why is FindRoot initial value far from the specified one?

I am trying to numerically find the root of a function that looks a bit like: 1/x - (SchurDecomposition[A[x]][[2]])[[1]], where ...
2
votes
0answers
182 views

Solving $k\log(n) - k\log(k) \leq u$

How can I persuade Mathematica to help me find solutions to the inequality $k\ln(n) - k\ln(k) \leq u$ over the reals, assuming $u \geq n \geq k \geq 1$. When I evaluate ...
2
votes
0answers
207 views

Am I entering this 2-variable recurrence correctly?

I'm trying to solve this recurrence with two variables. ...
2
votes
0answers
274 views

FindRoot - Speed and errors

I am using FindRoot[] to solve a complicated equation. It seems I get the correct answer even though I get errors about an ...
2
votes
0answers
400 views

Numerically/Analytically Solving a System of Equations

I have $6$ functions $f_i(x,y,z)$, $(i = 1, \ldots, 6)$ in three variables $x,y,z$, and I would like to find a simultaneous instance of these variables, say $(x_0, y_0, z_0)$, such that $f_i(x_0, y_0, ...
1
vote
0answers
72 views

Equation Solving

As I stated in this question I want to solve a series of equations while setting some rules such as q[2]->1 or q[2]->3*10^-4. However, what should I do if I want to set another rules to my input? I ...
1
vote
0answers
42 views

Apart for complex roots?

Let p[x] be a polynomial in x and consider the partial fraction decomposition of 1/p[x]. ...
1
vote
0answers
33 views

Resolve not accepting equations involving factorial

If I take , Resolve[ForAll[x,x!>=0],Integers] True But if I make slight modification, ...
1
vote
0answers
58 views

Help inteprete the solution given by reduce

i am trying to solve for the a2 in the below equation on Mathematica ...
1
vote
0answers
75 views

Using NMinimize to solve the Multi Commodity Flow problem

I have the following optimization problem and I was hoping to solve it with NMinimize. Unfortunately, I don't get the desired result. I would like to get an optimal vector for $d_{e}$. This is how ...
1
vote
0answers
45 views

Creating a matrix from a list of simultaneous equations

I have a list of n simultaneous equations of the form, w A[n] = (K[n,n+1]+K[n-1,n])A[n] - K[n,n+1]A[n+1] - K[n-1,n]A[n-1] ...
1
vote
0answers
80 views

Comparing two ways of fitting systems of ODEs to data

I'd like to ask you if you see any difference between the two following ways of estimating parameters. Given my system of ordinary differential equations that I want to use to find the values of all ...
1
vote
0answers
72 views

Stability region of midpoint method using Solve and RegionPlot

I've been trying to obtain the stability region in the $(\lambda_R, i\lambda_I)$ plane of the midpoint rule $$y_{n+1}= y_{n-1}+2h\dot{y}_n,$$ for the model problem $\dot{y}=\lambda y$. The roots for ...
1
vote
0answers
119 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 = ...
1
vote
0answers
82 views

Code for solving numerically an integro-differential equation

First of all, I want to greet the community. This is my first question, but I hope I will be able to help answering others members questions, although I am quite new working with Mathematica. I would ...