function P = makepos(P)
P = pnormdim(pl(P)+1e-10,1);
return