xa = np.linspace(ax, bx, inc)