I was hoping to fit a cylinder to some data points in 3d space and find the equation of the long axis of the cylinder.

Sample data (actual data much larger matrix):

I am trying to write a procedure that uses LinearAlgebra[SingularValues]. For some reason, the line U,S,Vt := LinearAlgebra[SingularValues](covariance, output=['U','S','Vt']); causes the procedure to give the error message: Error, (in LinearAlgebra:-LA_Main:-SingularValues) invalid output (U) specified However, outside the procedure the same line works perfectly. Any help will be much appreciated!