Computing wave propagation is of the utmost importance in computational optics, especially three-dimensional optical imaging and computer-generated hologram. The angular spectrum method, based on fast Fourier transforms, is one of the efficient approaches; however, it induces sampling issues. We report a Hybrid Taylor Rayleigh-Sommerfeld diffraction (HTRSD) that achieves more accurate and faster wave propagation than the widely used angular spectrum method.