Ra(𝐛;𝐳)=2B(b1,b2)0π/2(z1cos2θ+z2sin2θ)a(cosθ)2b11(sinθ)2b21dθ,