void djacg(long*,long,long,double[],double[],double*,long,double[],double[],
      long[],double[],long,long[],long);