1function C3 = C3f(epsi, C3x)
2%C3F  Evaluate C_3
3%
4%   C3 = C3F(epsi, C3x) evaluates C_{3,l} using Eq. (25) and the
5%   coefficient vector C3x.  epsi is a K x 1 array.  C3x is a 1 x 15 array.
6%   C3 is a K x 5 array.
7
8  nC3 = 6;
9  C3 = zeros(length(epsi), nC3 - 1);
10  mult = 1;
11  o = 1;
12  for l = 1 : nC3 - 1
13    m = nC3 - l - 1;
14    mult = mult .* epsi;
15    C3(:, l) = mult .* polyval(C3x(o : o + m), epsi);
16    o = o + m + 1;
17  end
18end
19