写了一个python的函数，不知道如何弄了

from scipy.stats import norm
u=float(0.68155),s=float(0.180658),TEa=float(0.02)
#u表示平均值，s表示标准偏差，TEa规定的允许偏差
def Pe (SE):
    xs=round([u+uniform(-r/2 ,  r/2 )+SE]，r)
#四舍五入保留r位小数
    b=1-norm.cdf(((u+TEa)-xs)/s)+norm.cdf(((u-TEa)-xs)/s)
    print(b)
#得出一个关于SE的Pe函数