LOZAN@TAUNOS.BITNET (01/08/90)
Hello
I found a bug in my "det" program for symbolic matrices and also a typing
error.I also added a zerochecking so that det multiplied by zero will not be
computed .
Here is the corrected version.
det [D81F] (Lowercase are important !!)
<< DUP SIZE -> A B
<<
IF B 1 == ; det of a 1x1 matrix
# THEN A 1 GET 1 GET
ELSE
IF B 2 ==
THEN A DET2 ; call basic case of 2x2 determinant
ELSE 0
1 B FOR I
A 1 GET I GET
# -> E << IF E 0 SAME THEN ELSE
# A 1 I MINOR det E ; recursive call
* -1 1 I + ^ * +
# END >>
NEXT
END
END
>>
>>
Lines with a # are line with changes
I'd like to get feedback if someone find these programs useful.
Thats all for today.
Eliel Louzoun <LOZAN@TAUNOS>