class matrix
{
   double a[N][M];

public:

   double * operator[](size_t r)
   {
      return a[r];
   }

   double const * operator[](size_t r) const
   {
      return a[r];
   }
};