急求各位高手下面这个有源程序的错误（fortran）

sunqingsuo

这个程序是用fortran编写的，我初学有限元和fortran，所以请教各位，帮助指点一下

1.       SUBROUTINE DECOG(A,N,NX)
   
2. C  DECOMPOSE A SYMMETRIC MATRIX INTO AN UPPER TRIANGULAR MATRIX
   
3. C
   
4.       DIMENSION A(NX,NX)
   
5. C
   
6. C         
   
7.       IF(A(1,1))1,1,3
   
8.     1 WRITE(6,2)
   
9.     2 FORMAT('ZERO OR NEGATIVE RADICAND')
   
10.       GO TO 200
    
11.     3 A(1,1)=SQRT(A(1,1))
    
12.       DO 10 J=2,N
    
13.    10 A(1,J)=A(1,J)/A(1,1)
    
14.       DO 40 I=2,N
    
15.         I1=I-1
    
16.         D=A(I,I)
    
17.         DO 20 L=1,I1
    
18.    20 D=D-A(L,I)*A(L,I)
    
19.       IF(A(I,I))1,1,21
    
20.    21 A(I,I)=SQRT(D)
    
21.       DO 40 J=I2,N
    
22.         I2=I-1
    
23.         D=A(I,J)
    
24.         DO 3O L=1,I2
    
25.    40 A(I,J)=D/A(I,I)
    
26.    30 D=D-A(L,I)*A(L,J)
    
27.       DO 50 I=2,N
    
28.         I1=I-1
    
29.         DO 50 J=1,I1
    
30.    50 A(I,J)=0.
    
31. C
    
32.   200 RETURN
    
33.       END

复制代码