function f = mufunc(x2,theta2w)
% This function computes the non-linear part of the utility (mu_ijt in the %%Guide)

% Written by Aviv Nevo, May 1998.
 

global ns vfull dfull
[n k] = size(x2);
j = size(theta2w,2)-1;
mu = zeros(n,ns);
for i = 1:ns
    	v_i = vfull(:,i:ns:k*ns);
      d_i = dfull(:,i:ns:j*ns);
 		mu(:,i) = (x2.*v_i*theta2w(:,1))+x2.*(d_i*theta2w(:,2:j+1)')*ones(k,1);
end
f = mu;