A059831 Determinant of Wilkinson's eigenvalue test matrix of order 2n+1.

0, -2, -4, -20, -252, -5610, -187944, -8760272, -540315160, -42535905530, -4158250120140, -493938370048692, -70043599065844404, -11684931733886455730, -2264985487707963662992, -504752888883221450120000, -128137017404994234514023024

Offset

0,2

References

P. Jansen and P. Weidner, High-accuracy arithmetic software - some tests of the ACRITH problem-solving routines, ACM Trans. Math. Software, 12 (1986), 62-70.

MathWorks, MATLAB Reference Guide, 1996.

Links

Table of n, a(n) for n=0..16.

Example

The matrix of order 5:
2 1 0 0 0
1 1 1 0 0
0 1 0 1 0
0 0 1 1 1
0 0 0 1 2

Maple

with(LinearAlgebra): unprotect(D):
for n from 1 to 7 do N:=2*n+1; D:=Matrix(N, N, 0 );
for i from 1 to N-1 do D[i, i+1]:=1; D[i+1, i]:=1; od: for i from 1 to N do D[i, i]:=abs(n-i+1); od: lprint(Determinant(D)); od:

Prog

(MATLAB) for i = 0:7 a(i) = det(wilkinson(2*i+1)); end; a

Crossrefs

Cf. A059834.

Sequence in context: A120388 A061348 A127103 * A064493 A145614 A145626

Adjacent sequences: A059828 A059829 A059830 * A059832 A059833 A059834

Keyword

sign

Author

N. J. A. Sloane, Feb 25 2001

Extensions

More terms from David Wasserman, May 24 2002

Status

approved