Hello,
how can I write the C++ code about 61-point Gauss-Kronrod Quadrature?

Take a look at function QK61 in the (Fortran, I'm afraid, but it's not terribly hard to read) code at

http://people.sc.fsu.edu/~burkardt/f_src/quadpack/quadpack.f90