<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def major_axis(points):
 """ Returns an approximate vector representing the major axis of points Parameters points: (n, dimension) float, points in space Returns axis: (dimension,) float, vector along approximate major axis """

    U, S, V = np.linalg.svd(points)
    axis = util.unitize(np.dot(S, V))
    return axis