Galīgo elementu metode

2.7. GEM UZDEVUMA RISINĀŠANAS ALGORITMS 39
1830 PRINT #2,”NODE#, APPLIED LOAD”
1840 FOR I=1 TO NL
1850 PRINT #2,ILD(I), AL(I) : NEXT I
1860 PRINT #2,”NODE#, DISPLACEMENT”
1870 FOR I = 1 TO NN
1880 PRINT #2, I, F(I) :NEXT I
1890 PRINT #2,”ELEMENT#, STRESS”
1900 FOR I=1 TO NE
1910 PRINT #2, I, STRESS(I) :NEXT I
1920 PRINT #2, ”SUPPORT#, REACTION FORCE”
1930 FOR I=1 TO ND
1940 PRINT #2,NU(I),REACT(I) :NEXT I
1950 CLOSE #2
1960 END
5000 N=NN
5010 REM ** FORWARD ELIMINATION **
5020 FOR K = 1 TO N - 1
5030 NBK = N - K + 1
5040 IF N - K + 1 > NBW THEN NBK = NBW
5050 FOR I = K + 1 TO NBK + K - 1
5060 I1 = I - K + 1
5070 C = S(K, I1) / S(K, 1)
5080 FOR J = I TO NBK + K - 1
5090 J1 = J - I + 1
5100 J2 = J - K + 1
5110 S(I, J1) = S(I, J1) - C * S(K, J2): NEXT J
5120 F(I) = F(I) - C * F(K)
5130 NEXT I: NEXT K
5140 REM ** BACK SUBSTITUTION **
5150 F(N) = F(N) / S(N, 1)
5160 FOR II = 1 TO N - 1
5170 I = N - II
5180 NBI = N - I + 1
5190 IF N - I + 1 > NBW THEN NBI = NBW
5200 SUM = 0!
5210 FOR J = 2 TO NBI
5220 SUM = SUM + S(I, J) * F(I + J - 1): NEXT J
5230 F(I) = (F(I) - SUM) / S(I, 1): NEXT II
5240 RETURN
Šo programmu var translēt ar jebkuru Basic kompilātoru, piemēram, Mi-