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