Experimental methods based on the magnetic resonance phenomenon belong to the most widely used experimental techniques for investigations of molecular and electronic structure. The difficulty with such experiments, usually a proper interpretation of data obtained from high-resolution spectra, opens new challenges for pure theoretical methods. One of these methods is density functional theory (DFT), that now has an advanced position among a whole variety of computational techniques. This thesis constitutes an effort in this respect, as it presents theory and discusses calculations of electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) parameters of paramagnetic molecules. It is known that the experimental determination of the magnetic resonance parameters of such molecules, especially in the case of NMR, is quite complicated and requires special techniques of spectral detection. On the other hand, paramagnetics play an important role in many areas, such as molecular magnets, active centers in biological systems, and defects in inorganic conductive materials. Therefore, they have spurred great interest among experimentalists, motivating us to facilitate the interpretation of spectral data through theoretical calculations. This thesis describes new methodologies for the determination of magnetic properties of paramagnetic molecules in the framework of DFT, which have been developed in our laboratory, and their applications in calculations of a wide range of molecular systems. The first two papers of this thesis deal with the theoretical determination of NMRparameters, such as nuclear shielding tensors and chemical shifts, in paramagnetic nitroxides that form core units in molecular magnets. The developed methodology is aimed to realize a high calculational accuracy for these systems. The effects of hydrogen bonding are also described in that context. Our theory for the evaluation of nuclear shielding tensors in paramagnetic molecules is consistent up to second order in the fine structure constant and considers orbital, fully anisotropic dipolar, and isotropic contact contributions to the shielding tensor. The next projects concern electron paramagnetic resonance. The well-known EPR parameters, such as the g-tensors and the hyperfine coupling constants are explored. Calculations of electronic g-tensors were carried out in the framework of a spin-restricted open-shell Kohn-Sham method combined with the linear response theory recently developed in our laboratory and allowing us to avoid by definition the spin-contamination problem. The inclusion of solvent effects, described by the polarizable continuum model, extends the possibility to treat molecular systems often investigated in solution. For calculations of the hyperfine coupling constants a so-called restricted-unrestricted approach to account for the spin polarization effect has been developed in the context of DFT. To examine the validity of the approximations implicit in this scheme, the neglect ii of singlet operators, a generalized RU methodology was implemented, which includes a fully unrestricted treatment with both singlet and triplet operators. The small magnitude of the changes in hyperfine coupling constants confirms the validity of the original scheme. / QC 20100923