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

learn more… | top users | synonyms (5)

4
votes
1answer
154 views

Does `Reduce` give complete and necessary conditions?

I want to obtain a solution for a system of five inequalities. For instance, ...
1
vote
1answer
74 views

using FindRoot with a Module that takes a parameter

Observed an error when calling a Module from FindRoot. I created a module which first solves a differential equation and then returns the value of that differential equation when the independent ...
1
vote
1answer
82 views

Is it possible to solve a certain set of nonlinear equations in Mathematica?

I read this question in a Chinese forum last year. At the end of the discussion, they provide a solution of the equations from a C++ program. Unfortunately, they did not get the right answer. ...
4
votes
1answer
163 views

What are the hidden specifications for FindRoot

The Help page of FindRoot says: "by default, FindRoot uses Newton's method (Newton-Raphson) to solve a nonlinear system" (or a nonlinear equation I ...
2
votes
1answer
56 views

Plot directional field with a initial vlue

I do not have enough reputation to comment on other questions, hence I have to ask my own one I want to solve the following nonlinear first order ODE, $$(1 + x) (0.5 y[x]^{-0.5} (x - y[x])^{0.5} - ...
1
vote
0answers
104 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 ...
2
votes
2answers
174 views

How can Mathematica help me to find a real radical expression for roots of this polynomial?‎

The polynomial $P(x)=x^‎‏4‏‎-‎‏4‏x^‎‏2‏‎-‎‏2‏x+‎‏1$‏‎ has ‎‏4‏‎ real roots (this can be clearly checked by plotting). But solving $P(x)=‎‏0‏‎$, using ...
0
votes
0answers
43 views

Finding common roots to several functions

I am trying to solve the following problem. Given $a\in\mathbb R^n$, $u\in\mathbb{R}^n$, $m\in\mathbb{N}^\star$, Find the/some common roots $(t_1,...,t_m)$ of the $\frac{m(m-1)}{2}$ ...
1
vote
1answer
115 views
49
votes
8answers
2k views

Updating Wagon's FindAllCrossings2D[] function

Stan Wagon's Mathematica in Action (second edition; I haven't read the third edition and I'm hoping to eventually see it), demonstrates a nifty function called ...
4
votes
2answers
147 views

how to find numerical zeros of $f(z)$ in a compact complex domain

I encountered a function $f(z)$ with several zeros in a region $1<\Re(z)<2, 3<\Im(z)<4$. These zeros are visible in a three dimensional plot like ...
2
votes
2answers
99 views

A peculiar behavior of NSolve[]

NSolve[Power[x + 1, x + 1] - Power[x, x + 2] == 0, x, Reals] returns: {{x -> 3.14104}} Great. However, ...
0
votes
1answer
59 views

Symbolically solve system of equations involving matrix

I need to solve the following system of equations symbolically for a_i, b_i, and c_i: ...
1
vote
0answers
89 views

Can Mathematica Find Closed-Form Solution this Equation

I have tries to solve an equation by using Solve[] to find pl as a function of other variables. However, it generates an Error! Can anyone help? ...
18
votes
5answers
16k views

Assign the results from a Solve to variable(s)

I understand Mathematica can't assign the results of a Solve to the unknowns because there may be more than 1 solution. How can I assign the 4 values of following result to variables?
0
votes
1answer
57 views

Question about the third argument of Solve [duplicate]

I'm learning about derivatives. This is my code: ...
2
votes
2answers
110 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 ...
19
votes
4answers
868 views

Hightlight all the self-intersections of a Lissajous figure

This graph–also known as a Lissajous figure–contains so many self-intersections.How can I highlight them? ...
3
votes
2answers
325 views

How to solve numerically an equation which contains numerical intergration?

F1[x_] := NIntegrate[E^Sin[y x], {y, 0, 1}] F2[x_] := Log[x] NSolve[F1[x] == F2[x], x] I want to solve the equation, but Mathematica don't solve it. ...
1
vote
1answer
106 views

Solving a complex logarithmic equation

I am trying to solve the following equation for $x$, with integer $n$. (1 - ((-1 + x)^E x^(1 - x) Log[x]^-E)^(1/(-1 + E)))/(-1 + x) == n
1
vote
1answer
332 views

Usage of Piecewise/RSolve to solve an equation based on condition

I am trying to compute price of call option on a stock using the finite difference method. I have a VBA code for it, and I am trying to find out how I can use Mathematica memoization functions and ...
2
votes
3answers
143 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
75 views

Reduce / Solve cannot solve simple equation

Why does Reduce[{a^b*c^b == c, 0 < b < 1, a > 0, c > 0}, c, Reals] give ...
0
votes
1answer
43 views

Solve and NSolve hitting iteration limit for simple system [closed]

I'm am trying to get Mathematica to determine the intersection of a line and a parabola, like so: Solve[{x + y == 6, x^2 == y}, {x, y}] But evaluating this gives ...
0
votes
1answer
52 views

Question about a catalytic equality?

I have a question on math.SE that explains how I came up with these equations. ...
1
vote
2answers
225 views

How to Fit experimental data with 4D ParametricNDSolve?

Let me first explain the background problem. I'm trying to solve a problem with a projectile motion in 4D such that the position vector is: $\boldsymbol{r}(t)=\boldsymbol{r}[x(t),y(t),z(t)]$ I have ...
0
votes
1answer
48 views
4
votes
5answers
252 views

On finding all the positive integral solutions of $x^2+y^2=z^2+1$

I am a new to Mathematica. My goal is to find many (if not all) positive integer solutions to the equation: $x^2+y^2=z^2+1$ using Mathematica. However the problem is that I can only find a ...
3
votes
2answers
186 views

Creepy numerical artefacts of a closed-form solution

I'm wishing you a nice day and here is my problem: I tried to obtain relation between chemical potential and temperature in a semiconductor and I arrived at this equation: $$ 2 e^{-\Delta/kT} \left( ...
0
votes
1answer
69 views

Solving a vector equation [duplicate]

I want a simple and intuitive way of solving for $t$ given that $\mathbf{i},\, \mathbf{j},\, \mathbf{k}$ are the unit vectors in $\mathbf{R}^3$. $$ 5 = (2 \mathbf{i} - 3 \mathbf{k} + (\mathbf{i} - ...
0
votes
0answers
62 views

How can I be sure Solve failed because there was no solution to be found?

Is there a way I can be sure that the reason Mathematica can't find an answer is because there isn't one, rather than because I made a mistake with the code? The system is as follows. ...
3
votes
1answer
86 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 ...
0
votes
0answers
51 views

Solve equations for derivatives

I want to solve some equations in General Relativity numerically with the Runge-Kutta method. Therefore, I need to solve for example $r^{''} \left(\tau\right) + t\left(\tau\right) + a = 3 $, where $ ...
0
votes
0answers
97 views

NSolve and Loop - getting error that does not occur with single instance

I am having difficulties getting this code to work. Code: ...
3
votes
1answer
163 views

Root finding: zeroes of Mathieu function

I am finding the roots of the Mathieu sine function, and find Mathematica and Maple do not agree on the solutions. For example, consider the solutions of ...
0
votes
0answers
57 views

Computing the coefficients for Gaussian integral

I wrote a small script to find the Gaussian integral points and weights: ...
0
votes
0answers
63 views

Plot zero contour of a function of two variables

I have F[x_,y_] = f(x,y). I would like to solve (numerically) F[x,y] == 0 to give x as ...
1
vote
1answer
74 views

Computing a recurrence table backward

I am trying to understand how to compute a recursive equation backward say n = 10 to n = 1. I am trying to solve a PDE for which ...
1
vote
1answer
96 views

Solve system of equations with implicit function [closed]

I am trying to solve a system of 3 (non-linear) equations with two variables $x,\ y$ and one unknown function $f(x,y)$ of these two variables. Also, the system includes the change of this function $f$ ...
1
vote
3answers
90 views

Getting a list of x coordinates of intersection points of two functions

I have two functions, f[x_] := 5 Sin[2 x] g[x_] := (5 x + 2)/(x + 2) I am trying to find the x coordinates that intersect ...
0
votes
1answer
177 views

RSolve with integrations

I want to solve this question with Mathematica. $$T_0(t)=c(t),\quad T_n(t)=a\int_0^t T_{n-1}(s)\, \mathrm{d}s+b(t).$$ Is it possible to solve $T_n(t)$ with ...
0
votes
0answers
45 views

RSolve with two variables

I want to find $a(m,n)$ with $~a(m,n)=\alpha * a(m-1,n) + \beta * a(m,n-1)$, where $\alpha$ and $\beta$ are constants. I expect an explicit $a(m,n)$ represented by $a(0,j),~j=1,2,\cdots,n$ and ...
0
votes
0answers
40 views

Deciphering output messages using Solve

I'm trying to solve the following system. The error messages I'm getting suggest that there may not be a solution [my desired result], but I'm not sure if there could be something else going on. The ...
1
vote
0answers
46 views

How to speed up geometric Resolve query involving $\exists$ and $\forall$?

I would want to test connectedness of semialgebraic sets with naive code like this: ...
0
votes
1answer
55 views

DSolve not outputting any result

I'm trying to solve a system of differential equations using DSolve, but I am not getting any output beyond a repetition of my input: ...
1
vote
0answers
42 views

Tracking solutions over a Brillouin zone

Thanks to @bbgodfrey I have a function that finds the two solutions to a set of equations bounded by a hexagonal region (corresponding to the Brillouin zone). It takes the angle ...
0
votes
1answer
179 views

Add assumptions to Solve command

I want Mathematica to solve the following by assuming the part within the Log function, which is ...
0
votes
0answers
71 views

Difficulty solving an equation involving a couple of different powers

I have this semi-complicated equation to which I want to find an analytical solution for $\epsilon$. The equation is a s following: $$a^2 (2 b k \epsilon \sqrt{1 - \epsilon^2/b^2} - 2 c k ...
0
votes
1answer
83 views