minpiece <- function(a,n,l) 1-((a-2*l)/a)^n m <- 10000 a <- 220*10^6 l <- (1:m)*(a/2)/m p <- mapply(minpiece, a, 24, l) min(l[p>0.025])