```
# expansion of S1(x) = S(x*sqrt(pi/2)), see reference[5] page 1-8
```
better?

Same for fresnelc