data4.txt data5.txt data9.txt data22.txt data30.txt data118.txt data1047.txt iter.m main.m make_J.m make_Y.m read_data.m var_pq.m