w := p -> 1-sum(binomial(20,k)*p^k*(1-p)^(20-k),k=0..6)+ sum(binomial(20,k)*p^k*(1-p)^(20-k),k=0..1); seq(evalf(w(j/10),3),j=1..9); evalf(w(0.01),3); evalf(w(0.001),3); plot(w(p),p=.001..0.999);