eb = [0 0.25];
yb = [0.5; 0.75];
K0 = [1.0009 1.3767; 0.9948 1.5583];
Check = [NaN eb; yb K0]
ebq = 0.1; % 'e/b' Value To Interpolate
ybq = 0.6; % 'y/b' Value To Interpolate
K0q = interp2(yb, eb, K0, ybq, ebq) % Interpolated 'K0' Value
figure
surfc(yb, eb, K0, 'FaceAlpha',0.9)
hold on
stem3(ybq, ebq, K0q, 'r', 'filled')
hold off
colormap(turbo)
grid on
xlabel('y/b')
ylabel('e/b')
zlabel('K0 (\theta = 0.20)')
.



