function Y=ber(X,n,f)
%BER Approximation de Bernstein sur [0,1]
% Y=ber(X,n,f)
% Renvoie l'approximation de Bernstein y=Pn(f)(X) de f(X).
% X est un vecteur quelconque a valeurs dans [0,1]
% Fait appel a StixBox pour les coef. binomiaux
%
% Exemple : X=[.1:.1:.9]; f=inline('x.^4.*sin(x)'); ber(X,100,f)
%
 m=length(X);
 D=repmat([0:n],m,1);  % Réplication de D
 B=bincoef(D,n);         % Calcul des coefficients binomiaux
 F=repmat(feval(f,[0,[1:n]/n]),m,1); % Calcul et réplication de f
 X=reshape(X,1,m);   % Remise en forme de X 
 XX=repmat(X',1,n+1); % Réplication de X
 P=(XX.^D) .* ((1-XX).^(n-D)); % Calcul des coefficients en X
 Y=sum((B.*F.*P)');  % Calcul de la valeur du polynôme