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

learn more… | top users | synonyms (5)

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
575 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
88 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
176 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
7k 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
130 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
159 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
103 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 e^{-c_{1,1}}-0....
4
votes
1answer
314 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
1k 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
53 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
286 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
554 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
172 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
492 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
249 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 there ...
3
votes
1answer
76 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
274 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 ...
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
183 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}\\ h_{3}(1-x_{1}-x_{...
1
vote
1answer
213 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
186 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
126 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
72 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
164 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
154 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
138 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
200 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
786 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
136 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
95 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
167 views

How to Solve[] for rational constants?

I am trying to solve a system equations for rational, constant coefficients: ...
1
vote
2answers
277 views

Trouble using Solve and NSolve with functions involving Erf

I have the following functions: R[k_, x_, t_] := -.5*(k - x)*(1 + Erf[-(k - x)/t]) L[k_, c_, x_, t_] := .5*c*(k - x)*(1 + Erf[(k - x)/t]) I'm interested in ...
3
votes
1answer
109 views

Inaccurate answer for a system of equations

When I solve the following system of 6 equations by NSolve, it gives me an answer set; however this answer won't make the 2nd equation equal to zero (i.e. when I put the answers into ...
-1
votes
1answer
82 views

Find a minimum value which satisfies an equation

I have to find a minimum value of "a" which makes x=-a+b+c, y=-a+b+c^2 positive where a>1, 0 ...