1 2 3 4 5 6 7 8 9 10 11 12 13 14
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