Prak Komstat Pert 4

norm<-function(x,miu,sigma){1/(sigma*sqrt(2*pi))*exp(-1/2*((x-miu)/sigma)^2)}
norm (750,800,50)

## [1] 0.004839414

dnorm(750,mean=800, sd=50)

## [1] 0.004839414

#p(X=15)
eks<-function(x,lamda){lamda*exp((-1)*lamda*x)}
eks(15,0.05)

## [1] 0.02361833

#P(X<_15)
cdfexp <- function(x,lamda = 0.05){
  pdf<- function(t){
    lamda * exp (-lamda*t)
  }
  hasil<-integrate (pdf, lower= 0, upper=x)
  return(hasil$value)
}
cdfexp(15,0.05)

## [1] 0.5276334

uniform <- function(x, a=0, b=1) {
  1/(b-a)
}

# P(X = 6), durasi tepat 6 menit
uniform(6, a=2, b=10)

## [1] 0.125

dunif(6, min=2, max=10)

## [1] 0.125