Heisenman-Drogulous isomoglobic algorithm