---

X = [0, 300, 650, 1000, 1300, 1700] % коорд. х
H1 = [1223, 1230, 1236, 1228, 1223, 1221] % пласт 1 кровля
H2 = [1229, 1235, 1241, 1232, 1226, 1224] % граница гл.пласт
H3 = [1235, 1239, 1246, 1235, 1229, 1228] % пласт 2 подошв

xx = 0:10:1700;
H12 = spline(X, H1, xx);
H22 = spline(X, H2, xx);
H32 = spline(X, H3, xx);
hold on

fill([xx, xx(end:-1:1)], [H12, H22(end:-1:1)], 'green')
fill([xx, xx(end:-1:1)], [H22, H32(end:-1:1)], 'blue')

plot(X, H1, 'o')
plot(X, H2, 'o')
plot(X, H3, 'o')
set(gca, 'YDir', 'reverse')