function p=cholest(x,y)

% x= cholesterol level (0= low, ..., 4= highly elevated)
% y= weight in pounds

lambda= 2 + (y-150)/75;
norm= 1 + lambda + lambda^2/2 + lambda^3/6 + lambda^4/24;
p= lambda^x/(norm * prod(1:x));
