( examples)
( values to be floating point)

3 #rows !
4 #columns !

#rows @  #columns @ ARRAY[] matrix
#rows @  ARRAY sequence        ( to be used to 'rearrange' order of equations)
#rows @  ARRAY x(i)            ( to hold solutions x0,x1 etc)

: zero_matrix #rows @ 0 DO  #columns @ 0 DO   0 J I matrix !   LOOP  LOOP ;

: load_matrix ( a0,a1,a2,a3,row...)
   0 #columns @ 1- DO DUP >R  I  matrix    ( a0,a1,a2,a3,addr...)
                      !   R>               ( a0,a1,a2,row...)
               -1 +LOOP DROP ;

: dsp_matrix #rows @ 0 DO CR  #columns @ 0 DO  J I matrix @ F. SPACE LOOP  LOOP ;

zero_matrix


3.0 FLOAT  2.0 FLOAT  4.0 FLOAT  4.0 FLOAT   0 load_matrix
1.0 FLOAT  0.0 FLOAT  5.0 FLOAT  0.0 FLOAT   1 load_matrix
1.0 FLOAT  1.0 FLOAT  6.0 FLOAT  2.0 FLOAT   2 load_matrix SP@ CR . CR dsp_matrix

EXIT

