10 REM ********** ΠΡΟΓΡΑΜΜΑ Π304 **********

20 REM * *

30 REM * Αντιστροφή πίνακα *

40 REM * μέθοδος απαλοιφής του Gauss *

50 REM * μερική οδήγηση *

60 REM * *

70 REM ***************************************

80 CLS

90 EPS=.0000001

100 DIM A(20,21)

110 INPUT "Δώσε τη διάσταση του πίνακα ";N

120 NC=N+N

130 FOR I=1 TO N

140 FOR J=1 TO N

150 READ A(I,J)

160 NEXT J

170 NEXT I

180 FOR I=1 TO N

190 FOR J=N+1 TO NC

200 IF J-I=N THEN GOTO 230

210 A(I,J)=0

220 GOTO 240

230 A(I,J)=1

240 NEXT J

250 NEXT I

260 GOSUB 400

270 CLS:LOCATE 5,5

280 PRINT "Ο αντίστροφος πίνακας είναι"

290 PRINT

300 FOR M=1 TO N

310 FOR J=N+1 TO NC

320 PRINT A(M,J),

330 NEXT J

340 PRINT

350 NEXT M

360 END

370 PRINT "Η ορίζουσα του Α ισούται με μηδέν "

380 PRINT "δεν υπάρχει αντίστροφος"

390 END

400 REM Απαλοιφή Gauss με μερική οδήγηση

410 FOR JR=1 TO N-1

420 W=0

430 FOR K=JR TO N

440 IF ABS(A(K,JR))>W THEN IC=K:W=A(K,JR)

450 NEXT K

460 IF IC=JR THEN GOTO 520

470 FOR IT=JR TO NC

480 F=A(IC,IT)

490 A(IC,IT)=A(JR,IT)

500 A(JR,IT)=F

510 NEXT IT

520 FOR I=JR TO N-1

530 IF ABS(A(JR,JR))<EPS THEN GOTO 370

540 P=A(I+1,JR)/A(JR,JR)

550 FOR J=JR+1 TO NC

560 A(I+1,J)=A(I+1,J)-P*A(JR,J)

570 NEXT J

580 NEXT I

590 NEXT JR

600 REM Πίσω αντικατάσταση

610 FOR II=N+1 TO NC

620 FOR M=N TO 1 STEP -1

630 IF M=N THEN GOTO 670

640 FOR L=N TO M+1 STEP -1

650 A(M,II)=A(M,II)-A(M,L)*A(L,II)

660 NEXT L

670 IF ABS(A(M,M))<E THEN GOTO 370

680 A(M,II)=A(M,II)/A(M,M)

690 NEXT M

700 NEXT II

710 RETURN

720 DATA 1,3,2

730 DATA 3,8,5

740 DATA 2,5,4