add w_k output in mbarpmf.m 4