function [H,H_naive]=temporal_entropies(sequenz,m1,m2,nbin);
%Als input werden zwei symbol sequenzen gleicher Länge benötigt,
%als Spalten in sequenz
%m1 bezeichnet die Vergangenheit die in Sequenz 1 berücksichtigt wird
%m2 die für Sequenz zwei. 
%Die Sequnez enthält die Symbole 0,1,...,nbin-1
ngesamt=length(sequenz);
vec1=nbin.^(m1:-1:0);
vec2=nbin.^(m2:-1:0);
vec3=nbin.^((m1+m2):-1:0);
m=max(m1,m2);
nm1=zeros(1,2^(m1+1));
nm2=zeros(1,2^(m2+1));
nm3=zeros(1,2^(m1+m2+1));
nm4=zeros(1,2^(m1+m2+1));
for l=1:(ngesamt-m)
     indx1=vec1*sequenz(l:(l+m1),1)+1;
     nm1(indx1)=nm1(indx1)+1;
     indx2=vec2*sequenz(l:(l+m2),2)+1;
     nm2(indx2)=nm2(indx2)+1;   
     indx3=vec3*[sequenz(l:(l+m2-1),2);sequenz(l:(l+m1),1)]+1;
     nm3(indx3)=nm3(indx3)+1;
     indx4=vec3*[sequenz(l:(l+m1-1),1);sequenz(l:(l+m2),2)]+1;
     nm4(indx4)=nm4(indx4)+1;
 end
 nm10=nm1(find(nm1>0));
 ngesamt=sum(nm10);
 H11=log(ngesamt)-sum(nm10*psi(nm10)')/(ngesamt);
 H_naive11=log(ngesamt)-sum(nm10*log(nm10)')/(ngesamt);
 
 nm1=nm1(1,2*(1:2^m1))+nm1(1,2*(1:2^m1)-1); 
 
 nm10=nm1(find(nm1>0));
 ngesamt=sum(nm10);
 H12=log(ngesamt)-sum(nm10*psi(nm10)')/(ngesamt);
 H_naive12=log(ngesamt)-sum(nm10*log(nm10)')/(ngesamt);
 
 nm20=nm2(find(nm2>0));
 H21=log(ngesamt)-sum(nm20*psi(nm20)')/(ngesamt);
 H_naive21=log(ngesamt)-sum(nm20*log(nm20)')/(ngesamt);
 
 nm2=nm2(1,2*(1:2^m2))+nm2(1,2*(1:2^m2)-1); 
 
 nm20=nm2(find(nm2>0));
 H22=log(ngesamt)-sum(nm20*psi(nm20)')/(ngesamt);
 H_naive22=log(ngesamt)-sum(nm20*log(nm20)')/(ngesamt);
 
 nm30=nm3(find(nm3>0));
 H31=log(ngesamt)-sum(nm30*psi(nm30)')/(ngesamt);
 H_naive31=log(ngesamt)-sum(nm30*log(nm30)')/(ngesamt);
 
 nm3=nm3(1,2*(1:2^(m1+m2)))+nm3(1,2*(1:2^(m1+m2))-1); 
 
 nm30=nm3(find(nm3>0));
 H32=log(ngesamt)-sum(nm30*psi(nm30)')/(ngesamt);
 H_naive32=log(ngesamt)-sum(nm30*log(nm30)')/(ngesamt);
 
 nm40=nm4(find(nm4>0));
 H41=log(ngesamt)-sum(nm40*psi(nm40)')/(ngesamt);
 H_naive41=log(ngesamt)-sum(nm40*log(nm40)')/(ngesamt);
 
 nm4=nm4(1,2*(1:2^(m1+m2)))+nm4(1,2*(1:2^(m1+m2))-1); 
 
 nm40=nm4(find(nm4>0));
 H42=log(ngesamt)-sum(nm40*psi(nm40)')/(ngesamt);
 H_naive42=log(ngesamt)-sum(nm40*log(nm40)')/(ngesamt);
 
 H(1)=H11-H12-(H31-H32);
 H(2)=H21-H22-(H41-H42);
 
 H_naive(1)=H_naive11-H_naive12-(H_naive31-H_naive32);
 H_naive(2)=H_naive21-H_naive22-(H_naive41-H_naive42);