In most cases, nonlinearities from magnets must be properly included in tracking and analysis to properly compute quantities of interest, in particular chromatic properties and dynamic aperture. One source of nonlinearities in magnets that is often important and cannot be avoided is the nonlinearity arising at the end of a magnet due to the longitudinal variation of the field at the end of the magnet. Part of this effect is independent of the longitudinal of the end. It is lowest order in the body field of the magnet, and is the result of taking a limit as the length over which the field at the end varies approaches zero. This is referred to as a 'hard edge' end field. This effect has been computed previously to lowest order in the transverse variables. This paper describes a method to compute this effect to arbitrary order in the transverse variables, under certain constraints.