dddddd.bas

J-O promgram 用于稀土研究工作者

8000 GOTO 8064
8002 PRINT
8004 FOR S = 1 TO N
8006  FOR T = S TO N
8008     IF A(T, S) < 0 THEN GOTO 8016
8010  NEXT T
8012  PRINT "NO UNIQUE SOLUTION"
8014  GOTO 8062
8016  GOSUB 8036
8018  LET C = -A(S, S)
8020  GOSUB 8048
8022  FOR T = 1 TO N
8024    IF T = S THEN GOTO 8030
8026    LET C = -A(T, S)
8028    GOSUB 8056
8030  NEXT T
8032 NEXT S
8034 GOTO 8062
8036 FOR J = 1 TO N + 1
8038   LET B = A(S, J)
8040   LET A(S, J) = A(T, J)
8042   LET A(T, J) = B
8046 NEXT J
8046 RETURN
8048  FOR J = 1 TO N + 1
8050    LET A(S, J) = C * A(S, J)
8052  NEXT J
8054  RETURN
8056  FOR J = 1 TO N + 1
8058   LET A(T, J) = A(T, J) + C * A(S, J)