How to numerically minimize sum of squared errors in regression analysis in Matrix form

Question

So i need to do a matrix calculation (which I did already, heres the code)

> incomeX=c(35,90,47,45,68,70,55)
> engelY=c(21,8,16,14,11,12,10)
> Y=enge