其中x = np.array([-1.5, 0, 1, 2]) y = np.array([0.125, -1, 1, 9])，满足边界条件S‘(-1.5)=0.75，S’(2)=14，完整代码展示，打印出满足要求的函数

根据提供的边界条件和点集，我们可以使用numpy的linalg模块来解决三阶线性系统找到样条插值函数的系数。下面是完整的Python代码：

import numpy as np

# 给定的数据点
x = np.array([-1.5, 0, 1, 2])
y = np.array([0.125, -1, 1, 9])

# 边界条件的导数值
dydx_left = 0.75
dydx_right = 14

# 我们可以将插值问题表示为线性方程组的形式
# 对于样条插值，一般形式是 A*c = b，其中 A 是系数矩阵，b 是常数项向量
# 这里 c 是我们寻找的插值系数数组，A[i,j] = (x[j] - x[i])^3，b[i] = y[i] + 3*(x[i+1] - x[i])*(y[i+1] - y[i]) / (x[i+1] - x[i])

# 系数矩阵A
A = np.zeros((len(x) - 1, len(x) - 1))
for i in range(1, len(x)):
    for j in range(i):
        A[i-1, j] = (x[j] - x[i]) ** 3

# 常数项向量b
b = y[:-1] + 3 * (np.diff(x) * (np.roll(y, 1) - y[:-1])) / np.diff(x)

# 添加边界条件的列向量
dydx_vector = np.zeros(len(x) - 1)
dydx_vector[0] = dydx_left
dydx_vector[-1] = dydx_right

# 扩展A和b来包含边界条件
A_ext = np.column_stack([A, dydx_vector])
b_ext = np.append(b, [0])

# 解决线性方程组
c = np.linalg.solve(A_ext, b_ext)

# 插值函数S(X)
def S(X):
    result = np.zeros_like(X)
    for i in range(len(c)):
        result += c[i] * (X - x[i]) ** (i+1)
    return result

# 打印满足要求的函数
print("满足要求的插值函数S(X):")
print(S(np.linspace(-1.5, 2, 100)))  # 在整个区间绘制插值函数

注意，由于实际问题中可能存在非唯一解（即存在多个三次多项式可以满足边界条件），这段代码返回的是一个满足条件的具体解。如果你的问题涉及到其他情况，例如要求最小化某种误差或有其他约束，可能需要采用更复杂的方法如勒贝格积分法等。