pvec3 = np.linspace(b, pmax, nd4 + 1).ravel()