Integral form of cdf in Python

I want to use a mathematical model of this kind. How do I find the exact integral form of cdf(x) which Python uses in this context?

from scipy import statsdef model(x, x0, s, A, c):    law = stats.norm(loc=x0, scale=s)    return A * law.cdf(x) + c