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
96 views

Solve or preferable Reduce equation

I am new to mathematica and trying to solve and reduce the following: Solve[1 - (1 - p^y)^(u (1 - p)) == a, y] or ...
4
votes
2answers
303 views

Why doesn't Roots work on a certain quartic polynomial equation?

In everything that follows I am using Roots to get the solutions to the equations. I start off with this equation: $$ \frac{A}{3}x^4 - x^2 + 2ax - b^2 = 0$$ Now ...
4
votes
2answers
157 views

Is there a giant AND symbol/How to solve tons of equations simultaneously

$$\sum_{k=1}^na_k=a_1+...+a_n$$ $$\prod_{k=1}^na_k=a_1\cdot...\cdot a_n$$ $$\large{?}\small{_{k=1}^na_k=a_1~\&\&...\&\& ~a_n}$$ What in mathematica would allow me to index logical ...
4
votes
1answer
229 views

Numerical errors/inaccuracies in ProductLog

Context In cosmology, a fairly accurate model to describe the gravitational potential, $\psi(r)$ of dark matter halos is given by $\psi( r)=\log(1+r)/r$. ...
4
votes
2answers
265 views

How can I eliminate equivalent equations from a list?

I want to find the numbers $a$, $b$, $c$, $d$ of the plane $ax + by +cz +d = 0$ which make the distance from the point $M(1,2,3)$ to the plane equal to 2. I tried ...
4
votes
2answers
169 views

Specifying and solving a PDE on an L-shaped domain

I wish to solve Laplace's equation with Dirchlet boundary conditions on an L-shaped domain as described here http://www.inf.ethz.ch/personal/tulink/FEM14/Ch1_ElmanSyvesterWathen_Ox05.pdf I know how ...
4
votes
2answers
142 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 ...
4
votes
2answers
216 views

How to set up and minimize a large system of equations to resolve the individual components

I am trying to set up a system of 400*N equations to deconvolute spectral data for four different compounds. I monitored 400 different wavelengths in N solvents (3 solvents in this case). In order ...
4
votes
2answers
122 views

Non-convergance of the solution of the 1+1 partial differential equation

I am stacked with a partial differential equation as follows: This is the $\LaTeX$ expression for the equation $$ \frac{\partial y}{\partial t}=\frac{\partial^2 y(x,t)}{\partial x^2} ...
4
votes
1answer
197 views

Check if an inequality holds among the solutions returned by Solve [duplicate]

Let's say I solve a system: Solve[{a == 3* c, b == 2 *a}, {a, b}] and then want to see if the values found for a and ...
4
votes
1answer
2k views

Finding the Minimum value of an interpolating function

I can't seem to use FindMinValue to find the min. value of a curve represented by an interpolating function. For instance the below code generates an interpolating ...
4
votes
1answer
79 views

Trouble using FindInstance with MatrixRank

I'm really confused about what's going wrong with my call to FindInstance: ...
4
votes
3answers
1k views

Backsubstituting solution into FindRoot

I am trying to solve a system of equations which is dependent to an parameter alpha. So I iterate over alpha via a ...
4
votes
2answers
139 views

Can't solve equation having complex coefficients

Why can't Mathematica solve $\quad\quad\frac{1-i}{\sqrt{2}}=e^{i \alpha } \tan \left(\frac{\beta }{2}\right)$ with the restrictions $\alpha \in [0, 2 \, \pi)$ and $\beta \in [0,\pi]$: ...
4
votes
1answer
156 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
1answer
1k views

What do the hashtags represent in the result of Solve, and how should I use them? [duplicate]

There is a function that I'm trying to get the inverse of. My function is a simple polynomial: Tn[T_] := 4/9 T^3/Tc^2 + 5/9 T Now I tried to find ...
4
votes
1answer
608 views

How do I use RSolve to solve a system of recurrence relations?

I am trying to solve a system of recurrence relations as follows. ...
4
votes
2answers
315 views

Both Sin[x]==0 && Cos[x]==0 as a solution

If I do FullSimplify[Reduce[Sin[p1] == 0 && Cos[p1] == 0, Reals]] I get Cos[p1] == 0 && Sin[p1] == 0 ...
4
votes
1answer
354 views

pH calculations, speeding up code involving Solve

I'm trying to calculate pH values of solutions containing water and acetic acid. The relevant equations have been defined and I now want to find a solution using Mathematica. For a specific acid ...
4
votes
1answer
135 views

How to solve a conditional inequality?

I am trying to use Mathematica to solve a conditional inequality. Find all (real) numbers $a$ and $b$ such that $|a| + |b| \ge 2/ \sqrt{3}$ and for any real $x$ the inequality $ |a\sin x + b \sin ...
4
votes
1answer
197 views

Asking Mathematica to produce a function of some variables

Given $$l_0=\sqrt{\frac{\lambda_1(\lambda_0+\mu_0)}{\lambda_0(\lambda_1+\mu_1-1)}}$$ and $$l_1=\sqrt{\frac{\lambda_0(\lambda_1+\mu_1)}{\lambda_1(\lambda_0+\mu_0-1)}}$$ ...
4
votes
1answer
220 views

How to demonstrate lack of stability with advection equation

I am trying to that using a coarse grid with an explicit method for, say, the advection equation leads to an unstable solution. The trouble is Mathematica avoids unstable solutions with good ...
4
votes
1answer
1k views

Use of $Assumptions [duplicate]

I don't really understand the behaviour and the usage of $Assumption. For instance, when I set $Assumptions = λ > 0, I expect ...
4
votes
1answer
239 views

Can mathematica resolve this equation?

Obviously,the equation $x^2-3 y^2=2 z^2$ doesn't have positive integer solutions. I tried ...
4
votes
1answer
358 views

Eliminating functions from system of PDE in Mathematica

Mathematica users, Let's see an example of three PDEs: How can I eliminate variables q1[x,y,z] and q2[x,y,z] from a system of ...
4
votes
1answer
1k 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 ...
4
votes
1answer
55 views

Sequential System of Equations

I’m trying to set up this system of equations that sequential feeds into itself. I have these 2 equations I’m running. I’m new to Mathematica so I have no idea how to even start this. The basic idea ...
4
votes
1answer
367 views

3D Plot of The Riemann Zeta Function

I am a new user and i would appreciate your help to visualize the 3D of the following function. Edit I tried this ...
4
votes
0answers
93 views

Synchrotron Radiation and ListDensityPlot

I'm trying to reproduce this plots: For this task i'm following this blog entry: Synchrotron radiation As you can see there, the task is to make a DensityPlot of a vector field essentially. This ...
4
votes
0answers
249 views

What is the algorithm behind Mathematicas Reduce in this equation?

Mathematica can use Reduce to arrive at an answer in this equation ...
4
votes
0answers
310 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 ...
4
votes
0answers
168 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( ...
4
votes
0answers
120 views

Symbolic solution to ODE, pure InverseFunction not evaluated

I'm new to Mathematica and I don't understand why in the solution of the following ODE, the #1 in the pure function is not immediately replaced by the corresponding ...
3
votes
2answers
284 views
3
votes
3answers
1k views

Using DSolve to solve for x[t,parameter]

I'm trying to solve an ODE with two independent variables (a cannon firing from a cliff incorporating wind resistance dependent on velocity). I've tried the following for the x-component: ...
3
votes
3answers
286 views

Finding intersections of two parametric curves

I am dealing with the root finding of two parametric data curves (data below). I have already checked Parametric Interpolation of 2D data, which I found useful for plotting, though it is not adequate ...
3
votes
3answers
143 views

Find $Y \in (0, 5)$ such that the roots of $4x^2 + 4xY + Y + 2 = 0$ are real?

The problem is to find all $Y \in (0, 5)$ such that the roots of $4x^2 + 4xY + Y + 2 = 0$ are real. I took a straightforward approach with Mathematica with the intent of listing the only solutions ...
3
votes
2answers
312 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. ...
3
votes
3answers
797 views

How can I solve Tan[t] - t == F[x] for t as a function of x?

How can I solve the equation Tan[t] - t = Ax, where A is a constant for t[x]? I know that ...
3
votes
2answers
3k views

How to convert a system of parametric equations to a normal equation?

For example, I have a system of parametric equations (R is a constant number) : ...
3
votes
3answers
248 views

How to solve the equation $A \sin x + B \cos x = 0$ for all $x$?

I have the equation: $A \sin x + B \cos[x] = 0$ that should be satisfied for all the $x$. Of course, the solution is obvious: $A=B=0$, but I need Mathematica to resolve it automatically because this ...
3
votes
1answer
510 views

How do I solve $1 - (1 - (A x)^2)^\frac{3}{2} - B(1 - \cos(x))= 0$?

Consider the equation $$ 1 - (1 - (A x)^2)^\frac{3}{2} - B(1 - \cos(x)) = 0 $$ where $A,B \in \mathbb{R}$ are constants. What is the analytic solution?
3
votes
2answers
252 views

Mathematica does not understand (R^3)^(1/3) is the same as R [closed]

In the output from a calculation in mathematica stands a/((R^3*c)^(1/3)), with c and a ...
3
votes
1answer
252 views

What properties make this equation difficult to symbolically solve?

I tried to solve the following equation with Mathematica's solve: Solve[K*(2*Tan[L/2*Sqrt[P/(EI)]]-L*Sqrt[P/(EI)])+4*P*Sqrt[P/(EI)] == 0, P] It gave the ...
3
votes
1answer
373 views

Simulating a quick return mechanism

Last time, I got a satisfying answer for simulating a toggle mechanism, just as expected. NDSolve is great. I tried playing with it (as a matter of fact, @MichaelE2 ...
3
votes
1answer
2k views

Finding intersection of two graphs

Here is a toy problem: I have a (directional) line that starts from {2,0} and goes through {1.8, 0.01}. Also I have a ...
3
votes
2answers
259 views

Computing the minimum distance in a contour plot

I have the following Mathematica code ...
3
votes
2answers
724 views

All possible solutions to the Matrix Equation (free variables appearing)

I am attempting to solve a system of linear equations using LinearSolve[] . In my case, the number of unknowns are more than constraints. I learnt that LinearSolve ...
3
votes
1answer
194 views

Very slow mathematica finite differences

I have a simple code which solves an equation by an explicit method (FTCS). It takes mathematica several minutes (mathematica 10.0.2) to finish the calculation while the same code in Fortran runs less ...