Listing of file='OP1D070A' on disk='vmedia/701-2046C.wvd.zip'
# Sector 457, program filename = 'OP1D070A' 0010 REM . OP1D070A, 00-00(06/28/75), 12002 0100 PRINT "TESTING MAT INV... "; : DIM A(5,5),B(10,10),C(10,10) : A(1,2)=2 : A(2,5)=8 : A(3,4)=5 : A(4,3)=4 : A(5,1)=1 : MAT B=INV(A),D : IF D<>-320THEN 840 : MAT C=A*B : MAT A=IDN : MAT C=C-A : FOR I=1TO 5 : FOR J=1TO 5 : IF C(I,J)<>0THEN 840 : NEXT J 0260 NEXT I : I=RND(0) : FOR M=1TO 10 : MAT REDIM A(M,1),B(M,M) : FOR I=1TO M : FOR J=1TO M : B(I,J)=10*RND(1)-5 : NEXT J : NEXT I : MAT C=INV(B) : FOR S=1TO M : FOR T=STO M : IF B(T,S)<>0THEN 410 : NEXT T : GOTO 840 0410 GOSUB 680 : GOSUB 750 : B(S,S)=1/B(S,S) : FOR T=1TO M : IF T=STHEN 490 : B=-B(T,S) : B(T,S)=0 : GOSUB 800 0490 NEXT T : NEXT S : FOR S=MTO 1STEP -1 : IF A(S,1)=STHEN 580 : FOR J=1TO M : B=B(J,S) : B(J,S)=B(J,A(S,1)) : B(J,A(S,1))=B : NEXT J 0580 NEXT S : FOR I=1TO M : FOR J=1TO M : IF C(I,J)<>B(I,J)THEN 840 : NEXT J : NEXT I : NEXT M : IF E(7)<>0THEN 850 : PRINT "OK" : GOTO 870 0680 FOR J=1TO M : B=B(S,J) : B(S,J)=B(T,J) : B(T,J)=B : NEXT J : A(S,1)=T : RETURN 0750 FOR J=1TO M : IF J=STHEN 780 : B(S,J)=B(S,J)/B(S,S) 0780 NEXT J : RETURN 0800 FOR J=1TO M : B(T,J)=B(T,J)+B*B(S,J) : NEXT J : RETURN 0840 E(7)=E(7)+1 0850 PRINT TAB(40);E(7);"ERRORS" : E9$(P9)= "NG" 0870 LOAD DC T#6,"OP1D080A" 0880 DEFFN'15 : COM CLEAR : LOAD DC T#6,"START"