1 subroutine rota(alfa) 2 implicit double precision (a-h,o-z) 3 common/rotxyz/amat(3,3),bmat(3,3),cmat(3,3),rmat(3,3),todeg 4 5c-- alfa expected in degrees 6 rada = alfa / todeg 7 amat(1,1) = dcos(rada) 8 amat(1,2) = dsin(rada) 9 amat(2,1) = -dsin(rada) 10 amat(2,2) = dcos(rada) 11 12 return 13 end 14