#ifndef Functions_h
#define Functions_h

void deriv_matrix(double ** mat, double ** der_mat, int direction, int low1, int high1, int low2, int high2, int odd);

void Jacobian (double **JJ, double dTC, double dmuBC, double angle1, double angle2);

double rr(double T, double muB);

double hh(double T, double muB);

void funcv(int n,double x[],double f[]);

double G (double R, double Theta);

double dGdmuBConT (double R, double Theta);

double dGdTConmuB (double R, double Theta);

double d2GdmuB2ConT (double R, double Theta);

double d2GdT2ConmuB (double R, double Theta);

double d2GdTdmuB (double R, double Theta);

double d3GdmuB3ConT (double R, double Theta);

double d4GdmuB4ConT (double R, double Theta);

double M(double R, double Theta);

double z(double R, double Theta);

double gEpsExp(double R, double Theta);

double gAsymp(double R, double Theta);

double CalcCorrLengthEpsExp(double R, double Theta);

double CalcCorrLengthAsympt(double R, double Theta);


double dMdhConr (double R, double Theta);

double GausFunc (double x, double y, double sigmax, double sigmay);

void GausFilter2D (double **data, double sigmax, double sigmay, int lx, int ly);



#endif /* Functions_h */