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

learn more… | top users | synonyms (5)

6
votes
2answers
1k views

Solve an equation in $\mathbb{R}^+$

How can I solve this equation (both numerically and literally) only in the positive reals $\mathbb{R}^+$? Solve[x == (v0 - (A CD t v0^2 ρ)/(4m)) Cos[θ] t, t] And ...
6
votes
3answers
2k views

Why won't Mathematica Solve a set of two equations for one variable?

For example, if I try to Solve this set of equations Solve[y == 3 x + 5 && y == -x + 7, {x, y}] Mathematica gives the right values for ...
6
votes
5answers
776 views

Can someone help me solve this x*y*z=180

I have one question I've been trying for a while but I can't get the result for my equation.. So I've been trying hard but I can calculate only by hand but the problem is that I am afraid to miss any ...
6
votes
2answers
490 views

How to work out the parameter in a definite integration which has an exact value while the integration doesn't have an analytical solution?

Here is the equation I'm trying to solve: NIntegrate[1/(E^(1/(λ T)) - 1), {λ, 200, 220}] == 1000 T is the parameter I'm ...
6
votes
3answers
495 views

Plotting $\omega(k_x, k_z)$ in $(\omega, k)$ plane with assumption of $k^2 = k_x^2 + k_z^2$

I'm currently studying atmospheric-gravity waves and their dispersion relations. For analysis I need to plot the dispersion function: \begin{equation} \omega^4 - \omega^2\cdot C_s^2(k_x^2 + ...
6
votes
4answers
454 views

Manipulating Equations

I would like to be able to manipulate algebraic equations. For example, let's say I enter 3x+5==2 I would next like to be able to manually subtract $2$ from both ...
6
votes
3answers
244 views

Labeling solutions of an Eigenvalue equation involving Bessel functions

I'm solving the Schrödinger equation for a particle in an annular geometry with hard wall boundary conditions and I've reduced it to the following equation: $$J_m(k\,R_1)\,Y_m(k\,R_2) - ...
6
votes
4answers
339 views

How to eliminate variables when using Solve[]

In certain problems, we need to solve systems of equations and get results in terms of just selected variables. For example, how could we solve eqn==0 below for ...
6
votes
2answers
2k views

Error entering equation in DSolve

I entered a command incorrectly as follows: DSolve[{y'[x]=y[x]},y[x],x] I am now experiencing: ...
6
votes
2answers
3k views

How can I get Mathematica to solve an equation with multiple variables?

How can I get mathematica to solve the non-homogeneous differential equation with undetermined coefficients listed below? Solve $\ y''(t)+3y'(t) = 2t^4 $, using: $\ y_p(t) = ...
6
votes
3answers
836 views

Arithmetic operation on the value returned by Solve

I want to modify the value that is returned by Solve. Example: SolC1 = Solve[Sin[ π x / L] == 1, x] {{x -> L/2}} I need to add n to x so that x becomes ...
6
votes
2answers
256 views

Solving a two variable two equation trig expression

I am trying to use Solve to get what two angles (A and B) are in a problem I have. ...
6
votes
1answer
185 views

Monitoring the progress of FindRoot during evaluation [closed]

I would like to know whether there is a way to monitor the evaluation of the norm of the error of FindRoot. Sometimes the evaluation of of this for a large system ...
6
votes
2answers
484 views

Reduce and Solve unable to solve equation

I have the following equation. Why aren't Solve and Reduce able to find the solution? ...
6
votes
3answers
393 views

How to get the correct answer when solve this equation?

I want to find the values of the parameter $m$ to the graph of the funtion $$x^4 - 2 m x^2 + m + 1$$ has three extremal points make an equilateral triangle. I tried ...
6
votes
2answers
2k views

How to invert an integral equation

There have been numerous times when I've needed to invert an integral equation, i.e. I have something like $$f(x) = g_1(x)\int_{0}^x g_2(x') dx'$$ for arbitrary functions $g_1$ and $g_2$, and ...
6
votes
1answer
173 views

Why is FindInstance failing when I relax a set of constraints?

I'm attempting to use FindInstance to generate coordinate sets for plausible triangles with edge length distance constraints. E.g.: ...
6
votes
2answers
238 views

How to reverse irreversible function in Mathematica?

How to reverse formula $y(x)=x (\frac{1}{sin \frac{\pi x}{2}})^\alpha$ i.e. express it as $x = x(y)$ in Mathematica? I did this way ...
6
votes
1answer
1k views

How Can I use Solve/Reduce Output

Suppose I want x and y to be rationals Solve[ x^2 + y^2 == 1, {x, y}, Rationals] I am ...
6
votes
2answers
551 views

Very long Refine/Solve batch run - is my code broken, or just complicated?

So I'm trying to run some Mathematica code in batch mode from my university cluster. Specifically, I'm trying to find the equilibria of a system of ordinary differential equations. Inspired by the ...
6
votes
2answers
348 views

NSolve missing solutions in Mathematica 10

Running Mathematica 8.0.4 and 10.0.0 on a Windows 8.1 machine. Processed the same code with both kernels: ...
6
votes
2answers
411 views

Why does Equal[“string”, False] not evaluate to False?

Why does Equal["string", False] not evaluate to False? I would have thought Mathematica could prove that these two things are ...
6
votes
3answers
71 views

Using Select on a 2d table with a criteria that accepts 3 parameters

I have a table sets := Table[{i, j, k}, {i, 1, 10}, {j, 1, 10}, {k, 1, 10}] and I want to use a function that accepts three parameters and returns a boolean, as ...
6
votes
1answer
734 views

Issue with NSolve

I defined the following function: sw[ θ_, ϕ_, h_] = 1/2 Sin[2 (ϕ - θ)] + h Sin[ϕ] Based upon whether I use Solve or NSolve, ...
6
votes
1answer
881 views

FindRoot options

I found an example of FindRoot usage that I couldn't find documented. It looks like: FindRoot[f[x] == 0.4, {x, 0, 1, 0, 1}] ...
6
votes
1answer
130 views

Why is Solve slow for one representation and fast for an equivalent form?

When I try to solve the following: Solve[(1 - (-1 + x)^- x)^8 == 95/100 && x > 1, Reals] solve chokes. If instead I use the equivalent form: ...
6
votes
1answer
243 views

Solving $\frac{dx}{dt} = A \frac{ (1-x)}{(t-t^2)} - \frac{(B*x -C*x^2)}{(t-t^2 )*(t-x)}$

I would like to solve the following equation: $\frac{dx}{dt} = A \frac{ (1-x)}{(t-t^2)} - \frac{(B*x -C*x^2)}{(t-t^2 )*(t-x)}$ I have already posted it here, but still it doesn't work for me (I ...
6
votes
4answers
791 views

How to find an approximate solution in a region when Solve or NSolve does not work?

Consider the relation x/Tan[Pi x] == y/Tan[Pi y] Solve or NSolve can't give a solution ...
6
votes
1answer
181 views

The third argument of the Root function

Consider the roots of an arbitrary indecomposable polynomial: sol = Solve[x^5 + 2 x + 1 == 0, x] The returned expression is ...
6
votes
2answers
1k views

Solution of equation with power series (perturbation)

So I need to use Mathematica to find the solution of $y=x- \epsilon \sin(2y)$ as a power series in terms of $\epsilon$. I'd assume I'd need to create an equation $f=x-y- \epsilon \sin(2y)$, then ...
6
votes
1answer
2k views

Efficient implementation of cubic equation of state

...In this case the Peng Robinson Equation of State. Equations of state are empirical equations with parameters derived from experimental data. They are used to predict pure component and mixture ...
6
votes
1answer
120 views

Closed form of two-layer recurrence

I am trying to find a closed-form expression for $A(0)$, described by the following recurrence forms, as a function of $L$ and $d$: \begin{align} A(\ell) &= 1 + 2^d A(\ell+1) + \sum_{g=0}^{d-1} ...
6
votes
2answers
855 views

Find all the integer numbers $a$, $b$, $c$, $d$, $e$, $f$, $k$ to this equation have three integer different solutions?

How to find all integer numbers $a$, $b$, $c$, $d$, $e$, $f$, $k$ that belong to interval $[-10,10]$, so that the equation $$\sqrt[3]{a x + b} +\sqrt[3]{c x + d} + \sqrt[3]{e x + f} =k$$ has three ...
6
votes
2answers
339 views

Time evolution/dynamics of circular plate with hole (Biharmonic equation and stiffness)

On this Mathematica.SE form, there exists information on how to use Mathematica to demonstrate the vibration of a circular membrane and the deflection of an orifice plate (the latter I had raised ...
6
votes
2answers
153 views

Why is ToRadicals[] not able to handle all cases? Is there a workaround?

The documentation for the function ToRadicals says: ♦ There are some cases in which expressions involving radicals can in principle be given, but ...
6
votes
2answers
170 views

NSolve erroneously gives no solution to a polynomial system

I have a polynomial system with three equations in three unknowns, the maximum degree is 26. Two equations are symmetric, i.e. eq1(x,y,z)=eq2(y,x,z). If I search ...
6
votes
4answers
601 views

FindInstance (and Solve, …) abysmally slow on a fully determined system of linear equations and inequalities: why?

I took a quick look for FindInstance and Solve-related questions but came up empty for the aspect that I am curious about. So ...
6
votes
1answer
383 views

Buggy Reduce/Resolve function?

It seems to me that this is a clear-cut case of Mathematica actually producing the wrong answer (with no warnings). I'm trying to express the fact that for any integer, there exists a factorization ...
6
votes
1answer
454 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 ...
6
votes
0answers
127 views

Heat convection differential equations from 1952 - Mathematica “fails to converge”

I am trying to solve a fundamental problem in analytical convective heat transfer: laminar free convection flow and heat transfer from a flat plate parallel to the direction of the generating body ...
6
votes
0answers
207 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( ...
5
votes
4answers
2k views

How to find the sum all even numbers of this sequence?

I have a sequence $(u_{n})$ $$u_1= 1, \quad u_2 = 2, \quad u_3 = 3, \quad u_{n}= -u_{n-3} + 3u _{n-2} +2 u_{n-1}, \quad \forall n \geqslant 4.$$ I want to list the first $20$ terms of this sequence ...
5
votes
4answers
237 views

Finding $y$ coordinates for $\sin(x+y)=x^2+y^2$

How can I find the $y$ coordinates of the equation $\sin(x+y)=x^2+y^2$ when $x=0.5$? I've tried NSolve[Sin[x + y] = x^2 + y^2 /. x -> 0.5, y]
5
votes
3answers
1k views

Solve for the intersection point of 3 surfaces

Basically, I have three 2-torus (tilted) in 3-dimensional Euclidean space. They are expressed by parametric equations: ...
5
votes
2answers
571 views

Order of a pole

Is there a simple way to determine the order of the poles of a rational function? I have a difficult function where I need Mathematica to find the poles. It would be interesting to also know what the ...
5
votes
2answers
318 views

DSolve returns a “wrong” solution

This wrong solution exists from mma9 to mma10.1: ...
5
votes
2answers
425 views

Equation solving with GroebnerBasis

I have a system of quadratic equations defined below ...
5
votes
2answers
180 views

Reduce applied to vectors doesn't work right?

Reduce[{x,0}=={0,0}, x] gives x==0 as it should, but Reduce[{x,0}!={0,0}, x] gives False, which indicates no solution ...
5
votes
2answers
296 views

How to solve this combination equation?

I want to solve the equation $$ 2n\cdot\binom{n}{0} + 5(n-1)\binom{n}{1}+13(n-2)\binom{n}{2}+\cdots + \left(2^{n-1} + 3^{n-1}\right)\binom{n}{n-1} = 1685. $$ or using alternative notion for binomial ...
5
votes
2answers
273 views

When does the real part of Zeta vanish on the critical line?

This seems to be a quite a simple problem but I cannot make it work. I am trying to find all values within a given range for which the real part of the Zeta ...