function H=temporal_entropies(feld,mmax);
%Zeilen in Feld Zeit
%Spalten in Feld Raum
%Feld soll zun?chst nur bin?re Variablen enthalten
[n1,n2]=size(feld);
ngesamt=n1*n2;
s1=sum(sum(feld));
n0=[s1,n1*n2-s1];
H=zeros(mmax+1,1);
if (s1*(ngesamt-s1) >0)    
    H(1)=log(ngesamt)-sum(n0*psi(n0)')/(ngesamt);
    m=mmax;
    nm=zeros(1,2^(m+1));
    vec2=2.^(0:m);
    for k=1:n2
        for l=1:(n1-m)
            indx=vec2*feld(l:(l+m),k)+1;
            nm(indx)=nm(indx)+1;
        end
    end    
    nm0=nm(find(nm>0));
    ngesamt=sum(nm0);
    H(m+1)=log(ngesamt)-sum(nm0*psi(nm0)')/(ngesamt);
    m=m-1;
    while (m >= 1)
        nm=nm(1,2*(1:2^(m+1)))+nm(1,2*(1:2^(m+1))-1);       
        nm0=nm(find(nm>0));
        H(m+1)=log(ngesamt)-sum(nm0*psi(nm0)')/(ngesamt);
        m=m-1;
    end
else 
    H(1:mmax+1)=0;
end