% calcolo della matrice delle prob a posteriori Pp
% px e' un vettore colonna 
function Pp=apost(px,Pc)
py=Pc'*px;
Pp=diag(px)*Pc*inv(diag(py));