#include <iostream>
using namespace std;


double **init_a(int n){
        double **a=new double*[n];
        for (int i=0; i<n; i++) 
           a[i]=new double[n];
        for (int i=0; i<n;i++){
            for (int j=0;j<n;j++) {
             a[i][j]= (i==j)?1:0;
             }
         }
             return a;
             }  


void print_a(double **a, int n) {
     for (int i=0;i<n;i++){
         for (int j=0;j<n;j++){
             cout << a[i][j]<<"\t";
         }
         cout<< endl;
    }
}


double **enlarge(double **a,int n, int m) {
        if (m<=n){
              return a;
        }else {
             a=(double**)realloc(a,sizeof(double *)*m); 
             for (int i=0; i<m;i++) {
             a[i]=(double *)realloc(a[i],sizeof(double)*m);
        }
        for (int i=0;i<m;i++) {
             for (int j=n;j<m;j++) {
                a[i][j]=(i==j)?1:0;
              }
        }
      } return a;        
} 


double **mmul(double **a, double**b, int n) {
            double **c;
            c=new double*[n];
            for ( int i=0;i<n;i++) {
                  c[i]=new double[n];
           }
           for (int i=0;i<n;i++) {
                for (int j=0;j<n;j++) {
                    for (int k=0;k<n;k++) {
                        c[i][j]+=a[i][k]*b[k][j];
                    }
              }
          }
         return c;
   }      

int main (void) {
int n;
    cout<< "dose n;" ;
    cin>>n ; 
    
   double **a,**b,**c;
  a=init_a(n);
  print_a(a,n);
   return 0;
}