function tot = sqsum(u,x,mu)

%used in err-fit.m

tot = u'*((x - mu).^2) ;