!-- showing how to uses Fortran for vectored operations !-- equations work on vectors, no need for loop program f08_a implicit none integer, parameter :: N = 7 real, parameter :: D(N) = [-0.2,1.0,1.5,3.0,-1.0,4.2,3.1] real,parameter :: H(N) = [2.1,2.4,1.8,2.6,2.6,2.2,1.8] real, parameter :: pi = 2.0* ACOS(0.0) real, allocatable :: V(:) V = (1.0/12.0)*pi*(D**2)*H print *, v end program f08_a