p <- makeFun(orings.model)(31)
1 - (1-p)^(1/6) -> q; q # P(damage to particular O-ring)
1 - dbinom(0, 6, q) # P(damage to >0 O-rings)
cbind(0:6, dbinom(0:6, 6, q)) # table of all probabilities
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.