PrakKomstat_Pert4

#ada kata tepat
#normal ada rata-rata dan standar deviasi
norm<-function(x,miu,sigma){
  1/(sigma*sqrt(2*pi))*exp(-1/2*((x-miu)/sigma)^2)
}
norm(750,800,50) #800 rata rata #750 yang ditanya # sd=50

## [1] 0.004839414

dnorm(750,mean= 800, sd=50)

## [1] 0.004839414

#eksponensial
#rata-rata waktu tunggu #rata-rata:20-> 1/20=0.05, #tepat 15 menit dan tepat setelah 15 menit
#untuk P(X=15)
eks<-function(x, lamda)(lamda*exp((-1)*lamda*x))
eks(15,0.05)

## [1] 0.02361833

#untuk P(X<=15)
cdfexp<-function(x,lamda=0.05){
  pdf<-function(t) {
    lamda*exp(-lamda*t)
  }
  hasil <-integrate(pdf,lower=0,upper=x)
  return(hasil$value)
}
cdfexp(15,0.05)

## [1] 0.5276334

#uniform dan ada interval
uniform<-function(x, a=0,b=1)(1/(b-a)) #tepat 20 #interval a=10, b=30
uniform(20,a=10,b=30) #harus pake yang ini

## [1] 0.05

dunif(20,min=10,max=30)

## [1] 0.05