
% okay, this isn't the most elegant program, but it'll do...
a0 = 1;
a1 = 1;
a2 = 1/gamma(2+1);
a3 = 1/gamma(3+1);
a4 = 1/gamma(4+1);
a5 = 1/gamma(5+1);
a6 = 1/gamma(6+1);
a7 = 1/gamma(7+1);
a8 = 1/gamma(8+1);
a9 = 1/gamma(9+1);
a10 = 1/gamma(10+1);

for i=1:21
  u = -1 + (i-1)*1/10;
  p10(i) = a0+u*(a1+u*(a2+u*(a3+u*(a4+u*(a5+u*(a6+u*(a7+u*(a8+u*(a9+u*a10)))))))));
  x(i) = u;
end
