The ‘y’ expression is not a function of ‘x’.
Correcting that (and assuming the independent variable is actually ‘ec’) —
ec = linspace(0, 0.035, 50);
Mt = rand; % Missing Value
h = rand; % Missing Value
y=Mt*(1./ec).*(1./sqrt(1-((ec-h)./ec).^2));
figure
plot(ec, real(y), 'DisplayName','Real')
hold on
plot(ec, imag(y), 'DisplayName','Imaginary')
hold off
ylim([-1 1])
legend('Location','best')
There were also some missing element-wise operators.
.
