pt_a = pnorm(0.38) - pnorm(-2.18) pt_b = pnorm(2.72, mean=1.35, sd=2.2) - pnorm(1.16, mean=1.35, sd=2.2) print(c(pt_a, pt_b))