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

learn more… | top users | synonyms (5)

2
votes
0answers
1k 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 ...
4
votes
3answers
502 views

Finding parameters making real part of eigenvalues vanish

I have the following $\;3\times3$ matrix: $\left( \begin{array}{ccc} 0.04 -0.4 b & 0 & 0.04 -0.4 b \\ 0 & -0.08-1.2 b & -0.06-0.9 b \\ 1.04 -0.4 b & 2.08 -0.8 b & 0 ...
0
votes
1answer
66 views

Help with Reduce, quadratic system of inequalities

I'm trying to obtain an equivalent system of inequalities to the one below but without the variables x and z. Somehow Reduce just repeats the command. ...
0
votes
4answers
161 views

RSolve fails on a two-equation non-linear system

I have a system of just two recurrence equations, I need to solve them and so I use RSolve: ...
1
vote
1answer
536 views

Why my differential equations become True? [duplicate]

I've been trying to solve a system of nonlinear differential equation, but the conditions are a bit weird. Two of the differentials equate to the same equation, but have different boundary ...
1
vote
1answer
312 views

Get Cos[2 x] in “Solve”

My goal is to get an analytical solution of x ( $\lambda$ in the problem), when $\epsilon$ (f[x] in the problem) is maximum (use inverse form). I found that if I can use trigonometric identity ...
0
votes
1answer
109 views

Equation in implicit form

This is simple example: v = {1.48622576, 1.53444276, 1.42299438} v/(v + 1) Out[1]:= {0.597784, 0.605436, 0.587288} In the above case expression v/(v + 1) is ...
0
votes
1answer
666 views

Solving a system of equations where two variables can only be one of several discrete values

I want to solve a three-equation system, but I need to constrain two of the variables so that they can only take one of several discrete values that I would specify (this is a circuit design problem ...
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 ...
-1
votes
1answer
232 views

Mathematica Compute Time Is Extremely Long (Integrate simple function with parameter scan)

I was wondering if you can help. I have a set of equations I need to solve for a list of parameters; entering each parameter value and solving anew each time is simply not economical. I found a way ...
2
votes
0answers
53 views

How can I select two lists from a list of lists so that I receive the system of equations has only rational solutions?

I want to make a sytem of equations has the form $$\begin{cases} a_1x^2+b_1y^2 + c_1xy+d_1x + e_1y+3=0,\\ a_2x^2+b_2y^2 + c_2xy+d_2x + e_2y+3=0. \end{cases}$$ so that the given system of equations has ...
1
vote
1answer
77 views

Solve never terminates with simultaneous equation set

I have the following set of simultaneous equations that I want to solve symbolically. The two unknowns are $a$ and $d$. $a=\frac{c*(d/10)}{1 + (f/130) + (d/50)} - \frac{a}{1 + (a/100) + (d/100) + ...
1
vote
1answer
425 views

Marking the points of x-intercept on a plot

OK, I have a nice plot of some oscillating motion. I can read some stuff off the graph, but I thought it would be cool to show a dot where the curve hits the x-axis and a little label showing the ...
1
vote
1answer
181 views

How do assign an interval of values to a variable in order to use them in a function?

From my code (below) I have K1 and K2 set to zero in order to test this numerical method. However, K1 is on an interval from 0 to 32, and K2 is on an interval from 0 to 8. I would like to know if ...
1
vote
1answer
114 views

A simple expression which Mathematica will not simplify

Sorry for the vague title, but I wasn't sure how to phrase this. I was using Mathematica to solve a fairly simple algebraic equation, but when I tried to substitute the solution into the original ...
5
votes
1answer
503 views

Solving a nonlinear, second order d.e. in Mathematica

I am trying to solve the following differential equation $$ s''(x) = \frac{8}{9}s(x) -2s^2(x) + s^3(x), $$ with boundary conditions $s(-30) = 4/3 $ and $s(30)=0$, in Mathematica by using the following ...
0
votes
3answers
606 views

How to solve correctly inside a DO loop and then plot

Here is the corresponding code First the equations setup ...
1
vote
0answers
23 views

Solving a pair of recurrence relations [duplicate]

Why does this RSolve[{y[t] == a*y[t - 1]}, {y[t]}, t] work, while this ...
11
votes
3answers
688 views

1D quasicrystal: points on a line nearest to points on a lattice

I have a simple lattice / line manipulation: ...
0
votes
1answer
275 views

Solving a system of two equations of third degree

I am trying to solve the following system of two equations for x and y. The solutions for x and y have to be positive and the parameters s and c are positive (to have a coherent interpretation of my ...
1
vote
1answer
231 views

Solution of Coupled Second Order ODEs with Boundary Conditions

D[Tu1[x], {x, 2}] qu - Tu1[x] (fu0 + I w p) == Tou[x] f1 D[T1[x], {x, 2}] q - T1[x] (f0 - b g0 + I w p)== To[x] (f1 - b g1) - g1 For T1 , ...
0
votes
0answers
340 views

Problem with Plot and ParametricNDSolve

This is my first question on this site, I hope it is written in an understandable way. Let's start. My aim is to plot a certain parametric region, with $x$ and $y$ coordinates given by the following ...
0
votes
1answer
103 views

Coefficient of function different than polynomial

My question arise from the problem that Mathematica can't solve the equation $$ \text{Solve}\left[e^{a-x} (0.0299584 a+0.1 b-0.0299584 x-0.152913)-0.9 a+0.9 b+0.9 x+0.18=0.0498614 e^{-x} \left(e^a ...
2
votes
1answer
649 views

How do I get all possible solutions in an underdetermined system? [duplicate]

I have two problems which I'd like to solve with Mathematica. If I have a system of two equations with three unknowns, how can I get to list all possible solutions for the unknowns? Here is what I ...
3
votes
1answer
1k views

Solving a differential equation with initial conditions

I want to get a solution to an ordinary differential equation with imposed initial conditions. I'm confused about why my following code doesn't work: ...
1
vote
1answer
390 views

Solution of an ODE in implicit form

For some non-linear ODEs there is only implicit form of solution using DSolve. For example DSolve[(y[x] + x - 1)*y'[x] - y[x] + 2 x + 3 == 0, y[x], x] gives ...
1
vote
2answers
595 views

Given an equation with one parameter, how do I use NSolve to solve it for many values of the parameter?

I was wondering if there was a way to use the Mathematica NSolve command to solve an equation that has a parameter with various values. An example of which is: ...
0
votes
1answer
44 views

Storing solution to x and y [duplicate]

If I do sol=Solve[{2x+2==y,3x+9==y},{x,y}] then Mathematica will tell me {{x -> -7, y -> -12}} But I would like to ...
1
vote
0answers
99 views

Finding undetermined coefficients with boundary conditions

I have 2 equations which are; To[x_] := c1 E^(-a1 x)/2 + c2 E^(a1 x)/2 + a2 Tu[x_] := c3 E^(-a3 x)/2 + c4 E^(a3 x)/2 Tu equals almost zero ...
2
votes
1answer
225 views

How to solve this numerically?

May I ask what is the best way to evaluate ...
2
votes
2answers
171 views

Solving for any one variable in a formula given values for all the others

Suppose I have an equation relating a number of variables and I know the values of all but one of them. Using the equation and known values as input, is there a simple command (built-in or otherwise) ...
0
votes
0answers
51 views
0
votes
1answer
253 views

Minimize: How to work with specific numeric types

This is a question about the use of symbols in Minimize and how to persuade Mathematica to treat them in a specific way (as integers or rationals) that I have not ...
2
votes
1answer
148 views

Equation involving hypergeometric functions

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

How to obtain all minima? [duplicate]

If you Solve an equation, you obtain all solutions. But If you NMinimize a function, you only get one minimum. The question ...
1
vote
2answers
299 views

Mathematica can't minimize a function

Mathematica seems not to be able to minimize this univariate function over integer arguments, $r>2, r \in \mathbb{Z}$. ...
5
votes
2answers
524 views

Finding Points of Tangency of a Circle and a Parametric Curve

I have the following Parametric plot Generated from the following: ...
0
votes
1answer
319 views

Numerical Methods - minimum of a function

Problem: The length of the longest ladder that can negotiate the corner depicted in Figure 2 can be determined by computing the value of theta that minimizes the following function: ...
1
vote
1answer
781 views

Trying to solve a transcendental equation involving Bessel functions

I've never used Mathematica before and am trying to numerically solve equation (12) from this paper. Ideally I'd be able to find the smallest value of $x_{n\nu}$ for $\exp(-kr\pi)$ close to 1, and ...
9
votes
1answer
250 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 ...
4
votes
3answers
319 views

FindRoot constrained to maximum

I am using FindRoot to find the local maximum op some function, by examining where its derivative is zero. However, keeping the starting value the same, sometimes ...
2
votes
1answer
522 views

Why does Mathematice dislike my boundary condition in NDSolve?

im trying to solve a PDE-System, with the "periodic boundary condition" c[t,0]==4*c[t,2pi] But mathenmatica says: NDSolve::bcedge: "Boundary condition ...
2
votes
2answers
889 views

Solving simultaneous exponential equations

In[1]: Solve[{(1 + k)^n == 0.6, (1 + k)^(n + 1) == 0.3}, {n, k}] Solve::ifun: Inverse functions are being used by Solve, so some solutions may not be found; ...
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 ...
0
votes
1answer
1k views

How to find the inverse a symbolic matrix that Mathematica thinks might be singular

I have a square matrix depending on a parameter "a". When "a" has a numeric value (eg 2.1 * 10 ^ 8) through the inverse function Inverse[] I can invert the matrix. If I write the matrix as a function ...
1
vote
2answers
263 views

Symbolic solution of irrational equation

Is there any a nice looking solution to the following equation? eq[x_] := 2 x q1 - 2 q2/x^3 - (4 x^3 q3 - 4 q5/x^5)/(4 Sqrt[x^4 q3 + q4 + q5/x^4]) Better than ...
1
vote
2answers
990 views

Why do Solve, NSolve, and FindRoot all fail on this simple, solvable actuarial equation?

If we express conditional mortality as a vector of annual probabilities of death, like so ...
1
vote
1answer
118 views

Having trouble assigning the result of a Solve to a function

Can anyone see what is wrong with my syntax? It looks like I have successfully assigned the result of Solve to a function. I want Mathematica to return me an numeric value of the function ...
1
vote
1answer
84 views

Invariant error plot for an arbitray system of ODE

i try to use InvariantErrorPlot command for a generic system of ODE but i find only example with pre-building equations from Mathematica Packages. This is one of ...
0
votes
0answers
215 views

Mapping the imaginary axis of the complex (s) plane

In this problem is very difficult to analyze the roots s. Firstly I thought that the problem is in precision, but I am not sure. The criterion which I want to obtain is that one of these roots has a ...