Application of the cluster variation method to interstitial solid solutions