anchor.diff:  		((>= k (get-states hstate)))
   	      (progn
   		(setf value (* (aref delta (- tm 1) k)
  -			       (aref (get-transition-matrix hstate) k j)))
  +			       (aref (get-transition hstate) k j)))
   		(if (> value max-value)
   		    (setf max-value value
   			  max-value-indice k))))
  -	    (setf (aref delta tm j) (* max-value (aref (get-measure-matrix hstate) j (aref o tm)))
  +	    (setf (aref delta tm j) (* max-value (aref (get-measure hstate) j (aref o tm)))
   		  (aref psi tm j) max-value-indice))))
         (setf (aref q (- tms 1)) 1)
         (do ((l 0 (1+ l)))
  
positive.code:     (progn
        (do ((i 0 (1+ i)))
  	  ((>= i (get-states hstate)))
  	(setf (aref alpha 0 i) (* (aref (get-initial-distributions hstate) i)
  				  (aref (get-measure-matrix hstate) i (aref o 0)))))
        (do ((tm 0 (1+ tm)))
  	  ((>= tm (- tms 1)))
  	(do ((j 0 (1+ j)))
  
negative.code: 	    (progn
  	      (setf (aref xi tm i j) (* (aref alpha tm i)
  					(aref beta (+ tm 1) j)
  					(aref (get-transition hstate) i j)
  					(aref (get-measure hstate) j (aref o (+ tm 1)))))
  	      (incf sum (aref xi tm i j)))))
  	(do ((k 0 (1+ k)))
  	    ((>= k (get-states hstate)))
  