We consider the operator L=−(d/dx)2+x2y+w(x)y,y∈L2(R), where w(x)=s[δ(x−b)−δ(x+b)],b=0, real, s∈C. This operator has a discrete spectrum: eventually the eigenvalues are simple and λn=(2n+1)+s2(κ(n)/n)+ρ(n), where κ(n)=2π1[(−1)n+1sin(2b2n)−21sin(4b2n)] and ∣ρ(n)∣≤C(logn)/(n3/2) If s=iγ, γ real, the number T(γ) of non-real eigenvalues is finite, and T(γ)≤[C(1+∣γ∣)log(e+∣γ∣)]2. The analogue of the above equations is given in the case of any two-point interaction perturbation w(x)=c+δ(x−b)+c−δ(x+b),c+,c−∈C.