function y=cwf(x)
t=x;
N=max(10,floor(abs(6*x)));
p=zeros(N+1,1);

g=ones(N,1);
i=1:N;
f(1)=1;
for j=2:N
   f(j)=f(j-1)*j;
   g(N-j+1)=f(j);
end;

pp=1./(N-i+1)./g';
p(1:N)=pp;
y=polyval(p,t);

