Path: utzoo!utgpu!news-server.csri.toronto.edu!cs.utexas.edu!sdd.hp.com!spool.mu.edu!uunet!orca!javelin.es.com!pashdown From: pashdown@javelin.es.com (Pete Ashdown) Newsgroups: comp.sys.handhelds Subject: Re: How do I row reduce matrices on a 48? Message-ID: <1991Feb14.173257.2538@javelin.es.com> Date: 14 Feb 91 17:32:57 GMT References: <1991Feb14.133413.12634@mathrt0.math.chalmers.se> Reply-To: pashdown@javelin.sim.es.com Organization: Evans & Sutherland Computer Corp., Salt Lake City, Utah Lines: 42 d0mo@dtek.chalmers.se (Mats Olsson) writes: >To solve a un-quadratic system of linear equations, one have to >reduce the matrix into echelon form. This is easily done by my >program ->RED which I will post here very soon.. I got this off of the mail server. Its amazingly small and efficient in comparison to other row-reducing programs. It works great! --------------------------------------------------------------------------- %%HP: T(3)A(R)F(.); @ RWRD - Matrix Row Reduction ( matrix - reduced matrix ) @ Robert T. Wilson - University of Iowa \<< DUP OBJ\-> OBJ\-> DROP \-> r c \<< r r 1 \-> m n s \<< 1 r START c \->ARRY DEPTH ROLLD NEXT 1 r START DEPTH ROLL NEXT 1 r START 0 's' STO 1 c FOR z IF DUP z GET 0 \=/ THEN z 's' STO c 'z' STO END NEXT IF s 0 > THEN DUP s GET IF 1 \=/ THEN DUP s GET / END 1 r 1 - START r ROLL IF DUP s GET 0 \=/ THEN DUP2 s GET * - 1 'n' STO END SWAP NEXT END r ROLLD NEXT r 'n' STO 1 c FOR z 1 r FOR y IF DUP z GET 1 == THEN n ROLLD 'n' DECR IF 1 == THEN r 'y' STO c 'z' STO END ELSE n ROLLD END NEXT NEXT 1 r START DEPTH ROLLD NEXT 1 r START DEPTH ROLL OBJ\-> DROP NEXT r c 2 \->LIST \->ARRY \>> \>> \>> --------------------------------------------------------------------------- >Mats Olsson >Chalmers University of technology >Sweden -- "Hi, we're 'Slaughter'. We'd just like to say how much we love our troops." Pete Ashdown pashdown@javelin.sim.es.com ...uunet!javelin.sim.es.com!pashdown