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

learn more… | top users | synonyms (5)

1
vote
2answers
37 views

FindRoot with Grad

I want to find a numeric solution using Grad. With Solve, I get ...
2
votes
1answer
206 views

Is there a package to find ALL exact roots of a polynomial, if they exist?

There are polynomials with roots not expressible with radicals but expressible as trigonometric or other functions, for which Solve[] only returns ...
1
vote
1answer
29 views

REAL Roots of a 2 variable Polynomial

I am looking to plot the real zeros of polynomials of two variables in Mathematica, but I have encountered an issue. To accomplish this task I have been using ContourPlot to plot the polynonial ...
0
votes
1answer
59 views

Solve homogeneous system [on hold]

Let x^2*y + x^2*z + x*y^2 + y^2*z + xz^2 + yz^2+ xyz =0 x^3-xyz + x^2*y-xyz + x^2*z-xyz + xy^2-xyz - xy*z+y^2*z - xyz+xz^2 - xyz+yz^2 =0 x^3+3*xyz+7*yz^2 + xy*z+y^3+9*y*z^2 + x^2*y+6*xyz+4*y*z^2 ...
0
votes
1answer
49 views

Solving for Polynomial roots

This simple Solve gives the roots of a quadratic: Solve[a x^2 + b x + c == 0, x] However, if I factor the polynomial in terms ...
33
votes
8answers
7k views

About multi-root search in Mathematica for transcendental equations

I have some questions for multiroot search for transcendental equations. Is there any clever solution to find all the roots for a transcendental equation in a specific range? Perhaps ...
0
votes
0answers
54 views

Equation cannot be solved with the methods available to Solve [duplicate]

I am a beginner to Mathematica. I am trying to get the solution of such equation as: Solve[-b*M + L*R - a*M*Log10[L*c + (L*d*(1 - M)*S)/R] == 0, R] Constants are ...
0
votes
0answers
34 views

Create a loop for assignment of Solve output [duplicate]

Is there a way to automatically assign the output of solve to a number of symbols? For example, assume that I have the solve statement as follows: ...
0
votes
0answers
38 views

The ` in Solve results [duplicate]

I assign the result of a Solve to the variable eq8. Then I attempt to use this result in the definition of a function as follows: ...
2
votes
2answers
44 views

Difficulty finding roots of an interpolation function from NDSolve

I have been trying to find the point at which one of the solutions of a system of two ODEs crosses zero. I used the method suggested in this answer to a previous question, which seemed to be the most ...
0
votes
0answers
44 views

How to make sure that Mathematica doesn't hog memory for intermediate steps while using functions like FindRoot? [on hold]

I am using Mathematica 10.4. I am trying to solve a set of around 10000 non-linear algebraic equations that I have computed upon discretization of a set of coupled non-linear differential equations. ...
4
votes
1answer
79 views

Trying to create an interactive plot (newbie here)

I'm new at Mathematica so please bear with me. I am trying to create an interactive plot that will end up looking like this: I want to input values for certain variables and then plot a curve. This ...
0
votes
2answers
58 views

Using result of Solve as plot parameter [on hold]

Recently I have faced a problem I would appreciate it, if you helped me I have the following code in Mathematica I want to plot a function let say z * 10 * t where ...
0
votes
0answers
17 views

Schmidt decomposition problem [migrated]

I'm having a problem in implementing the following problem: I have a quantum state so defined: $\left| \Psi\right>=\int ...
1
vote
1answer
40 views

How to define a Reduce function in order to define two parameters afterwards

I'm trying to evaluate numerically an equilibrium for a system of ODEs: dA = c A + b B - a A dB = d B + a A - b B So I define a list of known parameters: ...
3
votes
1answer
60 views

Solve with Reals gives different answer

I'm very new to Mathematica (started using ver 10.4) and I was just verifying my calculations. I have the following equation: ...
0
votes
1answer
68 views

Solution for Imaginary part of equation

I have a function Il12i[-0.001, -0.001, w, 0, 1, 0, 0.1002] which is depends on 1 variable w and it is real if w > -3666 (this number I found from graphs). This function is very complicated numerical ...
0
votes
0answers
3 views

Find all solutions to the equation over field Z_16 [migrated]

So for the first part of the question, I have to find the solutions to: $\ x^2+4x+3=0$ over $\ Z_(16)$ I have found these to be X=5,7,13 and 15, just by using standard method to solve a quadratic ...
1
vote
1answer
53 views

Testing equivalence of analytical expressions like $x^2 -x == x(x-1)$

On occasion I have two expressions that I may be equal but simplified differently, and I want to check if they are actually equivalent. Using the equation in the title, $x^2 -x == x(x-1)$, sometimes ...
0
votes
1answer
501 views

Problem with FindRoot applied to functions

I am having difficulty with the error "is not a list of numbers with dimensions..." when using FindRoot (and other numerical routines in Mathematica) to solve equations numerically when the argument ...
0
votes
2answers
37 views

Solving several inequalities in Mathematica

I am a very beginner of Mathematica and the following will probably be a very stupid question. I have a bunch of inequalities in 6 unknowns and I want to find the set of values satisfying them. Can I ...
0
votes
2answers
70 views

FindInstance is very slow on an apparently simple problem

I wish to use FindInstance to find a set of solutions of this problem: ...
4
votes
2answers
823 views

Solving a linear equation in Mathematica

This should be easy but I can't seem to find the right way to do it. I have an equation of the form $a x + b x + c y + a z + d z = 0$, and I'd like to solve for relations between the parameters ...
6
votes
2answers
245 views

Find a condition that b must satisfy so that Ax=b has solution

I'm new to Mathematica, so I'm sorry if this is really simple. I am trying to find the condition that vector b must satisfy so that Ax=b has solution. I would like to learn a general method, but I'll ...
1
vote
1answer
227 views

Why Integral does not converge?

I am trying to solve for $\Omega$ this nonlinear integral equation: $$1+\dfrac{z}{k^2}-\dfrac{z^2}{K_2(z)} \dfrac{\Omega}{k^3} \displaystyle\int_{1}^{\infty} \gamma^2\, \text{ArcTanh} ...
0
votes
1answer
48 views

Errors only when code is iterated

Hello I'm trying to have two for loops sweep two parameters and for each 2-tuple of those parameters I want to solve a differential equation. To do this I have the following code: ...
5
votes
4answers
169 views

Solving Recurrence Equation $a_{n+1}=a_n(1-a_n)$

How to solve following recurrence equation $$a_0=2$$ $$a_{n+1}=a_n(1-a_n)$$ I tried: ...
0
votes
1answer
50 views

Reduce unable to solve simple inequalities with quantifiers

I have try to solve following inequalities using Mamthematica. But it unable to solve that. ...
0
votes
0answers
27 views

Finding results using qualifiers?

I have set of recursive function definitions and assumptions $\forall n.t1(n+1)=t1(n)-Y*S(n)$ $\forall n.t2(n+1)=t2(n)+S(n)$ and assumption $\forall n.n<K\rightarrow L1=L2*t2(n)+t1(n)$ I want ...
4
votes
1answer
1k views

How to take the derivative w.r.t. an arbitrary function?

EDIT: The answer posted by @Jens works for built-in functions only. However I'm starting to wonder if this is related to my Mathematica install somehow. I have version 7. Here is an example in the ...
1
vote
1answer
75 views

Solving two equations with modified Bessel functions

I am trying to solve two equations with Bessel functions in them, 1) C1*BesselK[0, 3.7268*10^-4*x] == 1.3*10^-6 2) ...
2
votes
1answer
49 views

Substituting values of variables in solution of a system of linear equations

After solving a system of 21 linear equations in 21 variables (called Xi in my code) I got the values for all the 21 variables, in the following way: ...
-1
votes
0answers
37 views

Solving integral equation

I want to solve the following integral equation for $k$: $\left[\frac{\int_0^\infty e^{-ay}(1+y)^{C} dy}{\int_0^\infty e^{-ay}(1+y)^{C-1} dy}\right]-1=\int_0^\infty e^{-ay}(1+y)^{k} dy$ where $C$ ...
1
vote
0answers
23 views

Help Compiling Vectors using Iterations

So I have a nice code that runs properly, which I include below but, in my naiveté, I'm having issues iterating my procedure. Simply put, I have a function of two real variables $f(x,y)$ and I want ...
3
votes
1answer
73 views

Solve equation for a function as solution

I need to solve the following function: $$ h(z) = \sum^{40}_{k=0} \frac{k \, p(k)}{z- k \, h(z)} $$ where $p(k)$ is a probability distribution (in this case binomial with n = 40). Note that $h(z)$ is ...
0
votes
1answer
54 views

Arranging sides of an equation [duplicate]

I have to create my own function to solve equations (first and second degree). Assuming you have any input equation, how can I do to bring everything to the first side of the equation? Example. I ...
0
votes
0answers
36 views

Matrix with equations as entries - how to choose them and solve them

I have a matrix in Mathematica, which I have built using Normal[SparseArray]. Each entry of this matrix is an equation in one variable. What I need to do is 1) to select some of the entries of the ...
0
votes
1answer
31 views

Parameterize Subset of Complicated Zero Set

I recently asked this question: Evaluating Functions on Zero Sets (of other functions), and got fantastic assistance, but realized that there seems to be an extra issue when you're trying to ...
0
votes
1answer
53 views

Module output rules like Solve FindRoot etc built-in functions

How do I print rules {x -> values} as an output from my module? I'd like it formatted just like those from the built-in functions, so that users can use Out[ ] or % to directly apply the rules. ...
3
votes
1answer
208 views

How to change parameters to equations of physical properties they belong to in an equation?

First off, I want to apologize for my poor English as I am not a native speaker. In my physics research, I had to solve a system of equations, from which I would get a solution for the speed of an ...
0
votes
0answers
24 views

Using the “reduce” function with s lot of equations and inequalities

I am using "reduce" to obtain all positive integers (plus 0) which satisfy a group of equations and inequalities. I have to automatize this procedure because I have to use "reduce" with at least 100 ...
3
votes
2answers
103 views

Evaluating Functions on Zero Sets (of other functions)

I have a highly non-trivial function, call it $f(x,y)$, in which I want to look at its zero set in the $(x,y)$ plane, and plot/evaluate a different function $g$ precisely along this zero set. I have ...
1
vote
2answers
108 views

Find all the non-negative integer solutions of a linear equation

For the linear equation x[1] + x[2] + x[3] + x[4] + x[5] == n Is there any efficient way to find all the non-negative integer solutions of this equation when n ...
0
votes
1answer
87 views

Solve equation with two different sets of replacement rules

I want to solve eqn for c0 in the following code. I want to do it with both replacement rules. I realize I could just use solve ...
0
votes
3answers
85 views

Interactive Application using Manipulate for Fibonacci Series

I am working on creating interactive app using Manipulate function on fibonacci series using SeriesCoefficient. I want to plot the fibonacci numbers based on the slider of Manipulate function. I ...
14
votes
3answers
5k views

Solving a time-dependent Schrödinger equation

I want to solve the time-dependent Schrödinger equation: $$ i\partial_t \psi(t) = H(t)\psi(t)$$ for matrix, time-dependent $H(t)$ and vector $\psi$. What is an efficient way of doing this so that ...
5
votes
0answers
48 views

Is there a way to tell how much time FindRoot is spending calculating Jacobians?

I am solving a large number of trascendental equations by scanning FindRoot over a bunch of seeds and a bunch of different parameter values, and I would like to ...