function [nlp] = nlp_linear_exp(x, xi, yi)
a = x(1);
b = x(2);
s = x(3);
nlp = numel(yi)*log(s) + sum(abs(yi - a*xi - b))/s;
