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
1answer
76 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
96 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
109 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 ...
2
votes
1answer
184 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
196 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
242 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 ...
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
191 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
156 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 ...
2
votes
2answers
291 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
144 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
726 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
96 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
519 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 ...
2
votes
2answers
875 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
707 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
211 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
422 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
151 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$$ ...
2
votes
2answers
408 views

Why can't FindInstance find a solution to this equation?

FindInstance[680 == 1/3 2^m (-1 + 4^v), {m, v}, Integers] FindInstance seems to be unable to find the integer solution $ {m = ...
2
votes
1answer
673 views

Defining regions for `Exclusions`

I have the following function I want to plot in 3D (For simplicities sake we can say $a=5,b=-1,c=3$.) ...
2
votes
1answer
325 views

How to solve a symbolic equation with symbols in powers

I am trying to solve $\dfrac{\rho\ ((1-b)\ \lambda _f\ \lambda _s\ P _A)^{\rho }}{b-1}+\dfrac{\rho \left(\dfrac{b\ \text{ps}\ x}{1-\text{ps}}\right)^{\rho }}{b}=0$ in terms of $b$, I ...
2
votes
1answer
88 views

Putting as many expressions as possible equal to zero

Say I have a list of expressions involving many different variables. I would like to see how many of the expressions I can set equal to zero thereby solving for some of the variables in terms of the ...
2
votes
1answer
203 views

Constant Heat Flux Boundary Condition for the Differential Heat Equation

I am attempting to solve the transient heat equation for a very simple geometry and I am getting a junk solution. I have identified two distinct regions (wall and super-heated fluid layer) with their ...
2
votes
1answer
102 views

Solving a consistent overdetermined linear system

Not a duplicate question. I already tried all the Answers offered for this Question and they failed, so I believe this is a different problem with Solve[]. I have ...
2
votes
1answer
192 views

Solving SDE: $\frac{dy(t)}{dt}=(c+\sigma_wW(t))y(t)+\epsilon(t) $ in Mathematica

I want to solve this differential equation $\frac{dy(t)}{dt}=(c+\sigma_w W(t))y(t)+\epsilon(t) $. For details see ...
2
votes
1answer
120 views

Solve an inequality involving a sum with a parameter

Good day. I am new to Mathematica and I am looking for advice. Is it possible to solve $$\sum_{k=0}^{n} \frac{p^k}{k!}>0$$ for $n$ where $p$ is a parameter? When I do ...
2
votes
2answers
158 views

Solving a system of weird equations

I'm trying to have Mathematica solve a system of equations for a friend's work in Computer Vision. I'm essentially going in blind to this and am not very familiar with Mathematica syntax, so would ...
2
votes
1answer
97 views

Reduce / Solve cannot solve simple equation

Why does Reduce[{a^b*c^b == c, 0 < b < 1, a > 0, c > 0}, c, Reals] give ...
2
votes
1answer
100 views

Can't NSolve system of exponential and linear equations in several variables

At the moment I am stuck with the following code, when specified to use root , quadratic or log utility functions all it fine, only exponentials like: $\small \left\{-e^{-c_{0,1}}-0.25 ...
2
votes
1answer
96 views

How to solve for the closed form of this recurrence?

For integers $L,k \geq 0$: $$T(L,k) = \begin{cases} 0 & \text{if }k \geq d, \text{ever}\\ -a_{k} & \text{if }L=0 \text{ and } k<d\\ a_dT(L-1,k+1) - a_{k}T(L-1,0)& \text{if } L>0 ...
2
votes
1answer
143 views

Fixed Points Problem with 2D Mapping

I have been playing with the mapping given by ...
2
votes
1answer
246 views

How to figure out whether a nonlinear inequality holds or not

I have an inequality consisted only of seven parameters as follows: $(1 + g)^{-(m+w)}\big[\big\{a+(1-a)(1 + g)^m\big\}(1 + p q)Q - (1 + g)^w \big\{(1 + g)^m (1 + q)-(1-p)q\big\} Q\big] < 0$ ...
2
votes
1answer
225 views

How to solve three equations with four unknowns iteratively?

I have three equations of the form $\quad \eta = f(N1L, N2L)$ $\quad N1L = g(\eta, R)$ $\quad N2L = h(\eta, R)$ I want to eliminate $\eta$, so that I can get $N1L$ and $N2L$ for a range of $R$. I ...
2
votes
1answer
142 views

NSolve for system of polynomial equations

I would like to find all the solutions to the following system of 8 polynomial equations for the variables w1, w2, w3, w4, w5, w6, w7, w8 ...
2
votes
1answer
62 views

Solve[ ] for 2 variables not returning list of two variables

I type into mathematica: Solve[2 Sin[(Sqrt[3] x - 3 y)/2] == 2 Sin[(Sqrt[3] x + 3 y)/2], {x, y}] expecting it to return a list of pairs of solutions of the form ...
2
votes
1answer
148 views

Equation involving hypergeometric functions

I want to solve this equation but nor NSolve nor Solve are able to do this. ...
2
votes
1answer
174 views

Deriving least-squares equations in Mathematica

I define $Y = a + b*x + e$, and I want to find values for $a,b,$ by least squares, so I minimize $RSS = sum((Y_i - a - b*x_i)^2)$ I then take derivatives, ${d RSS \over da} = 0$ and ${d RSS \over db} ...
2
votes
1answer
222 views

Simplify equation using conditions/inequalities

If I have an equation like f[x](1 + a)==0 I'd like to be able to tell Mathematica that 1+a doesn't vanish in general, and get back ...
2
votes
1answer
139 views

Declaration of variables in large Linear Programming model with NMaximize

I want to solve a large LP problem using NMaximize. This is a small instance of my bigger problem: $$max \sum^{4}_{t=1}(\sum^{4}_{c=1} V_{c}*f_{c,t})\\ s.t\quad ...
2
votes
1answer
122 views

InterpolatingFunction with growing Domain

I want to find roots of a function that is very slow to calculate (the function itself involves a bunch of FindMaxima) so I interpolate the function in a region where I think it likely the the root ...
2
votes
1answer
158 views

Differentiating with respect to vectors of unspecified length

Given a function defined by f[x_,y_] := Sum[ x[[i]] (g[y[[i]]] - g[Total[y]] - Log[x[[i]]]), {i, 1, Length[x]}] I can easily calculate the derivative for a fixed ...
2
votes
1answer
1k views

NDSolve two PDE plus one ODE, together

I have a couple of PDE (time and space/length) that are correctly managed by NDSolve, but to which I’m adding an ODE on which the acc “accumulates”, through time, a ...
2
votes
1answer
658 views

How can I solve an exponential inequality in Mathematica?

I would like to solve the equation: A Exp[m t + s Sqrt[t] x - 1/2 s^2 t] >= K with respect to x (which is in the exponent) ...
2
votes
1answer
549 views

Inconsistent initial and boundary conditions

I have a 4th order nonlinear PDE with the following BCs and IC: ...
2
votes
1answer
656 views

Reduce/Solve an equation with symbols in powers

I am trying to solve the following equation with respect to x. eqn = x^z + y == a; But when I do ...
2
votes
0answers
67 views

Complicated implicit plotting - can it be tackled?

I'm quite new to Mathematica and have already spent tons of time on problems of this type, browsing through tutorials and topics here. I'm now coming to conclusion that it cannot be done and am about ...