C#NUMPAC#JACOBS
C.... EXAMPLE OF JACOBS ....
      REAL*8 P
      DIMENTION A(64,64),V(64,64),E(64)
      KA=
      EPS=
      P=DATAN(1.D0)*2.D0/DFLOAT(N+N+1)
C#２次元配列A(正方行列)の次数Nの値を入力
      N=
      
      DO 4 J=1,N
    4 E(J)=0.25D0/DSIN(DFLOAT(J+J-1)*P)**2
      CALL JACOBS(A,KA,N,EPS,V,IL)
      S=0.
      DO 5 J=1,N
    5 S=ABS(E(J)-A(J,J))+S
      S=S/FLOAT(N)
      WRITE(6,100) IL,S
  100 FORMAT(3X'ILL ='I6,10X'ERROR= 'E12.3/)
      WRITE(6,200) (A(J,J),E(J),J=1,N)
  200 FORMAT(10X2E20.6)
      STOP
      END