function [A] = array_rstar(m)

if nargin < 1 || isempty(m), m=3 ; end


v1= [1:m,m-1:-1:1];
v2= [1:m-1,0,m-1:-1:1];

A1 = diag(v1);
A2 = fliplr(diag(v2));
A = A1 + A2;