матрица

int i,s1,s2,n1,n2;
 
double sum(double a) { return a+1; }
 
template <typename T, size_t Rows>
struct avg : public std::unary_function<const T (&)[Rows], double>
{
        T operator()(const T (&arr)[Rows]) const
        {
                return std::accumulate(arr, arr + Rows, T()) / Rows;
            
        }
};
 
 
 
template <typename T, size_t Cols, size_t Rows>  
void go(T (&m)[Cols][Rows])
{
        double maxel;
std::vector<T> v(Cols);
 maxel=*std::max_element(v.begin(),
                std::transform(m, m +Cols , v.begin(), avg<T, Rows>()));
    std::cout << "\n Максимальный среднее арифметическое=" <<maxel ;
 
        
      std::transform(
                &m[s1][0], &m[s1][0] + Rows*n1, &m[s1][0],//Cols
                std::bind2nd(std::minus<T>(), maxel)
                );
     
 
      if (i==2) {
      std::transform(
                &m[s2][0], &m[s2][0] + Cols*n2, &m[s2][0],
                std::bind2nd(std::minus<T>(), maxel)
                );
      }
   }
 
 
template <typename T, size_t Rows>
struct X : public std::unary_function<const T (&)[Rows], void>
{
        X(std::ostream& os) : os_(os) {}
        
        void operator()(const T (&arr)[Rows]) const
        {
                typedef std::ostream_iterator<T> O;
                std::copy(arr, arr+Rows, O(os_, "\t"));
                os_ << std::endl;
        }
 
        std::ostream& os_;
};
 
 
template <typename T, size_t Cols, size_t Rows>
std::ostream& operator<< (std::ostream& os, const T (&m)[Cols][Rows])
{
        std::for_each(m, m+Cols, X<T, Rows>(os));
        return os;
}