OFFSET
0,3
LINKS
Wikipedia, Toeplitz Matrix
EXAMPLE
a(3) = 0:
[ 3, 2*i, i;
-2*i, 3, 2*i;
-i, -2*i, 3 ]
MAPLE
A359561 := proc(n)
local T, c, r ;
if n =0 then
return 1 ;
end if;
T := Matrix(n, n) ;
T[1, 1] := n ;
for c from 2 to n do
T[1, c] := (n-c+1)*I ;
end do:
for r from 2 to n do
for c from 1 to r-1 do
T[r, c] := -T[c, r] ;
end do:
for c from r to n do
T[r, c] := T[r-1, c-1] ;
end do:
end do:
LinearAlgebra[Determinant](T) ;
simplify(%) ;
end proc:
seq(A359561(n), n=0..25) ; # R. J. Mathar, Jan 31 2023
MATHEMATICA
Join[{1}, Table[Det[ToeplitzMatrix[Join[{n}, I Reverse[Range[n-1]]]]], {n, 19}]]
PROG
(Python)
from sympy import Matrix, I
def A359561(n): return Matrix(n, n, [(n+j-i if i>j else j-i-n) if i!=j else n*I for i in range(n) for j in range(n)]).det()*(1, -I, -1, I)[n&3] # Chai Wah Wu, Jan 25 2023
CROSSREFS
KEYWORD
sign
AUTHOR
Stefano Spezia, Jan 06 2023
STATUS
approved