program Calling
implicit none
integer i,n,ier
integer, parameter :: nmax=10
real*8, dimension(nmax,nmax) :: a
real*8, dimension(nmax) :: spectrum,dummy

! Solving:  A*C = E*C 

!Example for 2x2 matrix 

n=2

! Generating the A matrix

a(1,1)=-2.d0 ;a(1,2)=-1.d0 ;a(2,1)=-1.d0 ; a(2,2)=-2.d0 

! Calling the subroutine
 
call EISRS1(Nmax,N,spectrum,a,IER,dummy)

print*,' You should get the following eigenvalues '
print*,' -3.00000000000 -1.0000000000 '

print*,' But the subroutine returned: '
print*, (spectrum(i),i=1,n)

stop
end  program Calling