Развит полуэмпирический метод расчета компонентов тензора поляризуемости двухатомных гомоядерных молекул как функции межъядерного расстояния R. Метод базируется на установленной нами ранее функции поляризуемости молекулы при малых R, значениях поляризуемости и их производных в равновесном положении ядер молекулы и уточненной функции поляризуемости взаимодействующих ориентированных атомов при больших R, которая включает в себя учет изменения поляризуемости атомов при их сближении, а также учитывает мультипольные взаимодействия между ними. Учет мультипольного взаимодействия атомов проведен в рамках модели поляризуемости пары взаимодействующих диэлектрических шаров. Предложенный метод проверен на молекуле H2 и затем применен для расчета функций поляризуемости молекул N2 и O2.