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