Reduced row echelon form: Difference between revisions

m
Line 1,857:
=={{header|Maxima}}==
<lang maxima>rref(a):=block([p,q,k],[p,q]:matrix_size(a),a:echelon(a),
k:min(p,q),
for i:p thru 2 step -1 do (
for i thru min(p,q) do (if nota[i,i]=0 then integerp(k:for j thru q doi-1,return())),
for i:k thru 2 step -1 do (for ifj from a[i,j]#0-1 thenthru return(j))1 thenstep return-1 do a:rowop(a,j,i,a[j,i])),
for j from i-1 thru 1 step -1 do a:rowop(a,j,i,a[j,k])),
a)$
 
ma: matrix([1,212,-127,-4]36,[244,3,-1,-1159],[-2,0,-3,22]);
[26,41,-54,24,23],
rref(m);
[33,70,59,15,-68],
matrix([1,0,0,-8],[0,1,0,1],[0,0,1,-2])</lang>
[43,16,29,-52,-61],
[-43,20,71,88,11])$
 
rref(ma);
matrix([1,0,0,-80,1/2],[0,1,0,0,-1],[0,0,1,0,-1/2],[0,0,0,1,1],[0,0,0,0,0])</lang>
 
=={{header|Objeck}}==
Anonymous user