Вы можете передать аргумент ddof при выполнении np.std ()
a=[1,2,3,4,5]
vsum=0
for x in a:
vsum=vsum + (x - np.mean(a))**2
variance=vsum / (len(a) - 1)
print("Variance is {} and Actual STD Dev is {} ".format(variance,np.sqrt(variance)))
print("np.std is {} and np.std with parameter is {} ".format(np.std(a),np.std(a,ddof=1)))
То, что делает ddof = 1, заставляет numy использовать дельта степеней свободы до 1