I tested the new release and it is agreeing with other codes now. Michael, you may also wish to verify that the new version is working correctly.

I also think adding the cut off points as an additional x or y axis would be very useful to understand the trade off between sensitivity and sensibility.

Nice function, but I think it may have a bug, I am getting different results for significantly overlaping distributions when I compare to medcalc, and online calculators.

I made an example to use my code. Assume you have a distribution from -20 to 20 with mean at 0, and you want to calculate the area from -inf to 0
Example:
x=[-20:20] % your data points
sigma=length(x)/2/3.5 % standard div. normal dist. with is 3.5
mean=0 % mean of distribution
normaldistribution(0, mean, sigma,1)

Comment only