#ifndef __REVERSE_TASK_SOLVER__

#define __REVERSE_TASK_SOLVER__



class CReverseTask_solver

{

	public:

		CReverseTask_solver();

		~CReverseTask_solver();



		void solve_task();



		double (*phi1)(double, double);			//????????? ??????? ?? ????? ???????

		double (*phi2)(double, double);			//????????? ??????? ?? ?????? ???????

		double (*psi)(double, double);			//????????? ???????

		double (*f)(double, double);			//??????????????

		double (*analytic)(double, double);		//????????????? ???????



		double alpha1;		// alpha1 * Ut + betta1 * U = phi1(x, t)

		double betta1;

		double alpha2;		// alpha2 * Ut + betta2 * U = phi2(x, t)

		double betta2;



		double xb;		// x begin

		double xe;		// e end

		double maxT;	// maximum Time



		int N;			// ?????????? ????? ? ????????? ?? ??? X

		int K;			// ?????????? ????? ?? ????????? ?????????



		double stepT;	// tau - ??? ?? ???????



		double a;

		double b;

		double c;



		double accuracy;	// ????????

		double task_c;

		double task_lambda;

		double task_qmax;

		double task_tmax;



		double* q_result;	// ??????? ?????? q

		int iter_num;	//iteration number



		double** directTable;	//??????? ??????? ?????? (??????? 1-?? ????)



	private:

	protected:

		double get_discrepancy();	// ???????? ??????? (?????????? ??????? ???????????)

};



#endif