Abstract: We solve the linear GinzburgLandau (GL) equation in the presence of a uniform magnetic field with cylindrical symmetry and we find analytic expressions for the eigenfunctions in terms of the confluent hypergeometric functions. The discrete spectrum results from an implicit equation associated to the boundary conditions and it is resolved in analytic form using the continued fractions formalism. We study the dependence of the spectrum and the eigenfunctions on the sample size and the surface conditions for solid and hollow cylindrical superconductors. Finally, the solutions of the nonlinear GL formalism are constructed as expansions in the linear GL eigenfunction basis and selected by minimization of the free energy. We present examples of vortex states and their energies for different samples in enhancing/suppressing superconductivity surroundings.