# [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

end if

end for

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.

end for

i := i   1

end if

j := j   1

end while
Can you please guide me doing the program?
Thank you.

_________________________________________________________________
Easily edit your photos like a pro with Photo Gallery.
http://get.live.com/photogallery/overview
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.math.utexas.edu/pipermail/maxima/attachments/20080810/07c7f0fd/attachment.htm
```