  We compute the Bellman function of three integral variables associated to the
dyadic maximal operator on a subset of its domain. Additionally, we provide an
upper bound for the whole domain of its definition.
