(z±2-z∓2)⁢RD⁡(x2,y2,z2)=2⁢(z±2-a2)⁢RD⁡(a2,z∓2,z±2)-3⁢RF⁡(x2,y2,z2)+(3/z),