Mathematica Stack Exchange is a question and answer site for users of Mathematica. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I have following set of coupled transcendental equations.

$\tan^{-1} \left(\frac{\sigma+\beta}{K'-1}\right)+\tan^{-1} \left(\frac{\sigma-\beta}{K'-1}\right)=\frac{2 \sigma}{K'}\cos(\beta\tau)$
$\log\left(\frac{(K'-1)^2+(\sigma-\beta)^2}{(K'-1)^2+(\sigma+\beta)^2}\right)=\frac{4 \sigma}{K'}\sin(\beta\tau)$

Here K', $\sigma$ and $\beta$ are variables and $\tau=0.07$ is a constant. I want to obtain $\sigma$ vs K' curve (as shown in the attached figure) by eliminating $\beta$.

I tried to eliminate $\beta$ by using Eliminate function of Mathematica, but it gives:

Eliminate::ifun: Inverse functions are being used by Eliminate, so some solutions may not be found; use Reduce for complete solution information.

When a tried using Reduce, Mathematica keeps on running and does not give anything.

I also tried using FindRoot and NSolve to find values of K' and $\sigma$ that satisfy these equations for a range of $\beta$ going from $-2\pi/\tau$ to $2\pi/\tau$, but I keep on get either $\sigma\approx 0$ or very high.

Do[Print[FindRoot[{ArcTan[(sigma + beta)/(k - 1)] + 
  ArcTan[(sigma - beta)/(k - 1)] == ((2 sigma) Cos[0.07  beta])/k,
 Log[((k - 1)^2 + (sigma - beta)^2)/((k - 1)^2 + (sigma + 
     beta)^2)] == ((4 sigma) Sin[0.07  beta])/k}, {{sigma, 
 2.5}, {k, 2.0}}]], {beta, -2 Pi/0.07, 2 Pi/0.07, 0.01}]

I am relatively new to Mathematica, so please suggest me ways to handle such equations.enter image description here

share|improve this question
Please provide Mathematica code with/instead of latex equations. – Öskå Jun 3 '14 at 18:52
I don't have Mathematica at hand right now, but you may try to divide your equations by sigma (see my answer below), so that sigma=0 is no longer a solution to your system. – user10957 Jun 4 '14 at 7:41
@Jean-ClaudeArbaut Thank you very much for your answer. I used your suggestion of dividing both equations by (1/sigma). Though making surface plots for the functions and looking at their intersection at [0,0] plane didn't help in this case but it really helped for a similar set of equations. For this particular case, the surfaces do not intersect at any point other than sigma equal to zero. Anyways, I can't thank you enough for your help. – miki Jun 6 '14 at 18:22

Notice σ=0 is always a solution, so the problem is not well-posed.

It may help to have a look at the two equations separately:

ContourPlot3D[ArcTan[(\[Sigma] + \[Beta])/(k - 1)] + 
  ArcTan[(\[Sigma] - \[Beta])/(k - 1)] == ((2 \[Sigma])
  Cos[0.07 \[Beta]])/k,
  {\[Sigma], -5, 5}, {\[Beta], -5, 5}, {k, -5, 5}]

enter image description here

ContourPlot3D[Log[((k - 1)^2 + (\[Sigma] - \[Beta])^2)/
  ((k - 1)^2 + (\[Sigma] + \[Beta])^2)] ==
  4 \[Sigma] Sin[0.07 \[Beta]]/k,
  {\[Sigma], -5, 5}, {\[Beta], -5, 5}, {k, -5, 5}]

enter image description here

share|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.