Post date: May 09, 2018 7:6:8 PM
# pressurized tank fault tree using pbox.r, a probability bounds analysis library for R
source("C:\\Users\\ferson\\Google Drive\\workshop sra\\Arlington 2017\\pba.r")
cb = function(k,n) return(beta(interval(k, k+1), interval(n-k, n-k+1)))
s = cb(1,150)
ss = pmin(s,.1)
t = cb(0,2000)
tt = pmin(t,.005)
k2 = cb(3,500)
s1 = cb(0, 460)
k1 = cb(7,10000)
r = cb(0, 380)
e = t %|||% (k2 %|||% (s %|&|% (s1 %|||% (k1 %|||% r))))
ef = t %|% (k2 %|% (s %&% (s1 %|% (k1 %|% r))))
pmin(e, 0.03)
pmin(ef, 0.03)