function res=sommeC(p,n)

%	sommeC : calcul de somme (version C).
%
% *********************************************************
%
%	res=sommeC(p,n)
%	
%  calcul de la somme \sum_{i=1}^n  1/(i^p) en utilisant un calcul 
%  matriciel et en sommant dans l'ordre des indices décroissants.
%
%	variables d'entrées :
%   * p,n : deux entiers 
%
%	variables de sortie
%   * res: résultats de la somme
% 
% 
% ************ Fonctions auxiliaires utilisées ************
%
%       aucune
%
% *********************************************************



% Contrôles d'entrée
% nombre d'arguments
if nargin~=2
   error('nombre d''arguments de la fonction incorrect');
end
% autres tests éventuels 
   

% Corps d'algorithme
res=sum(1./((n:-1:1).^p));
