void cisi ( double x, double &ci, double &si );