We obtain bounds for the density of \(\left (x_{1,n}+W_{t}, x_{n+1}+{{\int }_{0}^{t}} |x_{1,n}+\right .\\\left .W_{s}|^{k}ds\right )\) where (Wt)t≥0 is a standard Brownian motion of ℝn, k∈ℕ∗ is even and x=(x1,n,xn+1)∈ℝn+1. This process satisfies a weak Hormander condition but the support of its density is not the whole space. Also, the Density has various asymptotic regimes depending on the starting/final points considered (which are as well related to the number of brackets needed to span the space in Hormander’s theorem). The proofs of lower and upper bounds are based on Harnack inequalities and Malliavin calculus respectively. The case of the joint law of Brownian motion and the integral of odd powers of its coordinates is also considered.