Phenotype phase plane analysis using interior point methods