% Poisson kernel on the unit disc

x=0;
y=1;

r = [0:0.05:1];
theta = [0:2*pi/100:2*pi];

X = r'*cos(theta);
Y = r'*sin(theta);

n = 1 - X.^2 - Y.^2;
D = 2*pi*((X-x).^2 + (Y-y).^2);
z = n./D;

colormap(hsv);
surf(X,Y,z);
shading flat; 


