#ifndef _EULER_H_
#define _EULER_H_
 
#include <vector>
//#include <cmath>
 
typedef std::vector< double > vec;
typedef double (*functPtr)( double, vec );
 
class Euler
{
private:
	double timeStep;
	std::vector< double > initialConditions;
	std::vector< double > solution;
 
public:
	Euler();
	vec solve(std::vector< functPtr > functions);
	void setInitial(std::vector< double > initial);
	void setTimeStep(double time);
};
 
#endif