clear; clf; format compact
c=4;

for nitems = 1:12;

    retint=2;

    clear dtim;

    nlists=1;

    if nitems==1; dtim=[0]; end
    if nitems==2; dtim=[1 0]; end
    if nitems==3; dtim=[2 1 0]; end
    if nitems==4; dtim=[3 2 1 0]; end
    if nitems==5; dtim=[4+120 3+120 2+120 1+120 0]; end
    if nitems==6; dtim=[5+120 4+120 3+120 2+120 1 0]; end
    if nitems==7; dtim=[6+120 5+120 4+120 3+120 2 1 0]; end
    if nitems==8; dtim=[7+120 6+120 5+120 4+120 3 2 1 0]; end
    if nitems==9; dtim=[8+240 7+240 6+240 5+240 4+120 3+120 2+120 1+120 0]; end
    if nitems==10; dtim=[9+240 8+240 7+240 6+240 5+120 4+120 3+120 2+120 1 0]; end
    if nitems==11; dtim=[10+240 9+240 8+240 7+240 6+120 5+120 4+120 3+120 2 1 0]; end
    if nitems==12; dtim=[11+240 10+240 9+240 8+240 7+120 6+120 5+120 4+120  3 2 1 0]; end

    dtim=dtim+retint;

    for n=1:nitems*nlists; dtim(n)=log(dtim(n)); end

    for n=1:nitems*nlists; dd(n)=0.0; end

    for n=1:nitems*nlists;
        for m=1:nitems*nlists;
            distance=abs(dtim(n)-dtim(m));
            dd(n) = dd(n) + exp(-c*(distance));
        end
        serpos(n)=1/dd(n);
    end

    results(nitems)=serpos(nitems*nlists)^3;

    plot(results)
    hold on

end
