vndinterp( fun, x, nbplot)

function vndinterp( fun, x, nbplot)
%VNDINTERP   trac de l'interpole polynmiale de la fonction fun
%            vndinterp( fun, x {, nbplot})
%            fun est le nom de la fonction  interpoler.
%            x est la base d'interpolation
%            nbplot est le nombre (facultatif) de points de trac
%            qui seront repartis rgulirement entre min(x) et ma(x)
%            si nbplot est absent il vaut par dfaut 10*length(x).

x = x(:);                           % on s'assure que x est un vecteur
sort(x);                            % et qu'il est class
n = length(x);
if nargin <= 2                      % combien de points faut-il tracer
   nbplot = 10*n;                   % si ce n'est pas prcis dans l'appel ?
end;

y = feval(fun,x);                   % calcul des ordonnes
                                    % feval value fun en tout lement de x
v = vander(x);                      % cration de la matrice de vandermonde

p = v\y;                            % calcul du polynme d'interpolation
x1 = linspace(x(1), x(n), nbplot);  % sur-chantillonnage
z = polyval(p',x1);
% traces
plot(x,y, 'go', x1, feval(fun,x1), 'r', x1, z, 'y');
title(['\bf Interpolation de degr ',int2str(n-1),' de ', fun, ...
      ' sur [',num2str(x(1)), ', ' num2str(x(n)), ']']);