program test
C
C  Solve A*U = b      
C        [2   -2    0    0   -1]
C        [0    3   -2    0    0]
C  A =   [0    0    5   -3    0]
C        [0    0    0   10    4]
C        [0    0    0    0   10] 
C
C The solution, U, is: 636.,619.,292.,74.,34.
C
        parameter (neqs = 5)
        real*8 au(7), ad(neqs), b(neqs), energy
        integer jp(neqs)
        data au/-2.,-2.,-3.,-1.,0.,0.,4./
        data ad/2.,3.,5.,10.,10./
        data b/0.,1.,0.,0.,0./
        data jp/0,1,2,3,7/
        call datri(0, au, ad, jp, neqs, neqs)
        call dasol(0, au, ad, b, jp, neqs, neqs, energy)
        write(6,*) 'solution', (b(i), i=1,neqs)
      end