plot tr(P)