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

20 REM * *

30 REM * Μέθοδος Jacobi *

40 REM * για γραμμικά συστήματα *

50 REM * *

60 REM ****************************************

70 CLS

80 DIM A(20,21),X(20),XN(20)

90 EPS=.0000005

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

110 NC=N+1

120 PRINT

130 FOR I=1 TO N

140 FOR J=1 TO N

150 READ A(I,J)

160 NEXT J

170 READ A(I,NC)

180 NEXT I

190 GOSUB 280

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

210 FOR I=1 TO N

220 PRINT "X(";I;")=";X(I)

230 NEXT I

240 END

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

260 REM Με τη μέθοδο του Jacobi

270 REM Ελεγχος της δαγώνιας υπεροχής

280 FOR I=1 TO N

290 PIV=0

300 FOR J=1 TO N

310 IF J=I THEN GOTO 330

320 PIV=PIV+ABS(A(I,J))

330 NEXT J

340 IF ABS(A(I,I))<PIV THEN GOTO 370

350 NEXT I

360 GOTO 390

370 PRINT "Δεν υπάρχει διαγώνια υπεροχή"

380 STOP

390 REM Επαναληπτική διαδικασία

400 FOR L=1 TO N

410 X(L)=0

420 NEXT L

430 FOR K=1 TO 100

440 FOR I=1 TO N

450 SUM=0

460 FOR J=1 TO N

470 IF I=J THEN GOTO 490

480 SUM=SUM+A(I,J)*X(J)

490 NEXT J

500 XN(I)=(A(I,NC)-SUM)/A(I,I)

510 NEXT I

520 IF ABS(XN(1)-X(1))<EPS THEN RETURN

530 FOR L=1 TO N

540 X(L)=XN(L)

550 NEXT L

560 NEXT K

570 PRINT "Σε 100 επαναλήψεις δεν έχουμε σύγκλιση"

580 END

590 DATA 4,1,0,14

600 DATA 1,4,1,12

610 DATA 0,1,4,6