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

learn more… | top users | synonyms (5)

7
votes
2answers
221 views

Why does Solve and Sum give different answers when according to Table it is the same sequence?

I have been thinking about a piece of code today, but I don't understand why I get different outputs from two in my mind similar Mathematical sums. First I observe that the Table command for this: <...
7
votes
1answer
327 views

Why does Mathematica not give simple solutions to simple system of equations with powers

I trt to get Mathematica getting me the simple solutions for x and y from this system of two power equations: $$\text{Eq. 1: } z=x^{m}y^{n}$$ $$\text{Eq. 2: }z=\tau x^{m-1}y^{1-n},$$ The simple ...
7
votes
3answers
576 views

NDSolve and {C, K, Slot} and other built-ins as a variable name

The following problem is an exploration of what causes the error "Input is not an ordinary differential equation" in Mathematica as it seems to have changed from version 8 to version 9. Specifically ...
7
votes
1answer
840 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, ...
7
votes
1answer
939 views

How to Simplify equations over a Ring with Mathematica?

For example, when we work over a ring, the equation x^3=0 does not imply x^2=0 or x=0, but ...
7
votes
1answer
251 views

How to plot inequalities taking into account that points must be integers and not reals?

I have solved with Mathematica 9 the following integer inequalities using Reduce: ...
7
votes
3answers
210 views

Ovals of plane curves (and esoteric surfaces in space)

Consider the curve $y^2 - x^3 + x=0$ with $x, y$ real. This curve has two connected components (one bounded, one not), and the question is: how do you get Mathematica to tell you that (if you do: <...
7
votes
1answer
284 views

Synchrotron Radiation and ListDensityPlot

I'm trying to reproduce the following plots: As described by this blog entry: Synchrotron radiation As you can see in the linked article, the task is essentially to make a ...
7
votes
1answer
391 views
7
votes
1answer
2k views

Solve Lagrange multipliers optimization problem

I have two nested solid figure, where $V(a,h,\tau)$ defines the volume and $A(a,h,t)$ defines the surface. The outer solid figure is parametrized in $a_s$,$h_s$ and $t_s$ (they share a common center). ...
7
votes
1answer
1k views

Using Mathematica to help to determine the consistency of and numerically solve systems of non-linear equations

The formal problem The question I am about to ask is applicable to many scientific and engineering problems. In short, I am bout to ask for pointers and practical advice for Mathematica beginners who ...
7
votes
1answer
336 views

How to find the positive region of this function?

I'm new-ish to Mathematica and I'm using it to learn the basics of derivative trading (for fun not for trading!). I've defined a simple function to calculate the pay off of a trading on a call. I have ...
7
votes
1answer
159 views

Recurrences of lists

Assume we have a multi-dimensional recurrence, e.g. $\qquad\begin{align*} a_1 &= (1,2) \\ a_n &= (1,2) + a_{n-1} \quad, n>1 \end{align*}$ with the easy solution $a_n = (n,2n)$. How ...
7
votes
0answers
82 views

Choosing appropriate WorkingPrecision when solving numerical system of equations [closed]

Consider a linear system of equations, which is conveniently written as A.x=y. The matrix A has dimensions ...
7
votes
0answers
151 views

How does Mathematica resolve symbolic systems of inequalities? [closed]

Before getting into my question, let me provide some context: I recently started using Mathematica to automatically perform stability analyses that I previously did by hand. A major part of the ...
7
votes
0answers
591 views

Creating and using an explicit piecewise function in a convenient way [closed]

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

Why don't I get numeric values from this Solve?

Solve[ a^2 + b^2 == c^2 && a < 100 && b < 100 && c < 100, {a, b, c}, Integers] gives me this result ...
6
votes
3answers
212 views

Find the number of $n$ such that $n!$ is a sum of three squares

I want to check the following theorem by using Mathematica: $\textbf{Theorem} $. $\text{The estimate}$ $\# \{n \le x:n! \text{ is a sum of three squares}\}=7x/8+O(x^{2/3})$ $\text{holds.}...
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
860 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
500 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
531 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 + k_z^...
6
votes
2answers
275 views

Find a condition that b must satisfy so that Ax=b has solution

I'm new to Mathematica, so I'm sorry if this is really simple. I am trying to find the condition that vector b must satisfy so that Ax=b has solution. I would like to learn a general method, but I'll ...
6
votes
2answers
4k 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) = t(A_0t^4+A_1t^3+A_2t^2+...
6
votes
3answers
896 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
350 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
262 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
496 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
427 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
1answer
597 views

Using external variable as a variable inside WhenEvent & NDSolve

I'm trying to use the new WhenEvent functionality of NDSolve in Mathematica 9, in order to perform an action when a variable in the system reaches a certain value. I would like this variable to come ...
6
votes
2answers
3k 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
2answers
170 views

NSolve behaving differently in V9, V10.3 and V10.4

EDIT In a very fast response from Wolfram they state that this is a known problem and give an example of it working for version 10.4.1. $Version "10.4.1 for Microsoft Windows (64-bit) (April ...
6
votes
1answer
175 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
268 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
2k 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
583 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
159 views

NDSolve and matrix formulation do not produce the same result. Why and which is more precise?

Let's assume a vector of variables $\vec a=(a_1(t),a_2(t),a_3(t))$ and let's assume a following DE$$C\ddot{\vec a}=B\dot{\vec a}+A\vec{a}+\vec F.$$ If $A$ and $C$ are both reversible, than the DE can ...
6
votes
1answer
255 views
6
votes
2answers
394 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
461 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
1answer
950 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
3answers
208 views

Approximate analytic solution to a polynomial equation

The analytic solution of the equation $$ 1-\frac{2M}{r}+\frac{Q^2}{r^2}-\frac{\Lambda}{3}r^2 = 0 $$ obtained by using the code ...
6
votes
3answers
75 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
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
166 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
250 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 don'...
6
votes
4answers
917 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
185 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
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 ...