% Prescribe arrays
K = ...;
X = ...;
L = ...;
Y = ...;
% Define function
fun = @(a,b,c,d,h) K.*log(a./(a-X)) + K/b.*log(X/c) + (d+L*h).*X - Y;
% Define initial values for a,b,c,d,h
a0 = ...;
b0 = ...;
c0 = ...;
d0 = ...;
h0 = ...;
% Call the optimizer
p = lsqnonlin(@(p)fun(p(1),p(2),p(3),p(4),p(5)),[a0,b0,c0,d0,h0]);
% Output fit parameters
a = p(1)
b = p(2)
c = p(3)
d = p(4)
h = p(5)
