pvec1 = np.linspace(pmin, a, nd4 + 1).ravel()