Stresses from arbitrary loads on a penny-shaped crack

An implementation of a solution to the problem of a penny-shaped crack in an infinite elastic solid with arbitrary normal and shear loads is described, and is used to generate the stresses corresponding to some simple crack loads. The program described is fast and stable, and is shown to give accurate results even if the crack loads are not of the desired polynomial form. Stress intensity factors are obtained directly from combinations of load constants.