ERF/ERFC误差函数
def erfc(e_val)local tal = 1. / (1. + .3275911 * abs(e_val))
local val = ((1.061405429*tal-1.453152027)*tal+1.421413741)*tal
val = 1. - ((val-.284496736)*tal+.254829592)*tal*exp(-e_val*e_val)
if e_val < 0.0 then
val = - val
end_if
erfc = 1. - val
end
def erf(e_val)
local tal = 1. / (1. + .3275911 * abs(e_val))
local val = ((1.061405429*tal-1.453152027)*tal+1.421413741)*tal
val = 1. - ((val-.284496736)*tal+.254829592)*tal*exp(-e_val*e_val)
if e_val < 0.0 then
val = - val
end_if
erf = val
end
页:
[1]