Pivotisation of a matrix - I try to change an existing file that works with another scheme
I have a Mathcad file that do the pivotisation of a matrix but I dont like it because is too indirect and obscure.
I want to implement another scheme that is more explicit and easily understanding but does not works by now. I attached my working file that also contains the scheme that i try to implement.
In your Piv(Aug) function, you let i run through all rows (from 1 to n), then you start k at i+1. That is bound to fail when i=n, because then k=n+1, which is out of range.
Potential remedy: in the first FOR statement, change:
"for i e 1..n" => "for i e 1..n-1"
It is also useful that function, but I am interested on Mathcad programming technics so that to use minimally the direct prefabricated functions.