#Distribusi Normal
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
#Distribusi Eksponensial

#Untuk P(X=15)
eks<-function(x,lambda){lambda*exp((-1)*lambda*x)}
eks(4,0.16)
## [1] 0.08436679
#Untuk P(X<=15)
cdfexp<-function(x,lambda=0.5)  {
  pdf<-function(t) {
    lambda*exp(-lambda*t)
  }
  hasil<-integrate(pdf, lower=0, upper=x)
  return(hasil$value)
}
cdfexp(15,0.05)
## [1] 0.5276334
#Distribusi Uniform
uniform<-function(x,a=10,b=30){1/(b-a)}
uniform(20,a=10,b=30)
## [1] 0.05
dunif(20, min=10, max=30)
## [1] 0.05

```