//---------------------------------------------------------------------------
#ifndef RecipesH
#define RecipesH

#include <recipes/nr.h>
#include <recipes/nrutil.h>

#define NUMREC 1

#define SIGN(a,b) ((b) >= 0.0 ? fabs(a) : -fabs(a))
#define SWAP(a,b) temp=(a);(a)=(b);(b)=temp;
//---------------------------------------------------------------------------
//int ndatat;
//float *xt,*yt,aa,abdevt;extern int ndatat;extern float *xt,*yt,aa,abdevt;
void    medfit(float x[], float y[], int ndata, float *a, float *b, float *abdev);float   rofunc(float b);
float   select(unsigned long k, unsigned long n, float arr[]);

int     ncom;
float   *pcom,*xicom,(*nrfunc)(float []);
float   fcircfit(float p[]);
void    fcircderiv(float p[],float xi[]);
void    sort2fi(unsigned long n, float arr[], int brr[]);
void    sort3fff(unsigned long n, float arr[], float brr[], float crr[]);

#endif