#include <math.h>
#include <iostream>

using namespace std;
int main()
{
    double a[6][6];
    double b[6][6];
    int c[6][6];
    double d[6][6];
    double to_swap;
    bool swapped;
    double e[6][6];
    for(int j=0; j<6; j++)
    {
        for(int i=0; i<6; i++)
        {
            a[i][j]=((i+1)+(j+1))*(cos(j+1)/sin(j+1));
        }
    }

    for(int j=0; j<6; j++)
    {
      for(int i=0; i<6; i++){
        b[i][j]=a[i][j];
      }
    }
    for(int j=0; j<6; j++){
      swapped == false;
      do{
        swapped = false;
        for(int i = 1; i < 6; i++){
          if(b[i-1][j] < b[i][j]){
          swapped = true;
          to_swap=b[i][j];
          b[i][j] = b[i-1][j];
          b[i-1][j] = to_swap;
          }
        }
        } while(swapped == true);
    }

    for(int j=0; j<6; j++){
      int x=1;
      for(int i=0; i<6; i++){
        if(i<j)
          c[i][j]=0;
        else{
          c[i][j]=x;
          x++;
        }
      }
    }

    for(int i = 0; i < 6; i++)
      for(int j = 0; j < 6; j++)
        for(int k = 0; k < 6; k++)
          {
        d[i][j]+=b[i][k]*c[k][j];
          }

    for(int i=0; i<6; i++){
      for(int j=0; j<6; j++){
        if((i<2) && (j<2)){
          e[i][j]=d[i+2][j+2];
        } else if((i<4) && (i>1) && (j<4) && (j>1)){
          e[i][j]=d[i-2][i-2];}
        else{
          e[i][j]=d[i][j];
        }
      }
    }

    int i_max_0=e[0][0];
    for(int i=0; i<6; i++){
      if(e[i][0] > i_max_0){
          i_max_0 = e[i][0];
        }
    }
    int i_max_4=e[0][4];
    for(int i=0; i<6; i++){
      if(e[i][4] > i_max_4){
          i_max_4 = e[i][4];
        }
    }

    cout << endl << endl << endl << "matrix A" << endl;
    for(int j=0; j<6; j++){
      for(int i=0; i<6; i++){
        cout << a[i][j] << ' ';
      }
      cout << endl;
    }

    cout << endl << endl << endl << "matrix B" << endl;
    for(int j=0; j<6; j++){
      for(int i=0; i<6; i++){
        cout << b[i][j] << ' ';
      }
      cout << endl;
    }

    cout << endl << endl << endl << "matrix C" << endl;
    for(int j=0; j<6; j++){
      for(int i=0; i<6; i++){
        cout << c[i][j] << ' ';
      }
      cout << endl;
    }

    cout << endl << endl << endl << "matrix D" << endl;
    for(int j=0; j<6; j++){
      for(int i=0; i<6; i++){
        cout << d[i][j] << ' ';
      }
      cout << endl;
    }

    cout << endl << endl << endl << "matrix E" << endl;
    for(int j=0; j<6; j++){
      for(int i=0; i<6; i++){
        cout << e[i][j] << ' ';
      }
      cout << endl;
    }

    cout << "max on row 1" << i_max_0 << endl << "max on row 5" << i_max_4;
    return 0;
}