% Tomimatsu-Sato solution, delta = 3
@ p, q, x, y;
a = p*(x^2-1)^3*(x^3+3*x)+1i*q*(1-y^2)^3*(y^3+3*y)-p*q^2*
    (x^2-y^2)^3*(x^3+3*x*y^2)
    -1i*p^2*q*(x^2-y^2)^3*(y^3+3*x^2*y);
b = p^2*(x^2-1)^3*(3*x^2+1)-q^2*(1-y^2)^3*(3*y^2+1)
    -12i*p*q*x*y*(x^2-y^2)*(x^2-1)*(1-y^2);
u = a*a'-b*b';
v = x^2-1;
w = y^2-1;
t = (u*(a:x:x*b-a*b:x:x)*v+(2*x*u-2*v*(a'*a:x-b'*b:x))
    *(a:x*b-a*b:x))-(u*(a:y:y*b-a*b:y:y)*w
    +(2*y*u-2*w*(a'*a:y-b'*b:y))*(a:y*b-a*b:y));
? t;
? t!;
