extern fint dpfpfl_c(fint *,double *,double *,fint *);