Esum := proc() ndraws*AVEbox; end; SEsum := proc() sqrt(ndraws)*SDbox; end; SEave := proc() SDbox/sqrt(ndraws); end; SEpercent := proc() SDbox*100/sqrt(ndraws); end; su42 := (42-40)/SEpercent(); ndraws := 1000: SDbox := sqrt(.4*.6): 'su42' = evalf(su42,3); chance := (100 - 80.64)/2.; ndraws := 5000: 'SEpercent'= evalf(SEpercent(),3); su42 := (42-40)/0.691; chance := (100 - 99.63)/2.; CF := proc() sqrt( (N - ndraws)/(N-1) ); end; ndraws := 1000: N := 6000: correction_factor := evalf(CF(),4); NewSE := evalf( SEpercent()*correction_factor, 4); NewSU42 := (42-40)/1.414; NewChance := (100 - 84)/2.; ndraws := 5000: correction_factor := evalf(CF(),4); NewSE := evalf( SEpercent()*correction_factor, 4); NewSU42 := (42-40)/0.2828;