48/50 Linear System Solver

[Han]

(01-16-2015 09:52 AM)Werner Wrote:  

(01-15-2015 05:35 PM)Han Wrote:  Edit: So I did the LU decomposition of the matrix, and got a different solution: [[84.9666...],[114.3666...],[-300]].

How did you arrive at this?
And the '/' operator does LU-decomposition.

Werner

My stack was:

2: [[-166.2],[157.2],[-3]]
1: [[12.25,1.25,4.5],[1.25,12.25,4.5],[4.5,4.5,3]]

I typed in the following in the command line; remove the comments after the @'s

Code:

LU @ b L U P
 4 ROLL * @ L U P*b
ROT INV SWAP * @ solve Ly=P * b by y=L^-1 * P *b
SWAP INV SWAP * @ solve Ux = y by x = U^-1 y