Двойной коммутатор орбитального момента

Для орбитального момента $\mathbf L=\mathbf r\times\mathbf p$ докажите
$$[L^2,[L^2,r_j]]=2\hbar^2(r_jL^2+L^2r_j).$$
Здесь $r_j$ — любая декартова компонента координаты; действуют канонические соотношения $[r_i,p_j]=i\hbar\delta_{ij}$.