lambda <- 1/6
eks <- function(x, lambda) (lambda * exp((-1) * lambda * x))
cdfexp <- function(x, lambda) {
pdf <- function(t) {
lambda * exp(-lambda * t)
}
hasil <- integrate(pdf, lower = 0, upper = x)
return(hasil$value)
}
eks(4, lambda) # 0.08556
## [1] 0.08556952
cdfexp(4, lambda) # 0.4865829
## [1] 0.4865829
uniform<-function(x,a=0,b=1){1/(b-a)}
uniform(6,a=2,b=10)
## [1] 0.125
norm = function(x,miu,sigma){1/(sigma*sqrt(2*pi))*exp(-1/2*((x-miu)/sigma)^2)}
norm(4,5,1)
## [1] 0.2419707
dnorm(4,mean=5,sd=1)
## [1] 0.2419707