The scaled error function in quartic pricing model that encaptulates
both scaled erfi and erfc functions
into a single representation. This is used to provide an elegant expression for
the MGF and local option prices, \(L_{c,p}\).
When sgn=-1, it is \(\sqrt{\pi}e^{-x^2} erfi(x)\), which twice of Dawson function.
When sgn=1, it is \(\sqrt{\pi}e^{x^2} erfc(x)\).
ecd.erfq_sum is the summation implementation with truncation rule set forth
in the quartic pricing model. It achieves high precision when x > 4.5.