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

learn more… | top users | synonyms (4)

7
votes
0answers
362 views

Creating and using an explicit piecewise function in a convenient way

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

What Method options does FindInstance accept?

While answering Trouble getting FindInstance to return multiple results for certain constraints it became clear that for this problem FindInstance would be better ...
5
votes
0answers
72 views

Using quantifiers with free indices as declarative conditions in equation solving

This is a toy example: Is there a way to rewrite ...
5
votes
0answers
79 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 ...
4
votes
0answers
653 views

Solve huge symbolic system of linear equations

I have a system of 76 symbolic linear equations (i.e. some coefficients are symbolic) with a sparse coefficient matrix. However, neither Solve[] nor ...
4
votes
0answers
120 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( ...
3
votes
0answers
81 views

Solving huge equations

I have this expression: $d=-\frac{\sqrt{\frac{780 (f+130) \left(5 d^2 (13 c+f+130)-5 d (c (f+1300)-325 (f+160))+125 \left(2 (2 c+725) f+5 f^2+104000\right)+13 d^3\right) \left(5 d^2 (52 c+17 ...
3
votes
0answers
126 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 ...
3
votes
0answers
285 views

FindRoot gives a wrong solution which obviously should not be there

I got stuck on FindRoot and I didn't see any similar problem posted, so let me explain what I am trying to do and what problem I meet here. I try to find roots of a particular function, which in the ...
3
votes
0answers
101 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
0answers
135 views

RSolve not returning known solution

I have the following recurrence relationship for which I know the solution is given by a binomial distribution: ...
3
votes
0answers
773 views

Solving a system of nonlinear equations self-consistently

I am trying to solve a set of three non-linear equations in Mathematica. I need help with them. The Mathematica code (in plain text format) is attached below. If you copy & paste the code below ...
2
votes
0answers
70 views

Deciding whether a semi-algebraic set is empty. Which command is best?

[Update] I was able to prove that the "larger" set $ 1\ge x \ge y,w \ge 0$ and $p(x,y,w)>0 $ is empty with another computer algebra system in my slow laptop in few minutes but in Mathematica even ...
2
votes
0answers
112 views

Question about Solve

I don't understand what Mathematica is doing to get these solutions. I ask it to solve two equations, the left hand side of both expressions are the same, the right hand sides differ by a ± sign. I ...
2
votes
0answers
93 views

Equation Solving: Solution in terms of predefined functions

Is there a way that I can have Mathematica give me the solution of a problem it has solved in terms of functions that I have predefined? The specific problem I'm trying to solve is: ...
2
votes
0answers
141 views

Using WhenEvent for derivative of discontinuous function

I have a discontinuous function ($u(t)$, a square wave) and I would like WhenEvent to trigger when the signal goes high/low, i.e. when the value of $u(t)$ changes. ...
2
votes
0answers
160 views

Problems with Simultaneous Equation Solving

Let me first explain what I am doing and then I will explain the problem that I can't figure out. I am basically trying to solve 2 equations for 2 unknowns (a1 and ...
2
votes
0answers
253 views

I can't solve this set of equations

I am having trouble solving a set of equations in mathematica. I do not get any errors, but i don't get any answers as well. I have some restrictions on my parameters and variables in the equations ...
2
votes
0answers
126 views

Nest for large value of n

I am working on a problem to find the Feigenbaum constant. In this problem I have to get $f^{2^n}(1/2)$ with $f(x) = kx(1-x)$. I use Nest to get $f^{2^n}(1/2)$ but ...
2
votes
0answers
198 views

Conditional statements in intial conditions?

This is potentially a daft question, but I thought I'd ask it; I have some material free to diffuse in a boundary between rn and ro; I've been able to get it working nicely for neumann type boundary ...
2
votes
0answers
216 views

differential equations with implicit functions

I have an ordinary differential equation that contains coefficient functions that depend implicitly on the independent variable via an algebraic equation. I am trying to go ahead and use NSolve to ...
2
votes
0answers
124 views

Why is FindRoot initial value far from the specified one?

I am trying to numerically find the root of a function that looks a bit like: 1/x - (SchurDecomposition[A[x]][[2]])[[1]], where ...
2
votes
0answers
171 views

Solving $k\log(n) - k\log(k) \leq u$

How can I persuade Mathematica to help me find solutions to the inequality $k\ln(n) - k\ln(k) \leq u$ over the reals, assuming $u \geq n \geq k \geq 1$. When I evaluate ...
2
votes
0answers
164 views

Am I entering this 2-variable recurrence correctly?

I'm trying to solve this recurrence with two variables. ...
2
votes
0answers
226 views

FindRoot - Speed and errors

I am using FindRoot[] to solve a complicated equation. It seems I get the correct answer even though I get errors about an ...
2
votes
0answers
289 views

Numerically/Analytically Solving a System of Equations

I have $6$ functions $f_i(x,y,z)$, $(i = 1, \ldots, 6)$ in three variables $x,y,z$, and I would like to find a simultaneous instance of these variables, say $(x_0, y_0, z_0)$, such that $f_i(x_0, y_0, ...
1
vote
0answers
23 views

Correct use of Mathlink's MLGetRealList

I am attempting to use Mathlink from C. I have managed to return single values from mathlink to my c program, but I am having trouble with lists. Consider the Mathematica operation: ...
1
vote
0answers
65 views

FindRoot in range

I am trying to find all roots in a range, using y /. FindRoot[Sin[y] == 0, {y, #}] & /@ Range[7] which gives ...
1
vote
0answers
25 views

How can I select two lists from a list of lists so that I receive the system of equations has only rational solutions?

I want to make a sytem of equations has the form $$\begin{cases} a_1x^2+b_1y^2 + c_1xy+d_1x + e_1y+3=0,\\ a_2x^2+b_2y^2 + c_2xy+d_2x + e_2y+3=0. \end{cases}$$ so that the given system of equations has ...
1
vote
0answers
50 views

Finding undetermined coefficients with boundary conditions

I have 2 equations which are; To[x_] := c1 E^(-a1 x)/2 + c2 E^(a1 x)/2 + a2 Tu[x_] := c3 E^(-a3 x)/2 + c4 E^(a3 x)/2 Tu equals almost zero ...
1
vote
0answers
31 views

Equation involving hypergeometric functions

I want to solve this equation but nor NSolve nor Solve are able to do this. ...
1
vote
0answers
54 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 ...
1
vote
0answers
93 views

Solve[] unexpectedly giving empty solution sets for extremely large problems

I was calculating the set of lattice vectors in $\mathbb{Z}^4$ which have taxicab norm $r$. For example, executing ...
1
vote
0answers
66 views

No response from DSolve

Using optimal control theory, I am trying to find the optimal paths of d and k such that they will maximize profit as given by lc: ...
1
vote
0answers
205 views

Differential equations: How do I solve a delay differential equation with 0-1 variables

How do I solve the following non-linear, delay-differential equations using Mathematica under the assumptions $$p(t) = 1 ,\; y(t) > k$$ $$p(t) = 0,\; y(t)<= \ k$$. $$ \frac{dw}{dt} = ...
1
vote
0answers
161 views

Why DSolve and NDSolve give different answer

Consider these code, which solve the same equation using DSolve and NDSolve, why do they give different answer? Which one is ...
1
vote
0answers
196 views

Solving asymptotic equations

Let's say I have several equations that have solutions given as a negative series $$y(t,x)=\sum_{n=-1}^\infty a_n(t) x^{-n}=a_{-1}x+a_0+a_1x^{-1}+\mathcal{O}(x^{-2})$$ My equations look like ...
1
vote
0answers
92 views
1
vote
0answers
165 views

Having Mathematica solve or help solve a set of 8 equations

There is a set of DDE's that arise in leukemia dynamics and I'd like to find the equilibrium points of the system by setting all of the derivatives equal to zero. The system can be described in ...
1
vote
0answers
107 views

Making simultaneous equation solving more efficient

I am trying to solve a system of several simultaneous equations that have been generated by a series of modules and functions, using this module: ...
1
vote
0answers
197 views

Why is a bump function making NDSolve take forever to solve?

I am attempting to solve a system of relatively complicated elliptic PDEs via a relaxation technique. In particular, I am trying to solve $\textrm{div} LW = S$ where $S$ is some vector and $L$ is the ...
1
vote
0answers
86 views

Solving for vector elements in a function to which the vectors are not passed explicitly

I want to define a function which takes in two integers (indicating the lengths of 2 vectors), and solves a simple set of expressions at a set of points to find all the values in both the vectors. so ...
1
vote
0answers
170 views

Automatic inverse for Lambert’s Azimuthal Equal-Area Projection

I’m interested in finding closed-form inverses for map projection equations in general (of course, this is not always possible). I’m hoping Mathematica can help me with some of the trickier ones. As ...
1
vote
0answers
445 views

Calculating the mean curvature of a surface - suggestions

I am attempting to calculate the mean curvature [1, 2] of a surface defined by a function of x and y. My function is rather ...
1
vote
0answers
89 views

Symbolic FindInstance

I'm wondering if there is any way of doing something like the following: FindInstance[x < y, {x}, Reals] and getting, as a possible return value, ...
0
votes
0answers
9 views

What is the efficient way to code Successive Over-relaxation (SOR) method in Mathematica?

I have written a SOR method code using C-style procedural loops. However, I think there might be (much) better ways to achieve the same end in Mma avoiding these loops. So my question is: Is that ...
0
votes
0answers
48 views

Equation in implicit form

This is simple example: v = {1.48622576, 1.53444276, 1.42299438} v/(v + 1) Out[1]:= {0.597784, 0.605436, 0.587288} In the above case expression v/(v + 1) is ...
0
votes
0answers
26 views

Minimize solution to linear system with singular system matrix

I am trying to solve a time-dependent linear system with a singular system matrix, $A$, thus an infinite number of solutions, ${\bf x}$(t), exist that satisfy the linear system. The solution I am ...
0
votes
0answers
43 views

Solving an exponential inequality with parameters

Can I use Mathematica for solving the inequality given below? b^(x/2 - 3/2) - b^(ax) - 2c > 0 where b, a, c are the parameters. I want to find out for what ...
0
votes
0answers
39 views

Solving a single equation with multiple parameters

If I wanted to solve equations like Cos[x]^2 = A Cos[Bx + C] + D what function should I use? The Solve function seems to only solve for as many variables as there are equations. Please note that my ...