We present an ab initio approach for the computation of the magnetic s
usceptibility chi of insulators. The approach is applied to compute ch
i in diamond and in solid neon using density functional theory in the
local density approximation, obtaining good agreement with experimenta
l data. In solid neon, we predict an observable dependence of chi upon
pressure.