function w = train(t,y)
% w = estimate(t,y)
% t -> the training set
% y -> the empirical values
% w == estimated weights (LSE)
t1s = [ones(size(t,1),1) t];
w = pinv(t1s) * y;
