$${dF \sim r^{–2} I_1 I_2 dl_1 dl_2 (3 cos \phi_1 cos \phi_2 – 2 cos \varepsilon)}$$