function f=gram(l1,l2);

%             f=gram(l1,l2); 
%
% Programme auxiliaire dans les programmes 
% franklin ou nfrankl. Calcule la matrice de 
% Gram de la base de Franklin pour l'intervalle 
% (1,2^(l2)+1) à l'échelle 2^(l1). 

N=2^(l2-l1)+1;
A=zeros(N,N);
a1=2^(l1)/3;
a2=2^(l1-1)-(2^(l1)/3);
A(1,1)=a1;
A(1,2)=a2;
A(N,N-1)=a2;
A(N,N)=a1;
for k=2:(N-1);
A(k,k-1)=a2;
A(k,k)=2*a1;
A(k,k+1)=a2;
end
f=A;