  We present two kinds of Bochkov-Kuzovlev work equalities in a two-level
system that is described by a quantum Markovian master equation. One is based
on multiple time correlation functions and the other is based on the quantum
trajectory viewpoint. We show that these two equalities are indeed equivalent.
Importantly, this equivalence provides us a way to calculate the probability
density function of the quantum work by solving the evolution equation for its
characteristic function. We use a numerical model to verify these results.
