subroutine dasol(al,au,ad,b,jp,neq,energy)
      real*8 al(*),au(*),ad(*),b(*),zero,energy,bd,dot
      do 100 is=1,neq
         if(b(is).ne.zero) go to 200
 100  continue
      return
 200  if(is.lt.neq) then
      endif
      do 400 j = is,neq
         energy=energy+bd*b(j)
 400  continue
      if(neq.gt.1)then
      endif
      end