double dcpval(double[],long,double);
double dsi(double);