function funkval(x,y,la,N)
!
!     Tar inn en x og y verdi, regner saa ut tilhoerende 
!     funksjonsverdi, skal nyttes for aa bestemme steglengde
!     i steepest descent iterajonen.
!
      use Interface, ONLY: leg_val
!
      double precision :: funkval
      integer, intent(in)          :: N
      double precision, intent(in) :: x,y
      double precision, intent(in) :: la(0:N-1,0:N-1)
!
      integer i,j
      double precision test,xval(0:N-1),yval(0:N-1)
!
      xval = leg_val(x,N)
      yval = leg_val(y,N)
!
      test = 0.d0
        do i = 0,N-1
         do j = 0,N-1
          test = test+la(i,j)*xval(i)*yval(j)
         end do
        end do
!
      funkval = test
!
      end function funkval