//DFP.H
#if !defined(_DFP_H)
#define _DFP_H
#include "matrix.h"
#include"differ.h"
#include<fstream.h>

class CDfp
{
	public:
		CDfp(){}
		CMatrix operator()(double(*p)(CMatrix&),
			CMatrix& X0,double e);               //DFP
	private: 
	    double(*pfun)(CMatrix&);                 //pfun---目標函數指針
		CMatrix X;                               //X------迭代點(列向量)
        CMatrix S;                               //S------迭代方向
		double error;                            //error--誤差
		int n;                                   //n------維數，

		double explore(double a0,CMatrix TiDu);  //一維探索
		double explore0(double a,double b);      //黃金分割法       
		double f(double a);                      //構造求a的目標函數
        double find_a(CMatrix TiDu);             //求步長a 
};
#endif