A numerical simulation for body surface potential mapping was investigated using Boundary Element Method (B.E.M.). After forward simulation, we investigated the accuracy and stability of ECG inverse solution from the viewpoint of Minimum Description Length Principle (MDL). It became obvious that the solution containing up to 16th expansion components is optimal for ECG inverse problem.