In inverse electrocardiography, the goal is to estimate cardiac electrical sources from potential measurements on the body surface. It is by nature an ill-posed problem, and regularization must be employed to obtain reliable solutions. This paper employs the multiple constraint solution approach proposed in Brooks et al. (IEEE Trans Biomed Eng 46(1):3-18, 1999) and extends its practical applicability to include more than two constraints by finding appropriate values for the multiple regularization parameters. Here, we propose the use of real-valued genetic algorithms for the estimation of multiple regularization parameters. Theoretically, it is possible to include as many constraints as necessary and find the corresponding regularization parameters using this approach. We have shown the feasibility of our method using two and three constraints. The results indicate that GA could be a good approach for the estimation of multiple regularization parameters.