function y = f2(x)
% functional form for Example 9.3

y = [ones(1,length(x)); x; x.^2; x.^3];