Distribusi peluang kontinue

 #Uji Normal pdf
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

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

## [1] 0.02361833

 #untuk P
 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 
uniform<-function(x,a=0,b=1){1/(b-a)}
uniform(20,a=10,b=30)

## [1] 0.05

dunif(20,min=10,max=30)

## [1] 0.05