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

20 REM * *

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

40 REM * Με τη μέθοδο των Gauss-Jordan *

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

60 REM * *

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

80 CLS

90 EPS=.0000005

100 DIM A(20,20),B(20,20)

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

120 NC=2*N

130 PRINT

140 FOR I=1 TO N

150 FOR J=1 TO NC

160 IF J>N THEN GOTO 190

170 READ A(I,J)

180 GOTO 230

190 IF ABS(I-J)=N THEN GOTO 220

200 A(I,J)=0

210 GOTO 230

220 A(I,J)=1

230 NEXT J

240 NEXT I

250 GOSUB 390

260 PRINT

270 PRINT"--------------"

280 PRINT "Αποτελέσματα"

290 PRINT"---------------"

300 FOR I=1 TO N

310 FOR J=N+1 TO NC

320 PRINT A(I,J),

330 NEXT J

340 PRINT

350 NEXT I

360 END

370 REM Υποπρόγραμμα που βρίσκει τον αντίστροφο

380 REM πίνακα με τη μέθοδο των Gauss-Jordan

390 FOR K=1 TO N

400 REM Μερική οδήγηση

410 W=A(K,K):LC=K

420 FOR IT=K+1 TO N

430 IF ABS(A(IT,K))<ABS(W) THEN GOTO 460

440 W=A(IT,K)

450 LC=IT

460 NEXT IT

470 REM Αλλαγή γραμμών

480 FOR LT=1 TO NC

490 PROP=A(K,LT)

500 A(K,LT)=A(LC,LT)

510 A(LC,LT)=PROP

520 NEXT LT

530 REM Κανονικοποίηση της οδηγού γραμμής

540 FOR J=1 TO NC

550 B(K,J)=A(K,J)/W

560 NEXT J

570 FOR I=1 TO N

580 IF I=K THEN GOTO 620

590 FOR J=1 TO NC

600 B(I,J)=A(I,J)-A(I,K)*B(K,J)

610 NEXT J

620 NEXT I

630 REM Αντικατάσταση του Β πίνακα με τον Α

640 FOR I=1 TO N

650 FOR J=1 TO NC

660 A(I,J)=B(I,J)

670 NEXT J

680 NEXT I

690 NEXT K

700 RETURN

710 DATA 1,3,2

720 DATA 3,8,5

730 DATA 2,5,4