Solution (source code)

= Solution

The canonical <Poisson bracket> is
$$
\{f,g\}=\sum_k\left(\frac{\partial f}{\partial q_k}\frac{\partial g}{\partial p_k}-\frac{\partial f}{\partial p_k}\frac{\partial g}{\partial q_k}\right).
$$
The coordinate derivatives are <Kronecker deltas>, so direct substitution gives
$$
\boxed{\{q_i,q_j\}=0,\qquad\{p_i,p_j\}=0,\qquad\{q_i,p_j\}=\delta_{ij}.}
$$