void dfft(double[],double[],double[]);
void drft1(double[],byte,long,long*,double[]);
void ermsg(STRING,long,long,STRING,byte);
void ierm1(STRING,long,long,STRING,STRING,long,byte);