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
2answers
404 views

How to explain NDSolve that functions must be positive only?

I have a simple system of differential equations to solve but both functions should be positive: x[t] > 0 and y[t] > 0 for all t. I was not able to explain NDSolve this fact so it goes x[t] < 0 ...
2
votes
1answer
279 views

Find solutions of equation involving Bessel functions

I'm new in Mathematica and I'm trying to find the solutions of this equation involving Bessel functions $$\eta \frac{ J_{n+1}(\eta a)}{J_n(\eta a)}+\chi \frac{ I_{n+1}(\chi a)}{I_n(\chi a)}=0$$ ...
2
votes
1answer
168 views

Simple differential equation to solve

Suppose we have this differential equation: $$ x^2 + y^2 + z^2 = \frac{\frac{\partial y}{\partial x}}{x+y+z}$$ I want to find $\frac{\partial y}{\partial x}$ at x=1,y=1,z=1. I tried this code but ...
2
votes
2answers
177 views

Graphical solution of a a two-variable system containing an equation and an inequality

I'm currently studying atmospheric-gravity waves and their dispersion relations. I want to have look at values of $(x,y)$ that solve the following system : \begin{equation}\label{}(1) \left\{ \...
2
votes
2answers
194 views

Quadratic Equations $\bmod p^k$

As part of a larger program, I need to solve $A x^2+B x+C \equiv 0 \pmod {p^k}$ for prime $p$. Right now I'm doing this by calling ...
2
votes
1answer
579 views

Resonance Curve using NDSolve

I am trying to make a resonance curve like this one below, using NDSolve. I tried to use code as below. Basically I try to hold solution for a w (driving frequency value) using NumericQ and then ...
2
votes
1answer
117 views

Plotting the solution of a 2nd order equation

I have got to solve a second order equation, where the coefficient depend on the parameter $k$: $$(A_k - \omega) (D_k- \omega) - B_kC_k = 0$$ so that also the solution $\omega=\omega(k)$ will depend ...
2
votes
1answer
488 views

Solving System of Inequalities Relatively

I have been trying to figure out how to solve a system of inequalities relatively. I don't care about values, I just want all possible ordering. An example system might be:$${b > r,\; p > r,\; c ...
2
votes
1answer
146 views

Compact form for domain ranges of many variables

When one sets up an equation to solve one often needs to specify conditions on the variables such as { x>0 && y>0 && z>0 } This is not ...
2
votes
2answers
798 views

How to make a Dynamic Calculator to solve a quadratic equation?

Taking my first steps in the world of Mathematica Programming. The aim of my program was to modify Wolfram's dynamic calculator (which performs an addition subtraction and a product - in response ...
2
votes
1answer
284 views

FindRoot equation-variable mismatch

I cannot figure out why FindRoot doesn't work and returns this error: The number of equations does not match the number of variables in ... My problem: drawing ...
2
votes
1answer
479 views

DSolve expression doesn't evaluate

Can Mathematica solve $\; d^2/dx^2{}+ (1/x) d/dx\;$ form of equations? For example, I am trying to solve the system of equations given below: ...
2
votes
2answers
2k views

How could I define a function as the solution of equations in Mathematica?

I have posed a specific question yesterday but it may be too trivial to answer. I think the question could be asked in a more general way so it is easier to answer. Then I could solve the original ...
2
votes
1answer
240 views

Finding the root of a nested function with small values

I have reduced an error in my program to this line of code: FindRoot[Nest[# (1 - #) k &, 1/2, 2^4] - 1/2, {k, 3.5}] It works for $2^1, 2^2, 2^3, 2^4,$ but ...
2
votes
2answers
2k views

Solving cubic equation for real roots

I'm looking to solve the following cubic equation for x: $\beta\, x^3 - \gamma \,x = c$. I have plugged in some sample values ($\beta = 2$, $\gamma = 5$ and $c = 2$). When I try to solve this ...
2
votes
2answers
564 views

How do I get the solutions of a cubic equation in trigonometric form?

I know $x^3-3 x+1=0$ has three roots that can be expressed in trigonometric form: $\{2\sin(10^\circ),\,-2\cos(20^\circ),\,2\cos(40^\circ)\}$. How can I get this result with Mathematica?
2
votes
1answer
144 views

What's the difference between the following two codes?

I have two codes to find the integer numbers $m$, $x$, $y$, $z$ in the system of equations. First code ...
2
votes
2answers
596 views
2
votes
1answer
199 views

Remove failed Solve[] attempts

I can't figure how to replace a failed evaluation with a null value. I am very new to Mathematica. I have an equation that are being processed using the Solve[] ...
2
votes
1answer
2k views

'not a valid variable' error when trying to solve 13 equations with 13 variables

I want to solve 13 variables by 13 equations. And here is my program: ...
2
votes
1answer
51 views

Finding an Integer, Unimodular Matrix that connects two given matrices

I have two symmetric, integer matrices, $K$ and $K_2$ which have the same determinant and the same signature (number of positive - number of negative eigenvalues). I want to find an integer valued ...
2
votes
1answer
122 views

Differential Equation with Numerically Integrated Boundary Conditions

I have a second order differential equation that I'd like to solve numerically, and then integrate its solution twice to get the parametric equations of a curve. The ODE is: $2 \kappa''(s)+\kappa^3(s)...
2
votes
1answer
48 views

Passing a list to a nested vector function

In order to apply a newton-rhapson algorithm to solve a set of nonlinear vector equations such as these, f[x_, y_] = {x^2 - 3*y + 6, (y + 1)^2 - x - 33} I need ...
2
votes
1answer
87 views

System of differential algebraic equations with one of the variables as an argument of the other varaibles

This is my first question here and I don't even know what to search for to find a solution. So I'd appreciate it even if someone can point me in the right direction. EDIT: My post had many problems, ...
2
votes
1answer
66 views
2
votes
1answer
78 views

Can not Solve FinancialDerivative for volatility

Although Solve can solve FinancialBond for an input parameter (e.g. par) given the price of a bond: ...
2
votes
1answer
132 views

How to handle a special Neumann-like boundary condition in NDSolve?

How to solve the following nonlinear ODE with two algebraic equations and one boundary condition? $$y''(x)=\dfrac{2\left((x+15)y'(x)-y(x)\right)\left(y'(x)^2+1\right)}{\left(y(x)^2+x(x+30)+236\right)...
2
votes
1answer
83 views

Can't obtain Real results using NSolve

I am trying to find the critical properties of the Carnahan-Starling equation of state, $$p=\rho R T\frac{1+b\rho /4+(b\rho/4)^2-(b\rho/4)^3}{(1-b\rho /4)^3}-a\rho^2$$ using this system of equations, ...
2
votes
1answer
102 views

Reduce not solving linear equation

I have a long equation in 7 variables but the equation is linear (affine) in x. But when I try to solve with Reduce, it keeps running forever. With Solve, I am able to get an answer. Usually I have ...
2
votes
1answer
110 views

How to plot a figure according to two equations?

There are two equations like this: $a=f(x,y)$ and $z=g(x,y)$, where $a$ is a constant. In principle, one can solve $a=f(x,y)$ for $y$ as $y=k(x)$ and substitute it into the second equation to obtain $...
2
votes
2answers
73 views

How to use the variables from Table in calculation

New to the site and Mathematica. I'm trying to use the finite element method to get a temperature profile. So if $q_0=f(T_0,T_1)$, $q_n=f(T_{n-1},T_n)$ and we also know that $T_1-T_0=T_2-T_1=T_3-T_2=\...
2
votes
1answer
186 views

How to properly interpret a seemingly contradictory result in Mathematica?

I have a problem with interpreting the results of my notebook outputs. It is important that I do not err with the interpretation. Here is the full input: ...
2
votes
1answer
202 views

NDSolve doesn't work for this simple differential equations

I try to use NDSolve to do this partial differential equations, but the result is wrong... The source term is a Gaussian pulse as time and an exponential decay with the depth z. The result should be ...
2
votes
1answer
279 views

Stability Boundaries via Routh-Hurwitz Criterion

I'm attempting to determine the stability boundaries of a 2nd order system via Routh-Hurwitz stability criterion. That is to say, I need to compute when a polynomial, which is in terms of variables "a"...
2
votes
1answer
108 views

Find tangency of InterpolatingFunction and circle

I am trying to use a system of equations to solve for the common tangent of a (polar) InterpolatingFunction and a circle. Unfortunately, both NSolve and FindRoot don't yield results, even when I "...
2
votes
2answers
198 views

How do I get all the roots within a certain range for an equation? [duplicate]

I wish to find ALL the roots of the following equation $ x \tan(ax) - c = 0$ within a range, where a and c are constants I control. I've been using ...
2
votes
1answer
348 views

Defining solutions to FindRoot as a function of parameters in Mathematica

This works: f[a_] := x /. FindRoot[x^2 - a == 0, {x, 1}] f[3] (* 1.73205 *) I need this to work, but it doesn't: ...
2
votes
1answer
159 views

Symbolic quadratic constrained maximization with non-negativity constraints?

Dear Mathematica users, I attempted to maximize the below quadratic function K, subject to x1+x2<=A, and a series of non-...
2
votes
2answers
326 views

How to find roots of a system of equation when number of equations and variables don't match?

I'm trying to check if a system of equations has a solution around a given point, the problem is that I'm using FindRoots and doesn't work when there are a different number of equations and variables, ...
2
votes
1answer
149 views

FindRoot over a domain, but give domain boundary if no root exists

I call r = x /. FindRoot[f[x] == 0, {x, 0, -5, 5}] f[x] is always monotonic, but sometimes it doesn't have a root in the ...
2
votes
3answers
778 views

Backsubstituting solution into FindRoot and Plot

I want to find roots from a transcendental equation using FindRoot as a function of a real parameter and then plot the solutions vs. the parameter. The following commands work perfectly: ...
2
votes
1answer
97 views

Strange behavior using variable in Exists

I came across a behavior that doesn't make sense to me when I tried to use Exists in Reduce. If I run commands: ...
2
votes
1answer
531 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 c[t,0]==4\...
2
votes
2answers
905 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; ...
2
votes
1answer
737 views

2d heat conduction equation: Boundary and initial conditions are inconsistent

I have the following code for a 2d heat c equation: ...
2
votes
1answer
130 views

Visualizing how the zeros for this function change as its input parameter changes?

I'm trying to explore how the zeros of a certain parameterized function change as its parameter changes. ...
2
votes
2answers
212 views

Plotting double inequality only for integer values

RegionPlot is plotting for all values. And using ListPlot I can only evaluate one inequality. How can I plot this type of ...
2
votes
1answer
476 views

Solving this system of equations produces an error message about badly conditioned matrix

I want to determine a balance distribution. Ü = {{0.4, 0.2, 0.3}, {0.3, 0.5, 0.2}, {0.3, 0.3, 0.5}}; and ...
2
votes
1answer
2k views

How to solve equations self-consistently

I want to solve the following equation self-consistently. So, H.u = e.u {{1, d}, {d, 1}}.{u1, u2} = e.{u1, u2} I guess an initial value for ...
2
votes
1answer
153 views

How to transform an expression using algebraical instead of pattern rules [duplicate]

I would like to transform rules algebraically. A very simple example would be: - k^2 - 2 k x + x^2 /. {2*k -> 1} This transforms to: - $$k^2-2 k x+x^2$$ ...