subroutine foo(bar)
real*8 bar(3,3),coefm
do ii=istart,iend
do i=1,21
bar(k,l)=4
enddo
do m=1,ne
do l=1,3
do k=1,l
bar(k,l)=bar(k,l)+(v3b-1.d0)
l = l*(v3b**(-coefm))
end