`def normpdf(x, mean, sd):
    var = float(sd)**2
    pi = 3.1415926
    denom = (2*pi*var)**.5
    num = math.exp(-(float(x)-float(mean))**2/(2*var))
    print("num: ",num)
    print(num/denom)
    return num/denom


  x_v1 = np.array([0.96450007,0.96450007,0.96450007,0.96450007,0.96450013,0.96450013,0.96450013,0.96450013,0.96450013,0.96450013,0.96450013,0.96450013])

x_v=np.sort(x_v1)

mean=np.mean(x_v)
sd=np.std(x_v)
print("mean: ",mean)
print("SD: ",sd)
prob=[]

for i in range(len(x_v)):
    x=x_v[i]
    vv=normpdf(x,mean,sd)
    prob.append(vv)
prob=np.array(prob)
print("prob_value: ",prob)

plt.plot(x_v,prob,marker='o')
plt.savefig("./fig")`