P⁡(a,z)=12⁢erfc⁡(-η⁢a/2)-S⁡(a,η),