A three-dimensional analytical technique for the representation of the cranial surface based on a spherical harmonic expansion is presented. The analytical surface representation is used in conjunction with a sample iterative ray-intersection algorithm to compute the dose calculation depth. The technique is compared to other conventional techniques and it is determined that the spherical harmonic approach of approximating the cranial surface decreases the depth calculation time by approximately a factor of 60.