Advanced Search

Question: LinearSolve gives wrong answer?

Question:LinearSolve gives wrong answer?

susan_likes_math 73 Maple
February 20 2006
0
Hello, all! I have a linear system of real numbers. I copied the example in the tutorial for using LinearSolve and converted my matrix to a float: with(LinearAlgebra); my_solve := proc(A::Matrix) local sz, local_A, B, sol; sz := Dimension(A); local_A := Matrix(A, datatype=float); B := Vector(1...sz[1], 1); sol := LinearSolve(local_A, B); end proc; And when it solved, it got the wrong answer!! It would produce a solution that simply didn't work. When I removed the float conversion step, and just used my original matrix, it worked perfectly. I would love to understand this better... is it because of rounding? Comments would be lovely! Take care, Susan
Please Wait...