@Rouben Rostamian  Thanks for the comment! You are right that the classification of the points was not done with the constraint in mind. The optimization package commands minimize and maximize were used in one of the versions of the script to classify points with constraints, but this was ommitted in the version posted. I have updated the document to simplify the calculator's operation to limit it to only giving the form of output that the command LagrangeMultipliers returns. Since the method of classifying the points using optimization (as was present in a previous version) would involve a convoluted and complicated method of catching warnings, I thought it would be better to focus the document around the LagrangeMultipliers command.