Decomposing a DFT into multiple FFT calls
The simplest thing to do if you have access to only a power-of-two FFT is to zero-pad your data going into the $N=2^p$ FFT and interpolate what comes out.  i think that's what MATLAB does for fft() for non-power-of-two DFT.