On Jun 29, 3:31 pm, Colin <[email protected]> wrote: > Dear Sympy list, > I'm trying to evaluate a numerical dispersion relation for a finite > element method we are working on, which requires construction of a 4x4 > Hermitian matrix A which is a function of l,k,w, and evaluation of > it's determinant. This gives a quartic polynomial which can be solved > for w for each l,k. I want to do this symbolically because then I can > trace the various root branches
Just a note of caution: sympy will give you a general solution to the quartic equation but it won't be valid for all values of l and k. You can see the quartic root solver to see the different solutions that are returned based on the various "determinant-like" quantities that are computed. If, when substituting in l and k values you get a w that doesn't satisfy the original equation you will know that you have to recalculate the roots with those specific values. -- You received this message because you are subscribed to the Google Groups "sympy" group. To post to this group, send email to [email protected]. To unsubscribe from this group, send email to [email protected]. For more options, visit this group at http://groups.google.com/group/sympy?hl=en.
