Now including code that correctly handles constraints on alpha, beta, and gamma for CD with energy fits. This new code also creates graphs of variation of SE, lambda, alpha, (and therefore beta) for a grid search over gamma.