We present the program SuSpect3 that calculates the masses and couplings of the Higgs and supersymmetric particles predicted by the Minimal Supersymmetric Standard Model (MSSM). The model is implemented in both its non-constrained version, the MSSM, and its constrained versions, such as the minimal supergravity and the gauge or anomaly mediated supersymmetry breaking models, in which the soft supersymmetry–breaking parameters obey certain universal boundary conditions at the high energy scale. The low energy parameters are then obtained using renormalization group equations and electroweak symmetry breaking, and all the dominant radiative corrections have been consistently implemented. SuSpect3 is a major rewrite, in C++ object oriented programming, of the FORTRAN code SuSpect. It includes all the features of the earlier code in an improved and updated manner, and involves new options such as compressed SUSY scenarios, an MSSM-inflation model and the possibility of using the observed Higgs mass as an input. The main features and the use of the program are explained. Program summaryProgram Title:SuSpect3CPC Library link to program files:https://doi.org/10.17632/55jv4kdrm8.1Developer's repository link:http://suspect.in2p3.frLicensing provisions: GPLv3Programming language:C++, compatible C++98, C++11, C++14, C++17Nature of problem: Supersymmetric models such as the MSSM, mSUGRA, GMSB, AMSB and others have specific parameter sets and boundary conditions. SuSpect3 translates the parameter sets of the models into predictions of the Higgs and supersymmetric particles masses. The mixing matrices of the physical states as well as the mixing angles are calculated in addition to the scale dependent parameters.Solution method: The spectrum of the Higgs and supersymmetric particles depends on the model, its supersymmetric parameter set and the Standard Model parameters. The evolution of the parameters as function of the energy scale is calculated by solving numerically the Renormalization Group Equations. Model dependent boundary conditions are applied at the appropriate scale. Electroweak symmetry breaking is calculated iteratively at the electroweak scale. The application of radiative corrections translate the scale dependent particle masses into the physical pole masses.Additional comments including restrictions and unusual features: The parameters and physical masses are defined to be real. Warnings are issued if intermediate results are unphysical.Compiles with gcc-4.8.5 and later (checked with gcc-8.1.0).