% file ygivenx.m
% Version of 12/11/96
% X, PX, and PY given X
% Determine joint dbn
X  = input('Enter VALUES for X  ');
PX = input('Enter PROBABILITIES for X  ');
Y  = input('Enter VALUES for Y  ');
py = input('Enter string variable for PY|X  ');
m  = length(X);
n  = length(Y);
P0 = zeros(m,n);
for i = 1:m
  pY = eval(py);
  P0(i,1:length(pY)) = PX(i)*pY;
end
P = rot90(P0);
disp('Results are in X, Y, P')