function pNx=parzen(N,h1,x)%h1對估計量有很大的影響
%Parzen
hN=h1/sqrt(N);
pNx=zeros(1,5000);
for u=1:5000
    for i=1:N
        pNx(u)=pNx(u)+exp(((x(u)-x(i))/hN).^2/-2)/sqrt(2*pi)/hN;
    end
    pNx(u)=pNx(u)/N;
end