## Codiguinho que calcula os coeficientes da
## aproximacao polinomial de grau k de 
## f^{(m)} (a) = m! \sum_{i=1}^n c(m-1,i) y_i
## sendo y_i a amostra tomada em x_i=-nm+i-1,
## para todo m <= k.
## No comeco se subtrai "a" de todos os x para
## simplificar a conta. 
## Notar que k deve ser <= np+nm
nm=4;np=0;k=4;a=0.;  ## dados
n=nm+np+1;
x=[-nm:1:np]-a;
M=ones(n,1);
for i=1:k
	M=[M,x'.^i];
endfor
c=M\eye(n)
## agora calcula as constantes de erro das derivadas
pol1=x.^(k+1)./factorial(k+1);
pol2=x.^(k+2)./factorial(k+2);
for j=1:k+1
## ruido
gamma(1,j)=norm(c(j,:));
## erro de consistencia -- derivada j-1
cerr(1,j)=-c(j,:)*pol1';
cerr(2,j)=-c(j,:)*pol2';
endfor
gamma
cerr