a0 = 1;
a1 = -1/(gamma(3+1)*3);
a2 = 1/(gamma(5+1)*5);
a3 = -1/(gamma(7+1)*7);
a4 = 1/(gamma(9+1)*9);

for i=1:11
  x(i) = 0 + (i-1)*1/10;
  u = x(i)^2;
  p8(i) = a0 + u*(a1 + u*(a2 + u*(a3 + a4*u)));
  p6(i) = a0 + u*(a1 + u*(a2 + u*a3));
end