A numerical stable method for integration of the multicomponent species diffusion equations