Freezing on "solve" and "vpasolve" when using more than one equation

Matlab freezes (won't stop running but doesn't print anything) on the last line

a1 = 26;
a2 = 1;
b1 = 3;
b2 = 1;

syms x y

p = 4*(a1+a2+b1+b2);
functionF = (p./((3*x+a1).^2 + (3*y+b1).^2 + (a1+b1))) - (2*p./((3*x-a2).^2 + (3*y+b2).^2 + (a2+b2))) - (3*p./((3*x+a1).^2 + (3*y-b2).^2 + (a1+b2))) + (4*p./((3*x-a2).^2 + (3*y-b1).^2 + (a2+b1)));
domain = [-10,10,-10,10];
partialX = diff(functionF, x)
partialY = diff(functionF, y)
sol= solve([partialX,partialY])

Any help would be v appreciated

[Star Strider]

Adding these assignments:

[nX,dX] = numden(partialX)
[nY,dY] = numden(partialY)

reveals that the numerators are 12-th degree polynomials. An analytic solution may not be possible.