function myerf, x, p

; P should be [mean, variance]

; p contains the mean and sigma values.

val = 0.5 * (1.d + erf( ( x - p[0]) / ( p[1] * sqrt(2) )))

return, val

end