Listing of file='OP1D070A' on disk='vmedia/wcs_diags.wvd.zip'
# Sector 450, program filename = 'OP1D070A' 0010 REM . OP1D070A, 00-00 (06/28/75), 12002 : 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 440 0140 MAT C=A*B : MAT A=IDN : MAT C=C-A : FOR I=1TO 5 : FOR J=1TO 5 : IF C(I,J)<>0THEN 440 : NEXT J : NEXT I : I=RND(0) 0180 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 250 : NEXT T : GOTO 440 0250 GOSUB 370 : GOSUB 390 : B(S,S)=1/B(S,S) : FOR T=1TO M : IF T=STHEN 280 : B=-B(T,S) : B(T,S)=0 : GOSUB 420 0280 NEXT T : NEXT S : FOR S=MTO 1STEP -1 : IF A(S,1)=STHEN 320 0300 FOR J=1TO M : B=B(J,S) : B(J,S)=B(J,A(S,1)) : B(J,A(S,1))=B : NEXT J 0320 NEXT S : FOR I=1TO M : FOR J=1TO M : IF C(I,J)<>B(I,J)THEN 440 : NEXT J : NEXT I : NEXT M : PRINT "OK" : GOTO 460 0370 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 0390 FOR J=1TO M : IF J=STHEN 410 : B(S,J)=B(S,J)/B(S,S) 0410 NEXT J : RETURN 0420 FOR J=1TO M : B(T,J)=B(T,J)+B*B(S,J) : NEXT J : RETURN 0440 E(2)=E(2)+1 0450 PRINT TAB(40);E(2);"ERRORS" 0460 LOAD DC T"OP1D080A" 0470 DEFFN'15 : LOAD DC T"START"50