Расхождения в графиках

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import scipy.stats as sts
%matplotlib inline
from scipy.stats import pareto
import math
#задали постоянные величины
xm=1
k=3
E=1.5
D=6.75
n=10
samp10=np.array([])#обозначим массив выборочных средних
i=1
while i in range (1,1001):#цикл генерации случайных выборок количества n из функции
    samp10i = pareto.rvs(n)#генерация выборки
    mean10i=sum(samp10i)/n#определение выборочного среднего (в.с.)
    samp10=np.append(samp10,mean10i)#добавляем в массив новое значение в.с.
    i+=1
x=np.linspace(1,3,10)
#определим параметры нормального распределения описывающие массив 1000 в.с. при n=5
D10=D/n#дисперсия для нормального распределения приближающего выборку средних
sigma10=math.sqrt(D10)#сигма для нормального распределения выборки средних
norm10 = sts.norm(E, sigma10)#определяем нормальное распределение с расчетными характеристиками Парето
pdf10 = norm10.pdf(x)#получаем плотность распределения
#строим гистограмму массива средних при n=5 и плотность описывающую этот массив нормального распределения
plt.plot(x, pdf10, label='PDF(n=10)')
plt.hist(samp10, bins =100, range=(1,3), normed=True)
plt.ylabel('number of samples')
plt.xlabel('$x$')
plt.legend(loc='best')