function [SSE, SSR, b] = modelo(x, y, campos)

K = length(x(1, :));

X = ones(length(x(:,1)),1);

for m = 1: length(campos)
   X = [X, x(:,campos(m))];
end;


b = inv(X'*X)*X'*y;

SSE = (X*b-y)'*(X*b-y);

yp = mean(y);

SSR = (X*b-yp)'*(X*b-yp);

