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

20 REM * *

30 REM * Γραμμικά Συστήματα *

40 REM * Μέθοδος Gauss-Jordan *

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

60 REM * *

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

80 CLS

90 DIM A(20,21),B(20,21)

100 EPS=.0000005

110 INPUT "Δώσε τον αριθμό των εξισώσεων";N

120 NC=N+1

130 PRINT

140 FOR I=1 TO N

150 FOR J=1 TO N

160 READ A(I,J)

170 NEXT J

180 READ A(I,NC)

190 NEXT I

200 GOSUB 300

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

220 FOR I=1 TO N

230 PRINT "X(";I;")=";A(I,NC)

240 NEXT I

250 PRINT "---------------"

260 PRINT

270 END

280 REM Το υποπρόγραμμα αυτό βρίσκει τη λύση

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

300 FOR K=1 TO N

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

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

330 FOR IT=K+1 TO N

340 IF ABS(A(IT,K))<ABS(W) THEN GOTO 370

350 W=A(IT,K)

360 LC=IT

370 NEXT IT

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

390 FOR LT=1 TO NC

400 PROP=A(K,LT)

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

420 A(LC,LT)=PROP

430 NEXT LT

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

450 FOR J=1 TO NC

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

470 NEXT J

480 FOR I=1 TO N

490 IF I=K THEN GOTO 530

500 FOR J=1 TO NC

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

520 NEXT J

530 NEXT I

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

550 FOR I=1 TO N

560 FOR J=1 TO NC

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

580 NEXT J

590 NEXT I

600 NEXT K

610 RETURN

620 END

630 DATA 1,-3,2,-5,6

640 DATA -3,2,-1,2,-3

650 DATA 2,-1,4,-2,11

660 DATA -5,2,-2,6,-15