Add stable_check.py

stable_check.py will do two things:

1.  Replicate the monte_carlo portion of the fitting algorithm (samples
a region for positive definiteness of the tangent stiffness), and

2.  Use logistic regression to find the decision boundary in strain
space between strain parameters which give stable models and strain
parameters which don't