double sq(double x);