Hi all, I have two questions about scipy.
1) When I was trying to solve a single variable equations using scipy, I found two methods: scipy.optimize.fsolve, which is designated to find the roots of a polynomial, and scipy.optimize.newton, which is used for Scalar function root finding according to the help(). I have tried both, and it seemed that both worked well, and fsolve ran faster. My questions is, which is the right choose ? 2) I have to solve a linear equation, with the constraint that all variables should be positive. Currently I can solve this problem by manually adjusting the solution in each iteration after get the solution bu using scipy.linalg.solve(). Is there a smart way ? Thanks in advance. Xiao Jianfeng -- http://mail.python.org/mailman/listinfo/python-list