function c = linL(theta);

global Winv g addobs population vars explinspec P

explinspec = exp(vars * theta');    % calculate probabilities according
P = explinspec./(1 + explinspec);   % to equation (15)

g = population - addobs * (1./P);   % g is phi(j) in equation (6)
c = g' * Winv * g;                  % objective of minimization problem (7)