function [theta]=train_linear_reg(X,y_train,lambda)
%number of columns
m=size(X,2);
theta=pinv(X'*X+lambda*eye(m))*X'*y_train;
end