The Hubbard model beyond the two-pole approximation: A composite operator method study