1function C4 = C4f(epsi, C4x)
2%C4F  Evaluate C_4
3%
4%   C4 = C4F(epsi, C4x) evaluates C_{4,l} in the expansion for the area
5%   (Eq. (65) expressed in terms of n and epsi) using the coefficient
6%   vector C4x.  epsi is a K x 1 array.  C4x is a 1 x 21 array.  C4 is a
7%   K x 6 array.
8
9  nC4 = 6;
10  C4 = zeros(length(epsi), nC4);
11  mult = 1;
12  o = 1;
13  for l = 0 : nC4 - 1
14    m = nC4 - l - 1;
15    C4(:, l+1) = mult .* polyval(C4x(o : o + m), epsi);
16    o = o + m + 1;
17    mult = mult .* epsi;
18  end
19end
20