%%need to finish

function[bound] = inferHMM2(cur,n)

global DST;

bound = 0;

obslik = eval_pdf_cond_multinomial(DST(cur).y{n}', DST(cur).eta');

if (cur==1)
    %need to fix xi
    [alpha, beta, DST(cur).s{n}, ll, DST(cur).xi{n}, gamma2]= fwdback( DST(cur).phi, DST(cur).Phi', obslik );

    bound = bound - computeH2(cur,n);

    bound = bound + entropy4(cur,n);
    
else
    
    %need to fix
    inferChain(cur);
    
    bound = bound + entropy4(cur,n);
 
end

