M代码的罚函数程序,包括DFP法以及进退法确定搜索区间。