Построение 3D графика

clear ,clc
 
M = [25, 9, 3.25, 206;
    8, 3, 3.25, 64;
    9, 50, 3.25, 62;
    12, 80, 3.25, 80;
    8, 8, 6.5, 75;
    12, 20, 6.5, 209;
    5, 42, 6.5, 26;
    3, 65, 6.5, 68
    14, 14, 9.75, 106
    12, 38, 9.75, 32
    3, 46, 9.75, 18
    12, 86, 9.75, 37];
 
x = M(:,1);
y = M(:,2);
z = M(:,3);
a = M(:,4);
[X Y Z] = ndgrid(linspace(min(x),max(x),20), linspace(min(y),max(y),20), unique(z)); % создаем трехмерную сетку
A = griddata(x,y,z,a, X,Y,Z,'cubic'); % интерполируем значения концентрации
slice(X,Y,Z,A,[],[],uinque(z), 'facealpha',0.5) % отображаем каждый слой по высоте
shading interp
colorbar

clear ,clc
 
M = [25, 9, 3.25, 206;
    8, 3, 3.25, 64;
    9, 50, 3.25, 62;
    12, 80, 3.25, 80;
    8, 8, 6.5, 75;
    12, 20, 6.5, 209;
    5, 42, 6.5, 26;
    3, 65, 6.5, 68
    14, 14, 9.75, 106
    12, 38, 9.75, 32
    3, 46, 9.75, 18
    12, 86, 9.75, 37];
 
x = M(:,1);
y = M(:,2);
z = M(:,3);
zu = unique(z);
a = M(:,4);
[Xm Ym] = meshgrid(linspace(min(x),max(x),20), linspace(min(y),max(y),20)); % создаем трехмерную сетку
A = zeros(size(Xm,1), size(Xm,2), length(zu));
X = A; Y = A; Z = A;
for i = 1:length(zu)
  A(:,:,i) = griddata(x(z==zu(i)), y(z==zu(i)), a(z==zu(i)), Xm, Ym);
  Z(:,:,i) = zu(i)*ones(size(Xm));
  X(:,:,i) = Xm;
  Y(:,:,i) = Ym;
end
slice(X,Y,Z,A,[],[],zu, 'facealpha',0.5) % отображаем каждый слой по высоте
shading interp
colorbar