3D FINITE ELEMENTS WITH HARMONIC BASIS FUNCTIONS FOR APPROXIMATIONS OF HIGH ORDER