Listing of file='20' on disk='vmedia/glbr_22.wvd.zip'
# Sector 241, program filename = '20' 0001 DIM A(5,5),B(5,5) 0005 PRINT "ENTER N" : INPUT N 0010 PRINT : PRINT "ENTER MATRIX" 0015 FOR I=1 TO N : PRINT : PRINT "COL";I 0020 FOR J=1 TO N : INPUT A(J,I),B(J,I) : NEXT J 0025 NEXT I 0030 D1,I1=1 : D2=0 0035 I3=I1 : S=ABS(A(I1,I1))+ABS(B(I1,I1)) 0040 FOR I=I1 TO N : T=ABS(A(I,I1))+ABS(B(I,I1)) 0045 IF S>=T THEN 50 : I3=I : S=T 0050 NEXT I : IF I3=I1 THEN 75 0055 FOR J=1 TO N 0060 S=-A(I1,J) : A(I1,J)=A(I3,J) : A(I3,J)=S 0065 S1=-B(I1,J) : B(I1,J)=B(I3,J) : B(I3,J)=S1 0070 NEXT J 0075 I3=I1+1 0080 FOR I=I3 TO N : S1=A(I1,I1)^2+B(I1,I1)^2 0085 S=(A(I,I1)*A(I1,I1)+B(I,I1)*B(I1,I1))/S1 0090 B(I,I1)=(A(I1,I1)*B(I,I1)-A(I,I1)*B(I1,I1))/S1 0095 A(I,I1)=S : NEXT I 0100 J2=I1-1 : IF J2=0 THEN 150 0110 FOR J=I3 TO N : FOR I=1 TO J2 0120 A(I1,J)=A(I1,J)-A(I1,I)*A(I,J)+B(I1,I)*B(I,J) 0130 B(I1,J)=B(I1,J)-B(I1,I)*A(I,J)-A(I1,I)*B(I,J) 0140 NEXT I : NEXT J 0150 J2=I1 : I1=I1+1 0160 FOR I=I1 TO N : FOR J=1 TO J2 0170 A(I,I1)=A(I,I1)-A(I,J)*A(J,I1)+B(I,J)*B(J,I1) 0180 B(I,I1)=B(I,I1)-B(I,J)*A(J,I1)-A(I,J)*B(J,I1) 0190 NEXT J : NEXT I 0200 IF I1<>N THEN 35 : I3=1 : J2=INT(N/2) 0210 IF N=2*J2 THEN 230 : I3=0 0220 D1=A(N,N) : D2=B(N,N) 0230 FOR I=1 TO J2 : J=N-I+I3 0240 S=A(I,I)*A(J,J)-B(I,I)*B(J,J) 0250 S1=A(I,I)*B(J,J)+A(J,J)*B(I,I) 0260 T=D1*S-D2*S1 : D2=D2*S+D1*S1 : D1=T 0270 NEXT I : PRINT 0280 PRINT "DETERMINANT:";D1;"+I*";D2 : PRINT 0290 END