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