void divset(long,long[],long,long,double[]);
void dnlsgb(long,long,long,double[],double[][2],double[],double[],
          void(*)(long,long,long,double[],long*,double[]),
          void(*)(long,long,long,double[],long*,double[]),long*,long,long[],long,long,double[]);
void drnsgb(double*,double[],double[][2],double[],double*,long[][2],long[],
	long,long,long,long,long,long,long,
	long,double[],double[]);