Code covered by the BSD License
30 Apr 2009
24 Jun 2014)
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;