% Dimensione
n = 5 ; 

% Allocazione con matrice nulla
A = zeros( n , n ) ; 

% Vettore dei valori diagonali
vi = 1 : n ; 

% Vettore degli indici lineari degli elementi diagonali
vk = 1 : (n+1) : n^2 ; 

% Assegnazione vettoriale dei termini diagonali
A(vk) = vi