We present a Bayesian method that simultaneously takes all available observables for a given star, their uncertainties and prior knowledge like the initial mass function into account to determine probability distributions of the stellar parameters based on grids of stellar evolutionary models. This allows to homogeneously analyse stars, to determine stellar parameters including robust uncertainties and to identify stars that cannot be reproduced by current stellar models. Our code is available through an easy-to-use web-interface.