The topic of this paper is to describe the 3-D current density in the windings of a 3-D coil, which fills the volume between two coaxial cylinders at a precisely defined distance from each other, and which serves to generate a magnetic field gradient in the center of the cylinder axis. The 3-D current density is considered an unknown input quantity, which is calculated from the known gradient magnetic field output. It is an inverse problem in mathematics, where the direct problems are the calculation of unknown output quantities based on known input quantities. Fourier series expansion methods in the context of cylindrical coordinates were used to describe the 3-D current density. In that case, Bessel functions are used as development components. The current densities, at each point in space, were lined up to represent current lines. Each power line is associated with a coil winding through which a current of a certain strength flows. After that, the principle of discretization of coil windings was applied. Each winding is divided into a large number of elementary segments that were considered as current elements, which create, based on Bio-Savar's law, an elementary magnetic field. In this way, the total, continuous magnetic field is broken into many elementary components, which come from different current elements. An important result of this process is that each current element can be controlled independently by a current source. This means that the output magnetic field of the gradient can be controlled by current sources, which are the input sizes, and this is what is at the core of the topic of this paper.