We present a systematic investigation of the distribution of normal forces at the boundaries of static packings of spheres. A new method for the efficient construction of large hexagonal-close-packed crystals is introduced and used to study the effect of spatial ordering on the distribution of forces. Under uniaxial compression we find that the form for the probability distribution of normal forces between particles does not depend strongly on crystallinity or inter-particle friction. In all cases the distribution decays exponentially at large forces and shows a plateau or possibly a small peak near the average force but does not tend to zero at small forces.