template<typename T, std::size_t N> 
struct gauss_coefficients {
  // types
  typedef std::tr1::array< T, N/2+1 > value_type;

  // public static functions
  static const value_type & weights() ;
};