This report documents a finite-element analysis procedure for computation of Mode I and Mode II stress intensity factors at a cracked fastener hole in a single row or a double staggered row of fastener holes on a tension panel. The panel is loaded by cosine bearing pressure distributions applied to each fastener hole. Different bearing forces may be specified for each hole, and any hole may be designated as the damaged hole. The physical structure, finite- element model, program flow and options, input/output conventions, execution times, and limitations are discussed. Results of some example analyses are presented.