double sum_square(double (*f)(double), int m, int n)
{
   .....