Assuming[ρ > 0, Integrate[(ρ Exp[-(ρ x)]) (ρ Exp[-(ρ y)]), {x, 0, Infinity}, {y, 0, x}]]