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