void dbi0k0(double,double*,double*,long,long*);
void dbi1k1(double,double*,double*,long,long*);