yfunction fn=f1(y0)
y=[n个数]
f=[];n=509;hn=707;
for i=1:1:509;
    f(i)=(1/(2*pi)^(1/2))*exp(-((y0-y(i))./hn).^2/2);
end
fn=1/(n*hn)*sum(f(i));
要求：这个计算核密度估计的函数对不对，则呢么对fn求积分？？