#ada rata rata dan sd
norm <- function (x,miu,sigma) {1/(sigma*sqrt(2*pi))*exp(-1/2*((x-miu)/sigma)^2)}
norm(750,800,50)
## [1] 0.004839414
dnorm(750,mean = 800,sd = 50)
## [1] 0.004839414
#punya ratarata waktu tunggu
#exponensial
eks<-function(x,lambda){lambda*exp((-1)*lambda*x)}
eks(15,0.05)
## [1] 0.02361833
#untuk P
cdfexp <- function(x, lambda = 0.05) {
pdf <-function(t) {
lambda * exp(-lambda * t)
}
hasil <- integrate(pdf, lower = 0, upper = x)
return(hasil$value)
}
cdfexp(15,0.05)
## [1] 0.5276334
#punya interval uniform
uniform<- function(x,a=0,b=1){1/(b-a)}
uniform(20,a=10,b=30)
## [1] 0.05
dunif(20, min=10,max=30)
## [1] 0.05
1/6