Based on the Volterra integral technique, a universal method to calculate the sorption uptake curves of multicomponent gas mixtures in porous solids under variable high pressure conditions with non-ideal behavior of the gas phase has been proposed. It takes into account a non-linearly behaving concentration dependence at the surface of the porous medium. The applicability of the introduced method presupposes an isothermic gas phase and a constant diffusivity matrix. For the description of the non-ideal gas phase behavior, equations of state were utilized.where Do stands for the self-diffusivity, E, T, Ii and a"" denote the activation energy of intraparticle diffusion, the temperature, the average sorbate concentration and saturation value of the sorbate concentration, respectively, and D