implicit real*8(A-H,O-Z)  
      include 'param.h'
      PARAMETER(MJ=20)
      DIMENSION JCGA(NA,MJ)
      DIMENSION CGA(NA,MJ)
      DIMENSION SOLX(NA)
      DIMENSION BVEC(NA),IHEAD(NA+1),ICOL(NNONZERO),A(NNONZERO)

      N=1000
      DO 1 I=1,N
        DO 2 J=1,20
          JCGA(I,J)=0
      CGA(I,J)=0D0
 2    CONTINUE
        SOLX(I)=0D0
 1    CONTINUE

      DO 10 I=1,N

        IF (I-1.GT.0) THEN
        JCGA(I,4)=I-1
        CGA(I,4)=3D0
        ENDIF

        JCGA(I,5)=I
        CGA(I,5)=110.3d0

        IF (I+1.LE.N) THEN
        JCGA(I,6)=I+1
        CGA(I,6)=3d0
        ENDIF


        IF (I-79.GT.0) THEN
        JCGA(I,3)=I-79
        CGA(I,3)=-16d0
        ENDIF

        IF (I+79.LE.N) THEN
        JCGA(I,7)=I+79
        CGA(I,7)=-16d0
        ENDIF

        IF (I-80.GT.0) THEN
        JCGA(I,2)=I-80
        CGA(I,2)=-2d0
        ENDIF

        IF (I+80.LE.N) THEN
        JCGA(I,8)=I+80
        CGA(I,8)=-2d0
        ENDIF


        BVEC(I)=mod(I,100)/100D0
c       BVEC(I)=I

 10     CONTINUE

        ERR=1D-12

        KA=0
        NONZERO=0
        IHEAD(1)=1 
        DO 20 I=1,N
	    KA=0
	    DO 30 J=1,MJ
	      IF (JCGA(I,J).GT.0) THEN
             ICOL(KA+IHEAD(I))=JCGA(I,J)
	       A(KA+IHEAD(I))=CGA(I,J)
             KA=KA+1
             NZERO=NZERO+1
	ENDIF
 30     CONTINUE
        IHEAD(I+1)=IHEAD(I)+KA
 20     CONTINUE	   



        CALL PICCGRP(IHEAD,ICOL,A,BVEC,N,SOLX,ERR)


      RESBMA=0D0
      DO 410 I=1,N
        TEMP1=0D0
        DO 411 J=1,MJ
          IF (JCGA(I,J).GT.0) THEN
            TEMP1=TEMP1+CGA(I,J)*SOLX(JCGA(I,J))
          ENDIF
 411    CONTINUE
        RESB=ABS(BVEC(I)-TEMP1)
        RESBMA=MAX(RESBMA,RESB)

 410  CONTINUE

      WRITE(6,*) 'ZANSA ',RESBMA

      END