function response = Hill( par, dose )
%% Return the predicted response for a Hill function
ymax = par(1);
EC50 = 10^par(2);
n = par(3);

response = ymax./(1 + (EC50./dose).^n);
end