posttestkomstat

# expo
eks <- function(x, lamda) {lamda * exp((-1) * lamda * x)}

# a. 
eks(4, 1/6)

## [1] 0.08556952

cdfexp <- function(x, lamda = 1/6) {
  pdf <- function(t) {
    lamda * exp(-lamda * t)
  }
  hasil <- integrate(pdf, lower = 0, upper = x)
  return(hasil$value)
}
# b.
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

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