function res=expression_part_symb(c,n)


%	expression_part_symb : calcul de l'intégrale de x^n*w(x) sur (a,b) (en symbolique)
%
% *********************************************************
%
%	expression_part_symb(c,n) : calcul en symbolique de \int_a^b x^n w(x) dx
%        
%
%	variables d'entrées : 
%    *  n : entier naturel.
%    * c définit le poids w : 
%            c=1 : Gauss-Legendre,
%            c=2 : Gauss-Tchebytchev,
%            c=3 : Gauss-Hermite,
%            c=4 : Gauss-Laguerre.
%	variables de sortie
%   *res=valeur de l'inétrale
% 
% 
% ************ 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 
if ~((c==1) | (c==2) | (c==3) | (c==4))
   error('c doit être un entier égal à 1, 2, 3 ou 4');
end

% Corps d'algorithme.
syms x a b;
switch c
   case 1
      a=sym(-1);
      b=sym(1);
      w=sym(1);      
   case 2
      a=sym(-1);
      b=sym(1);
      w=1/sqrt(1-x^2);      
   case 3
      a=-sym(inf);
      b=sym(inf);
      w=exp(-x^2);
   case 4
      a=sym(0);
      b=sym(+inf);
      w=exp(-x);
end     
res=simplify(int(x^n*w,a,b));
