X = [1.6 1.4 1.2 1 0.8]; % in terms of r
Y = [0.3 0.5 0.7 1]; % in terms of r
Z = [0.093 0.072 0.053 0.027 0.0005; ...
0.165 0.097 0.06 0.058 0.021 ; ...
0.226 0.134 0.11 0.089 0.043 ; ...
0.301 0.205 0.163 0.146 0.104 ];
surf(X,Y,Z)
xlabel('X/r')
ylabel('Y/r')
zlabel('Z')