Fs = 44100;
p = 200; %length of delay line
tmax = 3;

sr = rand(1,p); %initial shift register contents
[b,a] = butter(1,0.005); %modify pluck -- lower is softer
sr = filter(b,a,sr);

out = zeros(1,Fs*tmax);
ptrout = 2; %pointer to sr
ptrin = 1; %pointer to sr
factor = 0.495; %decay factor 0<factor<=0.5

for i=1:Fs*tmax
    out(i) = sr(ptrin);
    sr(ptrin) = (sr(ptrin)+sr(ptrout))*factor;
    %updata and wrap pointers
    if (ptrin==p) ptrin=1;
    else ptrin=ptrin+1;
    end
    if (ptrout==p)ptrout=1;
    else ptrout=ptrout+1;
    end
end

sound(out,Fs)