Looking back at your original xmcd file, the problem with P0(f) is that you are taking the determinate rather than absolute value. Unfortunately, Mathcad uses vertical bars for both so you can't tell which it's doing by looking at a printout. Point the mouse at the S21_0 term in P0(f) and right click, if the box for determinate is checked, change it to absolute value. Note that both the matrix toolbar and the calculator toolbar have a button for vertical bars, matrix toolbar is determinate, calculator toolbar is absolute value.
For your question about Q0, it will be easiest to work with the xmcd file. Please attach that to your next post.
Peter