Elliot Tanis
April 20, 2006
Exercise 5.3-1  (Exercise_5_3-1.mws)

>	X := 'X': XX := 'XX':

>	f := 4*x^3;

>	mu := int(x*4*x^3, x = 0 .. 1); evalf(%);

>	var := int((x - mu)^2*4*x^3, x = 0 .. 1); evalf(%);

Here are the random numbers from Table IX.

>	Y := [.7510,.0543, .4089, .8130, .0345, .3640, .5469, .1788, .3131, .3937, .8718, .8826, .3322, .2676, .4747, .4364, .9846, .9317, .5799, .8396];

>	for k from 1 to 20 do   XX[k] := Y[k]^(1/4): od: X := [seq(XX[k], k = 1 .. 20)];

>	xbar := Mean(X); mu := evalf(mu);

>	svar := Variance(X); var := evalf(var);

Here are the observations of X  rounded off to 4 places.

>	X := [.9309, .4827, .7997, .9496, .4310, .7767, .8600, .6503, .7480, .7921, .9663, .9693, .7592, .7192, .8301, .8128, .9961, .9825, .8726, .9572];

>	xbar := Mean(X);

>	svarx := Variance(X);

>