Matlab Code For Lagrange Fractional Delay

function h = lagrange(N, delay)
%LAGRANGE  h=lagrange(N,delay) returns order N FIR 
%          filter h which implements given delay 
%          (in samples).  For best results, 
%          delay should be near N/2 +/- 1.
n = 0:N;
h = ones(1,N+1);
for k = 0:N
    index = find(n ~= k);
    h(index) = h(index) *  (delay-k)./ (n(index)-k);
end