The three-dimensional (3-D) VMEC code has been modified to model an energetic species with a variant of a Bi-Maxwellian distribution function that satisfies the constraint B[nabla][script F]h = 0, and the 3-D TERPSICHORE stability code has been extended to investigate the effects of pressure anisotropy in two limits. The lower limit is based on a purely fluid Kruskal-Oberman (KO) energy principle (ignoring the stabilizing kinetic integral), and the upper limit is obtained from an energy principle in which the hot particle pressure and current density refrain from interacting with the dynamics of the instability because their diamagnetic drift frequency is considered much larger than the dominant growth rate. We have specifically investigated the instability properties of a Heliotron device with a major radius of 3.9 m and total <> [approximately equal to] 3.9%, where the energetic particle contribution <h> varies from 0 to 1.3% for T