Create calc_prob.m