Ubuntu-gr Forum

Η μέθοδος Jacobi με την ορίζουσα δεν είναι; Τέλος πάντων όταν έδινα πανελλήνιες και έπρεπε να λύνω καμιά 20ρια συστήματα τη μέρα είχα κάνει ένα πρόγραμμα για να μου δίνει το αποτέλεσμα με τη μέθοδο του Gauss. Το είχα γράψει σε ST BASIC και μάλλον το πρώτο κομμάτι είναι αυτό που χρειάζεσαι. Όπως και να χει δες το αν σε βοηθήσει:

Κώδικας: Επιλογή όλων

1 REM D
2 FULLW 2:CLEARW 2
3 DIM A(3,4)
10 FOR I=1 TO 3
20 FOR J=1 TO 4
25 PRINT "Εισάγετε το στοιχείο ";J;" της γραμμής ";I;"για σύστημα της μορφής ax+by+cz=d"
30 INPUT A(I,J)
40 NEXT J
50 NEXT I
60 K=A(2,1)/A(1,1):L=A(3,1)/A(1,1)
70 FOR M1=1 TO 4
80 A(2,M1)=A(2,M1)-K*A(1,M1)
90 NEXT M1
100 FOR M2=1 TO 4
110 A(3,M2)=A(3,M2)-L*A(1,M2)
120 NEXT M2
130 K=A(3,2)/A(2,2)
140 FOR M3=2 TO 4
150 A(3,M3)=A(3,M3)-K*A(2,M3)
160 NEXT M3
170 K=A(2,3)/A(3,3):L=A(1,3)/A(3,3)
180 FOR M4=3 TO 4
190 A(2,M4)=A(2,M4)-K*A(3,M4)
200 NEXT M4
210 FOR M5=2 TO 4
220 A(1,M5)=A(1,M5)-L*A(3,M5)
230 NEXT M5
240 K=A(1,2)/A(2,2):A(1,2)=A(1,2)-K*A(2,2):A(1,4)=A(1,4)-K*A(2,4)
250 X=A(1,4)/A(1,1):Y=A(2,4)/A(2,2):Z=A(3,4)/A(3,3)
260 FULLW 2:CLEARW 2
270 PRINT "X=";X;:PRINT "Y=";Y;:PRINT "Z=";Z
280 INPUT AN$