MathGroup Archive 2010

[Date Index] [Thread Index] [Author Index]

Search the Archive

root finding routine

  • To: mathgroup at
  • Subject: [mg113600] root finding routine
  • From: nt <sagittarius5962 at>
  • Date: Thu, 4 Nov 2010 03:59:29 -0500 (EST)


I'm trying to find the roots of the following fourth order equation in
w and k:

eq = w^4 -a w^2 k +b  k^4 /. k -> (kx^2 + ky^2)^0.5;

eq1 = eq /. {w -> w1, kx -> k1x, ky -> 0};
eq2 = eq /. {w -> w2, kx -> k2x, ky -> k2y};
eq3 = eq /. {w -> w3, kx -> k3x, ky -> k3y};

given w1 and teta, and that  the roots for w and k should satisfy
these linear conditions:

sol = FindRoot[{eq1 == 0, eq2 == 0, eq3 == 0, k3x == k1x - k2x,
   w3 == w1 - w2, k2y == k2x*Tan[teta], k3y == -k2y}, {w1,
  10}*), {k2x, .81}, {k2y, -.5}, {k3x, 1.83}, {k3y, -10}, {w3,w1/2},

The above system has four sets of solutions.The findroot routine gives
one of the acceptable sets of solutions but could not get the others.
I played alot with the initial values but it seems that the routine
cannot converge to the other roots even if I put the actual other root
values as the initial guess. I would appreciate it if anyone could
provide helpful suggestions on the problem.


  • Prev by Date: Re: EventHandler vs ControllerState
  • Next by Date: Manipulate with Matrices
  • Previous by thread: Re: Embed extra info?
  • Next by thread: Manipulate with Matrices