Numerical integration of the Gibbs–Thomson equation for multicomponent systems