function yhat=expect(prob) 
% probability of y, yhat is
% expected value of the output
%
  my=max(size(prob));
  yhat=0;
  for y=1:my
    yhat=yhat+y*prob(y);
  end