clear all;
clc;

A=[500,10,53.6          
   700,10,59.7
   900,10,63.5
   1100,10,66.2
   1300,10,67.2
   1500,10,70.6
   1550,10,70
   1600,10,69.4
   1700,10,69.3
   1900,10,69.8
   2100,10,71.2
   2300,10,67.8
   2450,10,68.2
   2500,10,67.1
   2550,10,68.2
   2700,10,67.1
   2900,10,67.1
   3100,10,66.8
   3300,10,66.2
   3500,10,67.9
   3700,10,69
   3900,10,70.2
   4000,10,71.8
   4050,10,72.8
   4100,10,74.4
   4300,10,74.6
   4500,10,74.5
   4700,10,74.1
   4900,10,75.7
   5100,10,75.3
   5300,10,75.3
   5500,10,74.2
   5700,10,74.4
   5400,50,75.6
   5300,60,75.8
   500,70,57.5          
   700,70,60.7
   900,70,64.9
   1100,70,68.1
   1300,70,69.8
   1500,70,72.4
   1550,70,71.7
   1600,70,73
   1700,70,72.3
   1900,70,73.1
   2100,70,70.6
   2300,70,71.4
   2450,70,73.2
   2500,70,72.5
   2550,70,70.1
   2700,70,69.3
   2900,70,68.6
   3100,70,67
   3300,70,67.2
   3500,70,68.3
   3700,70,71.2
   3900,70,72
   4000,70,71.9
   4050,70,73.7
   4100,70,73.4
   4300,70,76.5
   4500,70,75.1
   4700,70,75.2
   4900,70,76.6
   4550,80,75.7   
   4300,90,74.6
   4100,100,74.4
   500,110,56
   700,110,59.1
   900,110,64.4
   1100,110,67.8
   1300,110,67.2
   1500,110,70.3
   1550,110,69.8
   1600,110,70.5
   1650,110,70.9
   1700,110,70.2
   1900,110,70.5
   2100,110,70.3
   2300,110,69.4
   2450,110,70.4
   2500,110,70.2
   2550,110,69.1
   2700,110,70
   2900,110,67.4
   3100,110,67.9
   3300,110,67
   3500,110,68.8
   3700,110,70.2
   3900,110,70.9
   3700,120,71.5      
   3600,130,69.3
   3450,140,69.6
   500,150,57.9
   700,150,60.3
   900,150,65
   1100,150,67.3
   1300,150,69.7
   1500,150,70.6
   1550,150,70.3
   1600,150,71.5
   1650,150,70.8
   1700,150,70.4
   1900,150,69.8
   2100,150,70
   2300,150,70
   2500,150,70.8
   2700,150,69
   2900,150,66.8
   3100,150,67.3
   3300,150,68.4
   3200,160,69.3   
   3100,170,69
   3000,180,69.4
   500,190,58.5
   700,190,61
   900,190,65.3
   1100,190,67
   1300,190,70.3
   1500,190,70.1
   1550,190,70.1
   1600,190,70.9
   1650,190,71.1
   1700,190,71.2
   1900,190,70.2
   2100,190,71.4
   2300,190,71.3
   2450,190,70.9
   2500,190,71.6
   2700,190,70.5
   2900,190,68.5 
   2800,200,69.3
   2750,210,70.8
   500,220,59
   700,220,61.3
   900,220,65
   1100,220,67
   1300,220,70.4
   1500,220,70.5
   1550,220,70.1
   1600,220,71.7
   1650,220,71.3
   1700,220,71.6
   1900,220,70.9
   2100,220,71.8
   2300,220,71.6
   2450,220,71.9
   2500,220,72.3
   2550,220,71.3
   2700,220,69.7];

x=A(:,1);
y=A(:,2);
z=A(:,3);
stem3(x,y,z)
hold on
tx= min(x):1:max(x);
ty= min(y):1:max(y);
[XI,YI]=meshgrid(tx,ty);
ZI=griddata(x,y,z,XI,YI,'cubic');
surf(XI,YI,ZI);
hold off
grid on
set(gca,'CameraPosition',[70 70 100])
shading interp
colorbar EastOutside
colormap ('cool')
xlim([min(x) max(x)])
ylim([min(y) max(y)])
zlim([0 max(z)])