1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
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 enddo bar(k,l)=bar(k,l)+(v3b-1.d0) enddo enddo do m=1,ne do k=1,l l = l*(v3b**(-coefm)) enddo enddo enddo end