code: def Qk(self, k):
          
          return self.Q[:, :, self.index[self.Q_time_var_index, k]]
docs: function (k). Returns noise matrix of dynamic model on iteration k.
                  k (iteration number). starts at 0