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

learn more… | top users | synonyms (5)

3
votes
1answer
187 views

Finding closest point to implicitly defined curve

I'm facing the following problem: Given an implicit curve given by $f(x,y)=0$ and an initial guess $(x_0,y_0)$ which is hopefully close to the curve, I want to find the point $(x,y)$ on the curve ...
3
votes
1answer
587 views

What algorithms does NSolve use?

I've been using NSolve to find roots and now I've been wondering what underlying algorithm Mathematica uses. It's probabaly similiar to ...
23
votes
2answers
1k views

Using the Krylov method for Solve: Speeding up a SparseArray calculation

I'm trying to implement this Total Variation Regularized Numerical Differentiation (TVDiff) code in MMA (which I found through this SO answer): essentially I want to differentiate noisy data. The full ...
3
votes
1answer
212 views

Faster Solve for Fermat 4n+1 conjecture

Assuming that Fermat 4n+1 conjecture (each prime of the form 4n+1 is the sum of two squares) is true then I like to solve the ...
2
votes
1answer
519 views

Solving a transcendental equation

I want to solve the a transcendental equation, e.g. NSolve[-(1-x)*Log[1-x] - x*Log[x] == 0.5, x] Mathematica is not able to solve, wheres WolframAlpha -- is able ...
1
vote
1answer
86 views

Reduce Inequality

How can I reduce the inequality in terms of modulus of alpha. ...
2
votes
1answer
101 views

Implicit derivatives in one line. How do I do it?

I have this little script to derive implicitly $PrePrint = # /. {D[y_, x_, NonConstants -> {y_}] :> y'[x]} &; Example: ...
1
vote
0answers
168 views

How to solve a nonlinear overdetermined system of equations

i'm trying to solve a nonlinear overdetermined system of equations (124 equations, 123 variables) i tried this procedure (newton raphson) ...
2
votes
2answers
235 views

Finding integral points on a surface

Suppose I have a dimension formula (for a Lie algebra representation) that is as follows: $$ d(a,b) = {(a+1)(b+1)(a+b+2) \over 2} $$ Now consider the surface $F(a,b,n) = 0 = d(a,b) -n$ where $n \in ...
22
votes
6answers
6k views

Finding the intersection of a curve with an interpolation function

There are no issues trying to find the intersection points of two defined curves. ...
0
votes
2answers
123 views

Significantly different results between Mathematica and Wolfram Alpha with large radicals and small numbers [closed]

I have some code which involves tiny numbers being put to the power of very large numbers. The function I'm looking at is $\varphi = \omega(T) \left(1 - (1 - \epsilon)^{n_{e}(T)} \right)$ when ...
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
102 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 ...
4
votes
1answer
311 views

Numerical errors/inaccuracies in ProductLog

Context In cosmology, a fairly accurate model to describe the gravitational potential, $\psi(r)$ of dark matter halos is given by $\psi( r)=\log(1+r)/r$. ...
0
votes
1answer
67 views

Removing imaginary and negative results from Solve

I have been trying to solve a series of equations that I know have real positive values for the unknowns (they all have real world analogues and they are not linked to a frequency domain) but despite ...
15
votes
5answers
992 views

How to make a Line[] with no end?

I'm trying to do this: In this graph, the secant points are aproximated in order to become the tangent, it seems I need some kind of function which plots a line based on two points and it's points ...
1
vote
0answers
52 views

DSolve -> no results

I am trying to solve a system of differential equations. ...
0
votes
2answers
68 views

NSolve not solving [closed]

I've been working on a Newton Cooling equation problem for my diff equ class and cannot get NSolve to solve this equation. I've plugged it into Wolfram Alpha and my ...
0
votes
1answer
73 views

Elements inside arrays seem to be strings rather than signed numbers

When I have code like b = Array[Subscript[be, #] &, {3, 1}] b = b /. {Subscript[be, k_] :> Subscript[-be, k] /; k > 1} I am noticing that Mathematica ...
5
votes
2answers
284 views

When does the real part of Zeta vanish on the critical line?

This seems to be a quite a simple problem but I cannot make it work. I am trying to find all values within a given range for which the real part of the Zeta ...
4
votes
1answer
544 views

3D Plot of The Riemann Zeta Function

I am a new user and i would appreciate your help to visualize the 3D of the following function. Edit I tried this ...
2
votes
2answers
171 views

Solving for any one variable in a formula given values for all the others

Suppose I have an equation relating a number of variables and I know the values of all but one of them. Using the equation and known values as input, is there a simple command (built-in or otherwise) ...
0
votes
1answer
478 views

How to use `FindRoot` to solve an equation containing a parameter?

I'm trying to derive some of the results of the following paper: Electrodynamics of semiconductor-coated noble metal nanoshells, JT Manassah - Physical Review A In the paper there is matrix $\mathbf ...
0
votes
1answer
248 views

How to calculate the unknown quantity in an infinite series?

I'd like to calculate x value in this equation. Basically, I tried to 2 types of method which are FindRoot and NSolve. But, I have failed the calculation caused by these errors up to now. If ...
3
votes
1answer
70 views

Using result of Solve in further calculations

I'm trying to make a code for Gauss-Legendre quadrature. First I find the roots of the polyinomial with roots[n_,x_]:=Solve[LegendreP[n, x] == 0] Then I need to ...
0
votes
1answer
97 views

Solving non-linear ODE analytically in mathematica [closed]

I would like to solve the following differential equations $${dx\over dt}=\beta x-{{(C-x)^2}\over 4k} ,t\in[0, T]$$ where $\beta$,$C$ and $k$ are constant values. Also, the folowing terminal ...
1
vote
0answers
264 views

Solve set of non-linear equations with least-squares-fitting - constrain results?

I'm trying to solve a set of functions to determine the material properties from a set of measurement values. (To set up this method I just want to fit my model with some already calculated data). My ...
0
votes
0answers
39 views

Why does `NSolve` not solve this system of equations on a bounded domain?

I have three functions $E1,E2,E3$ of $t1$ and $t2$. The expressions are slightly complicated and involve trigonometric functions. For clarity purpose, they are given at the end of the question. I ...
1
vote
2answers
77 views

What is the role of the pattern object in this implementation?

In the documentation I found a nice idiom which allows me to do model fitting with interpolating functions. I think this kind of thing is amazing and a credit to the utility of an environment like ...
1
vote
1answer
178 views

Solve a system of non-linear equations

Is there any way to get an exact solution a system of non-linear equations? The system I would like to solve is: $\quad \quad \begin{eqnarray} h_{1}x_{1}=c_{1}\\ h_{2}x_{2}=c_{2}\\ ...
1
vote
1answer
209 views

Eliminating parameter from parametric equation [closed]

In Mathematica Online I tried: Eliminate[{x == t + t^3, y == t - t^3, z == 1 + t^4}, {t}] ...
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: ...
1
vote
1answer
120 views

How can I solve an equation independent of a variable? [duplicate]

Let's say I have the equation $$ ax + b = 0 $$ and I'd like to find a solution that doesn't depend on $x$, so here the solution is obviously $a=b=0$. I could solve this using ...
2
votes
0answers
71 views

Problem with plotting a function that calls FindRoot [closed]

I am having difficulties in plotting the solution of FindRoot to solve 4 variables. I defined a function and then tried tp plot it, but got an error message ...
2
votes
2answers
161 views

Why does Mathematica not evaluate the difference recurrence using `RSolve` or `RecurrenceTable`?

Mathematica does not evaluate the following recurrence of two variables using RSolve and only echoes the input: ...
0
votes
2answers
149 views

Solve returns empty list / running forever

Update: I made a mistake before and corrected the sign in code below in case people wonder what some comments mean. Thank you for those who answered. But I should have said that this is only part of ...
21
votes
4answers
4k views

How to get intersection values from a parametric graph?

I have this graphed: ParametricPlot[{2.4*Cos[t] + 1.6*Cos[3 t/2], 2.4*Sin[t] - 1.6 Sin[3 t/2]}, {t, 0, 4*Pi}] It is a star, and the lines cross each other ...
1
vote
1answer
135 views
0
votes
0answers
79 views

Inverse Superposition

There are 4 tables that, I have written them as bellow: (m, n, s, and r are orthogonal basis. r.r=n.n=m.m=s.s=1, and m.n=n.r=...=0). a's, b's, c's and d's are coefficients and are numeric values. ...
2
votes
1answer
196 views
2
votes
2answers
132 views

Solve[ ] doesn't work with transcendental equations

How to solve x^5-sinx=0 in mathematica? I try in this way: Solve[x^5 - Sin[x] == 0, x],but t doesn't work.
1
vote
4answers
5k views

How can I use assumptions with FindRoot?

I have something like FindRoot[f[x], {x, a}] Now I want FindRoot to constrain the solutions to 0 < x < 1. How can I ...
7
votes
6answers
761 views

Finding where two parametric paths cross?

I have been looking through other parametric equations questions and am still unsure where to start with this one. I am given two functions: ...
0
votes
3answers
135 views

Problems with finding the real solutions of equation

I am trying to find real solutions of equation using Mathematica 8.0 but there is some problem I don't know how to explain. For example, for the input ...
0
votes
1answer
56 views

Unexpected ReplaceAll behavior with NSolve

I ran into some unexpected behavior using ReplaceAll with an equation of motion using NSolve (or Solve, also). First solving for the initial velocity v0: ...
3
votes
2answers
94 views

Switching between solutions of a fourth order polynomial

I am plotting the roots of a fourth-order polynomial as a function of some parameter $r$ with the following code: ...
0
votes
0answers
33 views

why is NSolve giving {y,x} in place of {x,y}? [duplicate]

Why is NSolve giving an inverse set of solution? Is it a problem with my computer only? ...
2
votes
1answer
163 views

How to Solve[] for rational constants?

I am trying to solve a system equations for rational, constant coefficients: ...