M:matrix([c[p],c[p],c[x]+c[u]],[c[u],c[x]+c[p],c[u]],[c[u]+c[p],c[x],c[x]]);

determinant(M);

expand(%);