add w_k output in mbarpmf.m 2