piva

Просто добавить точки перед операторами деления, чтобы сделать поэлементное деление и получить матрицу Z, размер которой будет такой же как у матриц X и Y.

[X,Y]= meshgrid(1:.1:6, -5:.1:5);
Z= (1./X.^3 + 1./Y.^3)
surf(X,Y,Z); grid on; colormap;

Если есть координаты границ объекта и предположить что лицо - плоскость, то расстояние между ними можно вычислить просто как расстояние между двумя точками.

Если уже хотите учесть что лицо не плоское, то можете нарисовать на бумаге эллипс и посчитать расстояние между объектами, лежащими на его границе. И потом уже использовать эти расчёты чтобы определить расстояние внутри своей программы.

Что нужно для этого знать: пропорции головы (чтобы определить расстояние) + некий объект на фотографии размер которого Вам известен (который будет использоваться как коэффициент масшабирования для определения реального размера). Например, знаете что размер зрачка (15 мм). Расстояние между глаз - три таких зрачка. Т.е. 45 мм.

Просто потеряли знак точки и не получилось поэлементное деление.

epsilon=0.0382;
f0=451.36;

f=10:0.01:100;
eta=f/f0;
k_din=(sqrt(1-((epsilon*eta).^2)))./(sqrt(((1-(eta.^2)).^2)+((epsilon*eta).^2)));
plot(f,k_din);
grid on;
axis([10,200,0.98,1.05]);