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

learn more… | top users | synonyms (5)

2
votes
1answer
235 views

NDSolve fails for certain choices of parameters and solve range

I'm trying to solve a pair of coupled ODEs with NDSolve. I know roughly what the solution should look like (both should give periodic functions, pi/2 out of phase, the amplitude of which damp towards ...
1
vote
1answer
168 views

Fitting composed model to experimental data

I have mathematical model based on equation: $y(x) = 1 - \sum_{b=1}^{\inf} 2c_b \sin(c_b x/L) \exp(ax) $ (1), where $c_b$ are the positive roots of equation: $c_b \cot c_b + aL = 0$ (2). $L$ is a ...
7
votes
1answer
323 views

Why does Mathematica not give simple solutions to simple system of equations with powers

I trt to get Mathematica getting me the simple solutions for x and y from this system of two power equations: $$\text{Eq. 1: } z=x^{m}y^{n}$$ $$\text{Eq. 2: }z=\tau x^{m-1}y^{1-n},$$ The simple ...
2
votes
1answer
213 views

Numerically solving systems of (first order) linear delay differential equations

I was wondering, how to solve such a system with mathematica? I found the NDsolve function, but in the reference is not mentioned if it is possible to solve a system of equations.
2
votes
1answer
271 views

Find solutions of equation involving Bessel functions

I'm new in Mathematica and I'm trying to find the solutions of this equation involving Bessel functions $$\eta \frac{ J_{n+1}(\eta a)}{J_n(\eta a)}+\chi \frac{ I_{n+1}(\chi a)}{I_n(\chi a)}=0$$ ...
3
votes
2answers
255 views

What would be the inverse function for the following condition?

I tried the problem but i couldn't take above condition in question in terms of f(x)
2
votes
2answers
224 views

Extract data from NDSolveValue result

Okay this is a very newbie question. I have coded a sample problem that solves the heat diffusion equation on an annulus. Here's what I have: ...
2
votes
0answers
183 views

Solving a system of difference equations in an exponential steady-state growth setting

I have a system of difference equations in an exponential steady-state balanced growth setting; all variables grow at the same steady-state rate, say, $g$. I use the notations as follows: Variable $X$ ...
5
votes
2answers
289 views

How to get $x$ in terms of $y$ and $w$ from $w = a p x^{p-1} \left(a x^p+b y^p\right)^{\frac{1}{p}-1}$

I'm having trouble to solve the following equation for $x$ both manually and in Mathematica: $\quad \quad w = a p x^{p-1} \left(a x^p+b y^p\right)^{\frac{1}{p}-1}$ In Mathematica notation: ...
1
vote
1answer
83 views

m equations and n variables. Reducing using Mathematica

I have six equations and thirteen variables. I would like to reduce it to seven free parameters. I used Reduce to solve it but it is taking really long (one full day and no output yet) and I think ...
3
votes
1answer
514 views

Mathematica takes too long to solve this

I have the following system of equations: ...
2
votes
1answer
229 views

How to solve three equations with four unknowns iteratively?

I have three equations of the form $\quad \eta = f(N1L, N2L)$ $\quad N1L = g(\eta, R)$ $\quad N2L = h(\eta, R)$ I want to eliminate $\eta$, so that I can get $N1L$ and $N2L$ for a range of $R$. I ...
3
votes
3answers
186 views

Solving determinant of a Kronecker product of matrices gives a numerical error - why?

I am doing the following steps (code at the end of the post): I start with a 2x2 matrix (smatrix), which is a function of a single variable (u2). I want to set the determinant of this matrix ...
8
votes
3answers
348 views

Is there a relatively quick way to solve my integer-valued equation?

This is a very simple problem. It can be solved by nested Forloops that exhaustively test every possible case, throwing out all non-solutions. Find all the ...
4
votes
1answer
709 views

How to find the integer solution number of a linear system with inequalities?

I have the following system with mixed equality and inequalities defined as: ...
2
votes
0answers
199 views

Error when using EquationTrekker

I just got to know EquationTrakker because I need to do phase graphs for my classes. When I load it ...
0
votes
0answers
107 views

What methods does Solve use?

Could you tell me which method that Solve applies to solve linear and polynomials equations?
5
votes
1answer
170 views

NDSolve solution violates initial conditions

I have the following code: ...
9
votes
2answers
370 views

How can I only draw outer common tangent to two circles?

I have code worked for drawing common tangent lines to two circles(http://en.wikipedia.org/wiki/Tangent_lines_to_circles), now I want to delete the inner common tangent lines, how? ...
7
votes
1answer
324 views

How to find the positive region of this function?

I'm new-ish to Mathematica and I'm using it to learn the basics of derivative trading (for fun not for trading!). I've defined a simple function to calculate the pay off of a trading on a call. I have ...
3
votes
0answers
86 views

How to solve for exponents in an infinite product?

Say I have a power series, $f(x,y) \in 1 + x\mathbb{Z}[[x,y]]$ Then it admits a unique expansion $$f(x,y) = \prod_{i=1}^\infty \prod_{j=0}^\infty (1 - x^i y^j)^{n_{i,j}}$$ for some integers, $i, ...
7
votes
4answers
789 views

Solving $L=\frac{3}{2} \sqrt{4 \pi ^2 A^2+W^2}-\frac{\sqrt{5 W \sqrt{4 \pi ^2 A^2+W^2}+6 \pi ^2 A^2+3 W^2}}{\sqrt{2}}+\frac{3 W}{2}$ for $W$

When I solve the aforementioned equation for $W$ or $A$ on Mathematica I get a long and ugly equation in return, namely one of the solutions for $W$ is: (attempt to read at your own health) ...
0
votes
2answers
406 views

Find Root Iteration

I have an iterative sum from $k=0$ to $k=n$ where the resulting sum is a polynomial of degree $n$. I want to find the numerical root of this polynomial using FindRoot, starting from $x_0$ where the ...
4
votes
3answers
201 views

Regarding how to use Reduce

Say for example I am doing the following in Mathematica ...
1
vote
2answers
222 views

Need help solving a system of two 1st order nonlinear differential equations

The original system of equations reads: $\begin{cases} f'(r) + f(r) \left(a(r) - \frac{1}{r}\right) = 0,\\ f^2(r) + a'(r) + \frac{a(r)}{r} - 1 = 0\,, \end{cases}$ with boundary conditions $f(0) = ...
4
votes
2answers
243 views

Can't solve equation having complex coefficients

Why can't Mathematica solve $\quad\quad\frac{1-i}{\sqrt{2}}=e^{i \alpha } \tan \left(\frac{\beta }{2}\right)$ with the restrictions $\alpha \in [0, 2 \, \pi)$ and $\beta \in [0,\pi]$: ...
9
votes
4answers
768 views

How many times do the two given curves intersect at the origin?

I'm reading/studying on my own a book about algebraic curves. I'm trying to find a way for Mathematica help me with the following (the book, ISBN: 038731802X, allows me to use any CAS / program. ...
3
votes
1answer
305 views

Issue in ParallelTable after evaluating another function using NDSolve and FindRoot. What is wrong with this inverse?

I am trying to find the inverse of a function which is defined through NDSolve and NIntegrate. The question is pretty similar ...
3
votes
1answer
111 views

Already solved DE, now I need to rearrange and plot

I have solved for $z(w,x,y)$ in a differential equation: $$ 3\frac{\partial z}{\partial y} = 2(z-1) + (1-wy^2 )x $$ And I obtained the general solution: $z = f(w,x,y)$ Now putting in $x=0$, we have ...
0
votes
2answers
92 views

Solve simple differential equation - Error?

I get the error "The function Te appears with no arguments." when running this code. I'm not sure why. Is it because of the $Abs[\Gamma]$ ? ...
1
vote
0answers
212 views

Best way to determine polynomial coefficients in series expansion

I would like to solve the equation $$h'(\boldsymbol{x}_1)\left[B_1\boldsymbol{x}_1+g_1(\boldsymbol{x_1},h(\boldsymbol{x}_1))\right]=B_2h(\boldsymbol{x}_1)+g_2(\boldsymbol{x}_1,h(\boldsymbol{x}_1))$$ ...
0
votes
1answer
166 views

Solving a Trigonometric parameter

I have a similar case as given before : Solving a Trigonometric function with one parameter Let: m = 1; Dx = 100; ν = 0.25; b = 1; and, ...
2
votes
0answers
149 views

Quartic equation [closed]

After using Solve I am choosing the physical root of my equation and dealing further with this root. However, when I am changing parameters using Manipulate, I clearly see not smooth behaviour. Is ...
0
votes
1answer
876 views

Using Solve but getting no output

I am trying to solve a model using a constant-elasticity-of-substitution production function, but when I try to solve the first order conditions, I just keep getting back what I typed. ...
1
vote
1answer
229 views

Solving equations with QuantityVariable

I'm trying to understand how to easily solve equations that use QuantityVariable. Why does this fail? ...
0
votes
0answers
59 views

Resources/ best practices for exploratory economic theory work

This is a little general, so I will try to make it as specific as possible. I used Mathematica a fair amount to derive and check analytical solutions to problems in theoretical economics. Lots of ...
1
vote
2answers
247 views

How do I plot the differential of this function?

Suppose I have this function: $$z(x,y) = \left| \frac{ \frac{1}{3x +iy} -2x}{iy + \frac{1}{x}} \right| $$ I want the contour plot of $\frac{\partial z}{\partial x}$ with axes $(x,y)$. Tried this ...
1
vote
1answer
100 views

Whats the most efficient way to solve an equation numericaly (because it has no analitical solution)

So, I think the problem that Im having is simple but, still, Im not sure on how to do it. I have an equation with no analitical solution: $a_1 \sin \left(2 \theta \right)+a_2 \sin \left(2 ...
1
vote
4answers
435 views

How to extract only real solutions from the result of Solve

For Example, I have the following polynomial equations: $$ \left\{ \begin{array}{ll} \dfrac{2025 (208 x+5 (y-3446))}{52 (y+90)^2}+\dfrac{300 (8 x-21 y+1920) (y-80)}{7 (x+30)^3}+\frac{300 (80 x-21 ...
1
vote
3answers
486 views

Finding maximum value with position from Table of values

I have this table of values, how do I find the maximum value? i.e. find $z_{max}(x,y)$ ...
0
votes
1answer
94 views

Solving minimum length confidence interval problem

I'm trying to find the minimum length of a confidence interval for the variance from a sample iid with normal distribution. I tried ...
15
votes
3answers
525 views

Solving the magic backyard puzzle with MMA

I was trying to solve this puzzle using mma It looks like the generalization of the magic square. Anyway here my attempt:- There are four columns $(x,y,z,w)$ and nine rowa which are denoted by the ...
0
votes
0answers
41 views

How can Mathematica solve 7 equations and want to find 7 variables [duplicate]

I have 7 equations as below * e1 := (1 - a) p Subscript[A, 0] Subscript[s, K ] ! *SubsuperscriptBox[(s), (L), (-a)] = (1 - b) Subscript[B, 0] (1 - Subscript[s, K]) (1 - Subscript[s, ...
3
votes
1answer
911 views

How do I solve for one variable in terms of other variables in Mathematica?

I'm fairly new to Mathematica, and I'm trying to solve a complicated equation with multiple variables. I've tried the following... ...
0
votes
2answers
131 views

Solving an Equation (Recursively)

I have the following Values: $a[0]=1;a[1]=0;a[2]=0;a[3]=0;a[4]=\frac{g_2}{20};a[5]=0$. I need to compute the value of $b[n]$ in the following equation $$a[n]=\sum _{m=0}^{n} (m-2) \;\;a[m] ...
1
vote
1answer
125 views

How do I solve for the real and complex parts of an equation simultaneously?

If I were to have an equation, say something similar to... (1-Sqrt[x - I y])/(1+Sqrt[x - I y]) = A + I B Where I = Sqrt[-1], is there a way for Mathematica to ...
0
votes
1answer
62 views

Why can NSolve only solve this system of equations after trivial substitution?

The first system immediately yields {}, but the second, equivalent, system yields the solution {{k -> 4.96048, tl1 -> 1.64842, tl2 -> 2.07466}} : ...
0
votes
0answers
34 views

What is that Root[…] inside the answer represent? [duplicate]

So I am encountering this situation, in which mathematica doesn't give me the answer, but in a form of "Root". How do I interpret this ? what is the number after the "&" represent? How should I ...
2
votes
0answers
116 views

How to handle infeasible points in FindRoot?

I am calling FindRoot[f[x,y],{{x,xInit,xMin,xMax},{y,yInit,yMin,yMax}}] where for some points {x,y}, ...
0
votes
1answer
170 views

Simplifying a symbolic product of rational powers

I'm trying to symbolically evaluate a lot of equations that can consist of up to 10 parameters. The form of these expressions is a product/quotient of rational powers of the mentioned parameters, ...