[fd,SST,Sheetnames,Sheetpos] = xls_open('E:\LI_messungen_alle.xls')
Sheets = readxls('E:\LI_messungen_alle.xls')
typeof(Sheets)
s1=Sheets(1)
data=s1.value
x=data(1,2:8)
y=data(2:8,1)
z=data(2:8,2:8)
xx = linspace(0,200,200);
yy = linspace(0,150,200);
[xp,yp] = ndgrid(xx,yy);
zp = linear_interpn(xp,yp, x, y, z);
zm=min(z)
zM=max(z)
clf()
surf(xx,yy,zp)
f=get("current_figure");
f.color_map=jetcolormap(512);
h=get("hdl") //get handle on current entity (here the surface) 
h.color_flag=1;
colorbar(zm,zM)


scf(1);
subplot(221);
xx = linspace(0,200,400);
yy = linspace(0,150,400);
[xp,yp] = ndgrid(xx,yy);
zp = linear_interpn(xp,yp, x, y, z);
surf(xx, yy, zp);
f=get("current_figure");
f.color_map=jetcolormap(512);

scf(1);
subplot(222);
x=data(1,2:8);
y=data(2:8,1);
z=data(9:15,2:8)
xx = linspace(0,200,400);
yy = linspace(0,150,400);
[xp,yp] = ndgrid(xx,yy);
zp = linear_interpn(xp,yp, x, y, z);
zm = min(zp);
zM = max(zp);
colorbar(zm,zM);
surf(xx, yy, zp);