void dcalca(long,long,long,double[],long*,double[]);
void dcalcb(long,long,long,double[],long*,double[]);
void divset(long,long[],long,long,double[]);
void dnlsfb(long,long,long,double[],double[][2],double[],double[],
      void(*)(long,long,long,double[],long*,double[]),long*,long,long[],
      long,long,double[]);