Интересно, и я думал, что решение задачи как-то связано с нарушением исходной симметрии при обходе вокруг квадруполя. Теперь надо думать, как вычислить потенциал, по принципу суперпозиции он будет равен

и теперь надо вычислить эти расстояния так, чтобы потенциал трехполья выразился через сферические координаты.
Впрочем, если посчитать его квадрупольный момент (полагая четвертый заряд нулевым) и вычислить компоненты его квадрупольного поля, оно будет зависеть от двойных углов:



Таким образом, зависимость от двойного угла и здесь не исключается, - но, может, они убывают быстрее чем те, которые надо вычислить, и приближенно такая система подойдет.