POSTTEST 1 KOMSTAT

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

#untuk P(X = 15) pdf
eks <- function(x, lambda) {lambda*exp ((-1) * lambda*x)}
eks (4, 1/6)

## [1] 0.08556952

#Untuk P(X<= 15) cdf 
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

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