#Distribus normal (PDF) ada nilai rata",standar deviasi
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) #fungsi bawaan dari r untuk peluang(altarnatifcara cepat), d untuk nyari peluang tpeat,p kumulatif
## [1] 0.004839414
#distribusi eksponial#waktu tunggu
#Untuk P(X=15)
eks<-function(x,lamda){lamda*exp((-1)*lamda*x)}
eks(15,0.05)
## [1] 0.02361833
#untuk 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
#Distribus 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