% LMex100202_4th.m
% Example 10.2.2 
% Larsen & Marx (2006) Introduction to Mathematical Statistics, 4th edition
% Written by Eugene.Gallagher@umb.edu 11/21/10
syms y;
int(sym('6*y*(1-y)'),0,0.25)
p1=eval(int(sym('6*y*(1-y)'),0,0.25))
int(sym('6*y*(1-y)'),0.5,0.75)
p3=eval(int(sym('6*y*(1-y)'),0.5,0.75))
int(sym('6*y*(1-y)'),0.75,1)
p4=eval(int(sym('6*y*(1-y)'),0.75,1))
P=exp(gammaln(6)-gammaln(4))*p1*p3^3*p4