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

learn more… | top users | synonyms (4)

1
vote
1answer
93 views

Cannot Plot Output of Solve

I'm trying to plot the output of the following function t2[x_] in ...
1
vote
1answer
296 views

How to guess initial complex value for FindRoot

I have to solve a transcendental equation for a parameter, say $\beta$. Now, the $\beta$ has a range from $ik$ to $k$ where $i$ is the usual imaginary root $\sqrt{-1}$ and $k$ is a real number. ...
1
vote
1answer
292 views

Stop Mathematica from giving imaginary solutions

I have the following equation: $$D=\frac{1}{64} \pi A ^3 B \sin \left(C\right)-\frac{1}{2} \pi A B \sin \left(C\right)$$ which I want to solve for $A$. The equation is cubic in $A$ so this should ...
1
vote
2answers
816 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 ...
1
vote
1answer
179 views

Finding parameters for first positive root of solution to differential equation to be equal to a desired value

I have a differential equation with parameters. The objective is to vary the parameters such that the first positive root of the solution is equal to a desired value ...
1
vote
1answer
86 views

How to get a solution involving logs in the right form?

How do I get Solve[(1/9)^(r/2) == 1/n , r] to give me Log[n]/Log[3] and not ...
1
vote
1answer
170 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[] ...
1
vote
1answer
221 views

How to find 4 unknowns?

There are 4 unknowns constants: a,b,c,p. The variable is ...
1
vote
0answers
103 views

Nonlinear FEM and FindRoot

I'm trying to develop a kind of nonlinear FEM application using mathematica to solve a bvp like the following: $$ \gamma(u') ~u^{iv} + 2 \gamma'(u') u''' u''+ u''^3 = f(x) $$ where $u = ...
1
vote
0answers
33 views

Determine Functions Given Data Points [duplicate]

Is it possible for Mathematica to solve for f(x,y) = 0 such that {x,y} contains {{0,0.5}, {1,0.5}, {0.6,0.8}, {0.4,0.2}} ? Do we need to specify f(x,y) first ? ...
1
vote
0answers
114 views

Finding Values of All Nine Constants in q-hypergeometric Equation

I suspect that an identity of the form $$\sum_{j=-n}^{j=n}(-1)^jq^{j(5j+1)/2}=\sum_{r=0}^{n}f_{ar+bn+c}(-1)^{n-r}q^{(n-r)(dr+en+f)/2}\frac{(q;q)_{gr+hn+j}}{(q;q)_{n-r}}$$ for $a, \dots , j$ integer ...
1
vote
0answers
38 views

RSolve works slow and does not give an analytical solution for solving high order linear difference equations

RSolve@@{{a[n]-a[n-10]==101}~Join~Table[a[i+1]==101+10i,{i,0,9}],a[n],n} I want to get the exact form of general formulas, but ...
1
vote
0answers
114 views

Solving Thue equations

A Thue equation is a 2-variable homogeneous integer polynomial of degree at least 3. It's well-known that such equations have only finitely many solutions over the integers. I'm trying to solve some ...
1
vote
0answers
23 views

Minimizing roots of a transcendental equations [duplicate]

I'm trying to use the solution of FindRoot with parameters for further calculations. Annoyingly enough it always gives me trouble with replacing the parameters. As an example I'm attaching a simple ...
1
vote
0answers
30 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
125 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
29 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
1answer
43 views

How do assign an interval of values to a variable in order to use them in a function?

From my code (below) I have K1 and K2 set to zero in order to test this numerical method. However, K1 is on an interval from 0 to 32, and K2 is on an interval from 0 to 8. I would like to know if ...
1
vote
0answers
19 views

Solving a pair of recurrence relations [duplicate]

Why does this RSolve[{y[t] == a*y[t - 1]}, {y[t]}, t] work, while this ...
1
vote
0answers
60 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
80 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
1answer
231 views

Solving a System of four equations (with logs) [closed]

I am trying to solve a system of four equations in four variables. I have read a number of threads on similar issues and tried to follow the suggestions. But I think it is a bit messy because of the ...
1
vote
0answers
132 views

Calculating derivatives of roots of equation

I am interested in calculating how the roots of an equation change as one of the variables does, when the equation is evaluated numerically rather than symbolically. I have edited this code to be ...
1
vote
0answers
57 views

Gram Schmidt Process for functions [duplicate]

I want to implement the Gram Schimdt procedure to the space of functions in the following form \begin{align}‎ ‎R(y,x)=\left\{‎ ‎\begin{array}{ll}‎ f(x,y), ‎& y\leq x\\ \\‎ g(x,y), ‎& y>x‎ ...
1
vote
1answer
211 views

Unable to solve a set of equations with Solve

I'm trying to solve the following equation using Mathematica 9.0.1.0 : ...
1
vote
0answers
163 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
91 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
229 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
201 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
233 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
80 views

Simple Solve example fails: why? [duplicate]

Solve is an amazing tool, but sometimes it doesn't work in the way I expect it to. The particular instance that is confusing me right now is this: ...
1
vote
1answer
197 views

Gram Schmidt inner and outer products

I know the Gram-Schmidt orthogonalization generates an orthonormal basis from an arbitrary basis. I need help with: Write a program that inputs a list $\{b_1,\dotsc,b_n\}$ of linearly independent ...
1
vote
0answers
37 views

How can I access numerical values of variables calculated with Solve? [duplicate]

I am trying to calculate some thermodynamical cycles so I often use lists with equations. Here's a simple example: ...
1
vote
0answers
101 views
1
vote
0answers
170 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
65 views

Is it possible for Mathematica to give the step by step process of Reduce command? [duplicate]

Input: Reduce[3 x^2 - 2 a x + a^2 - 1 >= 0 && x > a, x] // TraditionalForm Output: ...
1
vote
0answers
132 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
224 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
59 views

Finding Maximum Root [duplicate]

How can I find the greatest root for a polynomial if the polynomial degree is very large. One way is using NSolve. This command can find all roots in numeric. Then, ...
1
vote
0answers
88 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
194 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
543 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
94 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
3answers
171 views

Solving One Equation with Two Variables

I am trying to solve an equation with two variables. It is the last step in the process of using the method of undetermined coefficients to solve a nonhomogeneous differential equation. The equation ...
0
votes
2answers
144 views

Unexpected complex roots from solving a cubic

I am new to Mathematica and there is a behavior that I don't understand. When I insert the code: Solve[-(x - 2)^3 + 10 x - 10 == 0, x] I obtain: ...
0
votes
1answer
533 views

Solving linear equation

I want to solve following linear equations: \begin{align*} 3 x - y &= 0\\ y + 2 z &= -1\\ x + y + z &= 3 \end{align*} The output will be ...
0
votes
2answers
70 views

Solve prints message saying that it cannot solve my equation [closed]

I want to solve the fallowing equation: Solve[-(2/x^3) - 1/x^(5/2) - 1/x^2 + (4 Log[Sqrt[E] x])/x^3 == 0, x] Mathematica complaints: This system cannot be ...
0
votes
3answers
118 views

How to solve correctly inside a DO loop and then plot

Here is the corresponding code First the equations setup ...
0
votes
2answers
89 views

Gaps in plots of equation solution from parameter

I'm trying to plot solution of system of equations as a function of parameter ...
0
votes
2answers
193 views

Constraining two points, a specified distance apart, to a line

Suppose that I have two points in the xy plane: pt1 and pt2init. pt1 is fixed in space, ...