function [U,D,V]=prolate(T,tau,Omega)

%              [U,D,V]=prolate(T,tau,Omega);
%
%Calcule le spectre de l'opérateur Phi(Omega,T) lorsque [-T,T] est 
%échantillonné avec un pas tau. Les valeurs propres lambda (positives) sont
%en ordre décroissant sur la diagonale de D. Les colonnes de U sont les
%vecteurs propres (normalisés) correspondants. Pour des raisons pratiques,
%2T/tau est choisi inférieur à 256.
 
t=-T:tau:T;
l=length(t); 
A=ones(l,l);
for k1=1:l,
  for k2=1:l,
    A(k1,k2)=tau*sincnorm(2*(k1-k2)*tau*Omega/pi)*(Omega/pi);
  end
end
[U,D,V]=svd(A);