|
|
MAPLE
|
with(LinearAlgebra):
a:= proc(n) option remember; local d, i, j, t, m, M;
d:= 2*n; m:= d^2-4;
M:= Matrix(m, shape=skewsymmetric);
for i to d-3 do M[i+1, i]:=1 od;
for i to d-2 do M[i, i+d-1]:=1 od;
for i from m-d+3 to m-1 do M[i, i+1]:=1 od;
for i from m-d+3 to m do M[i-d+1, i]:=1 od;
for i from d-1 to m-2*d+2 do M[i, i+d]:=1 od;
for i to d-2 do for j to d-1 do
t:=d*i+j-2; M[t, t+1]:= `if`(irem(i, 2)=1, 1, -1);
od od;
isqrt(Determinant(M))
end:
|
|
|
MATHEMATICA
|
a[1] = 1; a[n_] := a[n] = Module[{d, i, j, t, m, M}, d = 2*n; m = d^2 - 4; M = Array[0&, {m, m}];
For[i = 1, i <= d - 3, i++, M[[i + 1, i]] = 1];
For[i = 1, i <= d - 2, i++, M[[i, i + d - 1]] = 1];
For[i = m - d + 3, i <= m - 1, i++, M[[i, i + 1]] = 1];
For[i = m - d + 3, i <= m, i++, M[[i - d + 1, i]] = 1];
For[i = d - 1, i <= m - 2*d + 2, i++, M[[i, i + d]] = 1];
For[i = 1, i <= d - 2, i++,
For[j = 1, j <= d - 1, j++, t = d*i + j - 2; M[[t, t + 1]] = If[Mod[i, 2] == 1, 1, -1]]]; M = M - Transpose[M]; Sqrt[Det[M]]];
|
