First, here is my code (the last 3 lines are actually a single expression, despite being split over multiple lines below):

with(Statistics):
with(Student):
with(Student[NumericalAnalysis]):
Student[Calculus1][Roots]:
k := 8;
size:=5000;
h:=3.1;
Llimit := evalf((size-1)-10*sqrt(2*(size-1)));
Uplimit := evalf((size-1)+10*sqrt(2*(size-1)));
evalf(Int((Int(CDF(Normal(0, 1), h/sqrt((size-1)*(1/x+1/y)))*PDF(RandomVariable
(ChiSquare(size-1...