home *** CD-ROM | disk | FTP | other *** search
/ Usenet 1994 January / usenetsourcesnewsgroupsinfomagicjanuary1994.iso / sources / unix / volume10 / ifp / part01 / fproot / math / linear / LU < prev    next >
Encoding:
Text File  |  1987-07-05  |  415 b   |  27 lines
  1. (* 
  2.  * LU
  3.  *
  4.  * [L,U] decomposition of matrix written as single function.
  5.  * This function is functionally identical to [L,U]
  6.  *)
  7.  
  8. DEF LU AS
  9.    IF Singleton THEN [#<<1.0>>,id]        (* definition of L *)
  10.    ELSE 
  11.       [
  12.          Li1,
  13.          Aik | [EACH #0 END,LU],
  14.      U1k
  15.       ] | 
  16.       [
  17.      [
  18.         1,
  19.         2 | [1,2|1] | apndl
  20.      ] | ApndlCol,
  21.      [
  22.         3,
  23.         2 | [1,2|2] | ApndlCol
  24.      ] | apndl
  25.       ]
  26.    END;
  27.