#1. Distribusi Eksponensial 
#Untuk P(X=4)
eks<-function(x,lambda){lambda*exp((1/6)*lambda*x)}
eks(4,1/6)
## [1] 0.1862532
#Untuk P(X<4)
cdfexp <- function(x, lambda = 1/6) {
  pdf <- function(t) {
    lambda * exp(-lambda * t)
  }
  hasil<- integrate(pdf, lower = 0, upper = x)
  return(hasil$value)
}
cdfexp(4,1/6)
## [1] 0.4865829
#2. distribusi uniform
uniform<-function(x,a=0,b=1) {1/(b-a)}
uniform(6,a=2,b=30)
## [1] 0.03571429
dunif(6,min=2,max=30)
## [1] 0.03571429
#3. Distribusi Normal
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