曲线绘制：
ker_density2=function(x,h){  
  x=sort(x)  
  n=length(x);t=0;y=0  
  for(i in 1:n){  
    t[i]=0;
    for(j in 1:n)  
      t[i]=t[i]+exp(-((x[i]-x[j])^2)/(2*h*h))  
  }  
  for(i in 1:n)  
    y[i]=t[i]/(n*h*sqrt(2*pi))  
  z=complex(re=x,im=y)  
  hist(x,ylim = c(0,0.006),breaks = 50,freq=FALSE)  
  #  hist(x,freq=FALSE)  
  lines(z)  
}  
ker_density2(s,11.21)

具体某个点的值：
density2=function(x,h,tn){  
  x=sort(x)  
  n=length(x)
  t=0
  for(j in 1:n)  
    t=t+exp(-((tn-x[j])^2)/(2*h*h))
  
  print(t/(n*h*sqrt(2*pi)) )
}  
density2(s,11.21,0)

请问您解决了吗