▁A ( i , ▁j ) ▁- ▁L ( i , ▁ 1 : j ▁- ▁ 1 ) ▁* ▁U ( 1 : j ▁- ▁ 1 , ▁j ); ▁▁▁▁▁▁▁▁▁▁▁▁ ▁end ▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁ ▁▁▁▁▁▁▁▁▁ ▁for ▁j ▁= ▁i ▁+ ▁ 1 : R ▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁ ▁U ( i , ▁j ) ▁= ▁( A ( i , ▁j ) ▁- ▁L ( i , ▁ 1 : i ▁- ▁ 1 ) ▁* ▁U ( 1 : i ▁- ▁ 1 , ▁j )) ▁/ ▁L ( i , ▁i ); ▁▁▁▁▁▁▁▁▁▁▁▁ ▁end ▁▁▁▁▁▁▁▁ ▁end ▁▁▁ ▁end ▁ ▁References ▁▁ ▁Im plementation ▁using ▁functions ▁ ▁In ▁Mat lab ▁ ▁Category : Matrix ▁decom pos itions <0x0A> </s> ▁A ▁total ▁lun ar ▁eclipse ▁took ▁place ▁on ▁May ▁ 4 , ▁ 1 9 8 5 , ▁the ▁first ▁of ▁two ▁total ▁lun ar ▁e cli ps es ▁in ▁ 1 9 8 5 , ▁the ▁second ▁being ▁on ▁October ▁ 2 8 , ▁ 1 9 8 5 . ▁ ▁This ▁lun ar ▁eclipse ▁is ▁first ▁of ▁a ▁t et rad , ▁four ▁total ▁lun ar ▁e cli ps es ▁in ▁series . ▁The ▁last ▁series ▁was ▁in ▁ 1 9 6 7 ▁and ▁ 1 9 6 8 , ▁starting ▁with ▁an ▁April ▁ 1 9 6 7 ▁lun ar ▁eclipse , ▁while ▁the ▁next ▁was ▁in ▁ 2 0 0 3 ▁and ▁ 2 0 0 4 , ▁starting ▁with ▁a ▁May ▁ 2 0 0 3 ▁lun ar ▁eclipse . ▁