function [f] = part_rhs(w, p)

[f, pf_pw] = part_pf(w, p);