[Maxima] Gaussian Elimination [Help]
Angelique Sta maria
leque_sm at hotmail.com
Sun Aug 10 00:21:10 CDT 2008
How do you program Gaussian Elimination in Maxima? I have researched a pseudocode and here it is:
i : 1;
j : 1;
while (i = m and j = n) do
Find pivot in column j, starting in row i:
maxi := i
for k := i 1 to m do
if abs(A[k,j]) > abs(A[maxi,j]) then
maxi := k
if A[maxi,j] ? 0 then
swap rows i and maxi, but do not change the value of i
Now A[i,j] will contain the old value of A[maxi,j].
divide each entry in row i by A[i,j]
Now A[i,j] will have the value 1.
for u := i 1 to m do
subtract A[u,j] * row i from row u
Now A[u,j] will be 0, since A[u,j] - A[i,j] * A[u,j] = A[u,j] - 1 * A[u,j] = 0.
i := i 1
j := j 1
Can you please guide me doing the program?
Easily edit your photos like a pro with Photo Gallery.
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Maxima