void divset(long,long[],long,long,double[]);
void dnlsgu(long,long,long,double[],double[],double[],
          void(*)(long,long,long,double[],long*,double[]),
          void(*)(long,long,long,double[],long*,double[]),long*,long,long[],long,long,double[]);
void drnsg(double*,double[],double[],double*,long[][2],long[],long,
	long,long,long,long,long,long,long,
	double[],double[]);