PDH discriminant calculations
run in iPython environment to show local variables
ipython
run PDH.py eps_pdh