fixed loading matfile with .bak_* extension; in fig1 computation, N_vec now includes endpoints