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