Hi Jon, Your best bet is to use the interval arithmetic in Sage since it will give you proveable intervals for your results. Here's an example of how to use it:
sage: CIF = ComplexIntervalField(20) sage: A = random_matrix(CIF, 3, 3) sage: b = vector([1,1,1]) sage: A \ b ([-1.0000039 .. -0.99999618], [-1.0000287 .. -0.99998092], [1.9999847 .. 2.0000191]) sage: A*_ ([0.99995422 .. 1.0000535], [0.99990844 .. 1.0000802], [0.99998092 .. 1.0000229]) --Mike On Mon, Feb 25, 2008 at 3:53 AM, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > > Hi, > > I am interested in being able to find a numerical solution of a matrix > system of linear equations Ax=b, where A is a matrix (and b is a > vector) over ComplexField(prec), of a specified precision prec. The > important thing is to know the level of accuracy of the resulting > solution vector x. (Here I can assume the vector b is known exactly, > but the matrix A is only known up to the specified precision.) I was > wondering what can be done along these lines with SAGE, and the many > numerical packages it contains? I would like to use some pre-existing > implementation if possible. > > Thanks, > > -Jon > =) > > > --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-devel URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
