+ double betacf(const double, const double, const double);
+ double betai(const double, const double, const double);
+ double gammln(const double);
+ double gammp(const double, const double);
+ double gammq(const double, const double);
+ double gser(double&, const double, const double, double&);
+ double gcf(double&, const double, const double, double&);
+ double erfcc(double);
+
+ double ran0(int&); //for testing
+ double ran1(int&); //for testing
+ double ran2(int&); //for testing
+ double ran3(int&); //for testing
+ double ran4(int&); //for testing
+ void psdes(unsigned long &, unsigned long &); //for testing
+
+ void ludcmp(vector<vector<double> >&, vector<int>&, double&);
+ void lubksb(vector<vector<double> >&, vector<int>&, vector<double>&);
+
+ void ludcmp(vector<vector<float> >&, vector<int>&, float&);
+ void lubksb(vector<vector<float> >&, vector<int>&, vector<float>&);
+
+