double dpsi(double);
void dpsie(double*,long*);