function y = f(x)
% functional form for Example 9.2

y = [ones(1,length(x)); x; cos(pi().*x)];