x = 2.6:0.1:8; y = -3.6:0.1:1.8;
[x ,y] = meshgrid(x, y);
z = (conj(x).*(21 - 10i) - conj(y).*(10 + 21i) + (conj(x) - conj(y)*1i).^2.*(4 - 18i) + (conj(x) - conj(y)*1i).^3.*(10 - 10i) + (17 - 12i)).*(x*(21 + 10i) - y*(10 - 21i) + (x + y*1i).^2.*(4 + 18i) + (x + y*1i).^3.*(10 + 10i) + 17 + 12i) + (max(0, (x-5.3).^2+(y+0.9).^2-2.7^2)).^2 + (max(0, -y-0.9)).^2;

surf(x, y, z)
shading interp % эта команда сглаживает поверхность