We investigate a biaxial Bianchi IX model with positive cosmological constant, which is sometimes called the Lambda-Taub-NUT spacetime, whose exact solution is well known. The minisuperspace of biaxial Bianchi IX models admits two non-trivial Killing tensors that play an important role for deriving the Taub-NUT metric. We also give a brief discussion about the asymptotic behaviour of Bianchi IX models.