Gradient Descent or Normal Equation?
If you can, it is preferable to use the normal equation to estimate the coefficients for multivariate linear regression. Since the normal equation is a closed-form expression, it will be faster than gradient descent.

Given you have relatively few examples and features, inverting the matrix is not an issue.