Code covered by the BSD License
by Chebfun Team
30 Apr 2009
(Updated 05 Apr 2012)
Numerical computation with functions instead of numbers.
This file was selected as MATLAB Central Pick of the Week
function pass = grammatrix
% Construct the Gram matrix with chebfuns and
% compare it to Hilbert matrix.
A = chebfun;
x = chebfun(@(x) x,[0,1]);
A(:,n) = x.^(n-1);
G = A'*A;
pass = norm(G-hilb(4)) < 10*eps;