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