function sigmoid, x0, numPoints=numPoints, width=width, beta=beta ;returns a sigmoid according to the imput parameters if not keyword_set(width) then width=1 if not keyword_set(beta) then beta = 0.1 if not keyword_set(numPoints) then numPoints = 100 x = findgen(numPoints) result = 1.0/( 1 + exp((2/width)*alog10(beta/(1 - beta))*(x - x0))) return, result end