Intermixing in heteroepitaxial islands: fast, self-consistent calculation of the concentration profile minimizing the elastic energy