Skip to content

Instantly share code, notes, and snippets.

@Makazone
Last active November 17, 2015 08:50
Show Gist options
  • Save Makazone/5bd01fa815e7b3f0514c to your computer and use it in GitHub Desktop.
Save Makazone/5bd01fa815e7b3f0514c to your computer and use it in GitHub Desktop.
figure
hold on
% Построим область
for x = -2:0.1:2
for y = -2:0.1:2
z = x + 1i*y;
if (abs(z-(17-8*1i)/13) > 17/11)
plot(x, y, '*r')
newZ = z2(z);
plot(real(newZ), imag(newZ), '*g')
else
end
end
end
% Построим границу
for x = 0:0.01:2*pi
z = (17-8*1i)/13 + 17/11*exp(1i*x);
newZ = z2(z);
plot(real(z), imag(z), '*r')
plot(real(newZ), imag(newZ), '*g')
end
end
function [w7] = z2(z)
w1 = z - (17-8*1i)/13;
w2 = w1 * 17/11;
w3 = 1/w2;
w4 = 1i*(w3+1)/(1-w3);
w5 = w4^(24/43);
w6 = w5*exp(-1i*20*pi/43);
w7 = w6-(8*1i+9)/(16-8*1i);
end
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment