---

x = [5225, 6660, 8090, 9525]
y = [3375, 3950, 4525, 5100, 5675]
[xg, yg] = meshgrid(x, y) % сетка x, y
z = [0.53, 0.37, 0.32, 0.47, 0.67, 0.54, 0.57, 0.68,
     0.74, 0.66, 0.69, 0.75, 0.65, 0.65, 0.70, 0.72,
     0.42, 0.53, 0.62, 0.63];
figure
mesh(xg, yg, z); % каркасная сетка
xxg = 5225 : 125 : 9525
yyg = 3375 : 125 : 5675
[xx, yy] = meshgrid(xxg, yyg);
figure
vz = interp2(x, y, z, xx, yy, 'spline')
surf(xx, yy, vz) % сплайн кубический
figure
contourf(xx, yg, vz, 10) % контурный график, представление на пл-ти ф-ции 2-х переменных
[c, h] = contourf(xx, yy, vz, 10) % метки значения z(x,y) ф-ции
clabel(c, h)