void divset(long,long[],long,long,double[]);
void dnlagb(long,long,double[],double[][2],
          void(*)(long,long,double[],long*,double[]),
          void(*)(long,long,double[],long*,double[]),long[],long,long,double[]);
void drn2gb(double[][2],double[],double*,long[],long,long,long,
	long,long*,long*,long,double[],double[],double[],
	double[]);