A new method for the characterization of coherent laser beams is proposed. It is based on the non-iterative solution of the transport-of-intensity-equation. The phase to recover is decomposed into paraxial properties of laser beams and a set of lateral shifted radial basis functions, which allows for the derivation of a direct solution of the phase by a least-squares fit without the need of an initial guess. The method is tested with synthetic data to deduce an accuracy metric. Additionally, two real laser beams are characterized. Including the real light source in terms of the reconstructed field allows for a more holistic simulation of optical systems.