eks<-function(x, lamda)(lamda*exp((-1)*lamda*x))
eks(4,0.167)
## [1] 0.08562641
#untuk P(X<=4)
cdfexp<-function(x,lamda=0.167){
  pdf<-function(t) {
    lamda*exp(-lamda*t)
  }
  hasil <-integrate(pdf,lower=0,upper=x)
  return(hasil$value)
}
cdfexp(4,0.167)
## [1] 0.487267
uniform <- function(x,a=0,b=1) {1/(b-a)}
uniform(6, a = 2 , b =10)
## [1] 0.125
dunif(6,min=2,max = 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