We review the long term project of the ALPHA collaboration to compute in QCD the running coupling constant and quark masses at high energy scales in terms of low energy hadronic quantities. The adapted techniques required to numerically carry out the required multiscale non-perturbative calculation with our special emphasis on the control of systematic errors are summarized. The complete results in the two dynamical flavor approximation are reviewed and an outlook is given on the ongoing three flavor extension of the programme with improved target precision.