#include <iostream>
#include <cmath>

using namespace std;

int main(){
int n, v;//A - матрица
cout<<"v="; cin>>v;
cout<<"n="; cin>>n;
int **A=new int *[n];

cout << n << endl;


for (int i=0; i<n; ++i)
  A[i]=new int [n];

A[0][0]=v;


for (int i=1; i<n; i++)
  A[i][i]=A[i-1][i-1]+2;


for (int i=0; i<n; i++)
  for (int j=0; j<n; j++)
    if(j!=i)
      A[i][j]=0.01*A[i][i];

cout << "Первичная матрица:" << endl;
for (int i=0; i<n; i++){
  for(int j=0; j<n; j++)
    if(A[i][j] >= 10)
      cout<<"A["<<i<<"]["<<j<<"] = " << A[i][j] << "| ";
    else
      cout<<"A["<<i<<"]["<<j<<"] = " << A[i][j] << " | ";
  cout << endl;
}

int *h=new int[n];//h-вектор-столбец

for (int i=0; i<n; i++)
  h[i]=A[i][i];

int *b=new int [n];//умножение матрицы А на вектор-столбец

for(int j=0;j<n;j++){
  b[j]=0;
  for(int i=0;i<n;i++){
    b[j]=b[j]+h[j]*A[j][i];
  }
}

cout << endl << "b: ";
for(int j=0;j<n;j++)
    cout << b[j] << " ";
cout << endl << endl;


int **M=new int *[n];//Расширенная матрица

for (int i=0; i<n; ++i)
  M[i]=new int [n];

for (int i=0; i<n; ++i){
  for (int j=0; j<n; ++j){
      M[i][j] = A[i][j];
      M[i][n] = b[i];
    }
}

cout << "Расширенная матрица: " << endl;
for (int i=0; i<n; i++){
  for(int j=0; j<=n; j++)
    if(M[i][j] >= 10)
      cout<<"M["<<i<<"]["<<j<<"] = " << M[i][j] << "| ";
    else
      cout<<"M["<<i<<"]["<<j<<"] = " << M[i][j] << " | ";
  cout << endl;
}

//прямой ход метода Гаусса

double d, s;
double *x = new double [n];


/*
for (int k = 0; k < n; k++){
    for (int j = k + 1; j <= n; j++){
        d = M[k][j] / M[k][k];
        for (int i = k; i <= n; i++){
            M[k][i] = M[k][i] - d * M[k][i];
        }
        b[k] = b[k] - d * b[k];
    }
}



for (int k = n; k >= 1; k--){
  d = 0;
  x[k] = (b[k] - d) / M[k-1][k];
  for (int j = k; j <= n; j++){
    s = M[k][j] * x[j];
    d = d + s;
  }
}

cout << endl << "Korni sistemy: " << endl;

for(int i = 0; i <= n; i++)
    cout << "x[" << i << "]=" << x[i] << " " << endl;
*/
}