在 Maxima 中计算矩阵的简化行阶梯形式#

1
rref(a):=block([p,q,k],[p,q]:matrix_size(a),a:echelon(a),
2
    k:min(p,q),
3
    for i thru min(p,q) do (if a[i,i]=0 then (k:i-1,return())),
4
    for i:k thru 2 step -1 do (for j from i-1 thru 1 step -1 do a:rowop(a,j,i,a[j,i])),
5
    a)