dt = delaunayTriangulation(IM,Z50) ;
t = dt.ConnectivityList ;
p = dt.Points ;
figure (1);
plot3(IM,Z50,MnXdisp,'b*');
hold on;
trisurf(t,p(:,1),p(:,2),MnXdisp')
title('unstructured')
x = IM ; y = Z50 ; z = MnXdisp ;
N = 50 ;
xi = linspace(min(x),max(x),N) ;
yi = linspace(min(y),max(y),N) ;
[X,Y] = meshgrid(xi,yi) ;
Z = griddata(x,y,z,X,Y) ;
figure (2);
plot3(IM,Z50,MnXdisp,'b*');
hold on;
surf(X,Y,Z)
title('structured')