Matriz identidad de orden 'n'

#include<iostream>
using namespace std;
int main()
{
 int ord;
 cout<<"Ingrese el orden de su matriz cuadrada: ";
 cin>>ord;
 cout<<"**************************************"<<endl;
 cout<<"La matriz identidad de orden "<<ord<<" es: "<<endl;

 int a[ord][ord],i,j;

 for(i=0;i<ord;i++){
  for(j=0;j<ord;j++){
   if(i==j){
    a[i][j]=1;   
   }else{
    a[i][j]=0;   
   }
   cout<<a[i][j]<<" ";
    }
    cout<<endl;
   }
   return 0;
}