c c ------------------- c subroutine axvecb2vecc(a,b,c) c transforms vector B by A to get vector C real a(3,3),b(3),c(3) c do 100, i=1,3 c(i)=0. do 100, j=1,3 c(i)=c(i)+a(i,j)*b(j) 100 continue return end c c ------------------- c