void dcalcr(long,long,double[],long*,double[]);
void dcalcj(long,long,double[],long*,double[]);
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[]);