I need to solve a matrix equation of the form A u = q for u, where u and q are N-vectors and A is a symmetric-positive-definite matrix. For now, at least, A is a 2-D matrix. I would like to use Python and numarray. Are there Python modules that I could use to obtain a solution for u? Thanks, Craig