请教

homage79

SUBROUTINE GAUSS(A,B,N)
      DIMENSION A(N,N),B(N)
      DO 1 I=1,N
      I1=I+1
      DO 10 J=I1,N
10    A(I,J)=A(I,J)/A(I,I)
      B( I )=B( I )/A(I,I)
      A( I,I )=1.0
      DO 20 J=I1,N
      DO 30 M=I1,N
30    A(J,M)=A(J,M)-A(J,I)*A(I,M)
20    B( J)=B(J)-A(J,I)*B( I )
1     CONTINUE
      DO 40 I=N-1,1,-1
      DO 50 J=I+1,N
50    B( I )=B( I )-A(I,J)*B(J)
40    CONTINUE
      RETURN
      END                                                
这段程序是A是刚度阵，B是外加载荷
其作用是什么啊？
请指教

yp51920

就是GAUSS消取法解方程：K*X=P
K：刚度矩阵；X：结点位移；P：外加载荷

homage79

嗯，多谢版主，俺已经明白了
多谢多谢