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
1answer
261 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
376 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
988 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
317 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
158 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
80 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
139 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
582 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
4answers
436 views

Am I missing anything? Solving Equations

I want to know what $ a^4+b^4+c^4 $ would be when $ a+b+c=1, a^2+b^2+c^2=2, a^3+b^3+c^3=3. $ I tried ...
6
votes
3answers
444 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
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
837 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
498 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
514 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
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) = ...
6
votes
3answers
862 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
312 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
217 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
490 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
412 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
174 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
251 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
567 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
152 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
250 views
6
votes
2answers
373 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
440 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
72 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
912 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
150 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
246 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
845 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
184 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
123 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
501 views

Adding noise to nonlinear control system using NDSolve

I have a nonlinear control system ...
6
votes
2answers
886 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
368 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
157 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
209 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
638 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
140 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 ...
6
votes
1answer
390 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
513 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
141 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 ...