Extend linear regression to multiple features using matrix form ŷ = Xw + b and vectorized NumPy operations.