source("91.r")
# zbog funkcije histogram

oldFaithful = read.table("baze//oldFaithful.txt")
podaci = oldFaithful$V1
n = length(podaci)

histogram(podaci, "Sturdzis")

# naivna ocena gustine
h = 0.25

Ind = function(x) {
  if(abs(x) < 1)  {
    return(1)
  }
  else {
    return(0)
  }
}

range(podaci)
# vrednosti u kojima racunamo ocenu gustine
# malo siri interval vrednosti od opsega podataka
y = seq(1, 6, 0.01)

f = c()
for (i in 1:length(y)) {
  suma = 0
  for (j in 1:n) {
    suma = suma + Ind((podaci[j] - y[i])/h)
  }
  f[i] = suma/(2*n*h)
}

plot(y, f, type = 'l', xlab = "podaci", ylab = "", main = "Naivna ocena gustine")