Four subroutines compute one-dimensional and multi-dimensional Fourier transforms for real data, multi-dimensional complex Fourier transforms, and multi-dimensional sine, cosine and sine-cosine transforms. Subroutines use Cooley-Tukey fast Fourier transform. In all but one-dimensional case, transforms are calculated in up to six dimensions