def METHOD_NAME(x):
    f_order = 3
    Pth_order = 1
    constants = np.zeros((6, 6, f_order + 1, Pth_order + 1))
    S11_0 = x[0]
    dS11df = x[1]
    d2S11df2 = x[2]
    dS11dPth = x[3] * 1.0e-11
    d2S11dfdPth = x[4] * 1.0e-11
    S44_0 = x[5]
    dS44df = x[6]
    d2S44df2 = x[7]
    dS44dPth = x[8] * 1.0e-11
    d2S44dfdPth = x[9] * 1.0e-11
    S12_0 = (1.0 - 3.0 * S11_0) / 6.0
    dS12df = -dS11df / 2.0
    d2S12df2 = -d2S11df2 / 2.0
    dS12dPth = -dS11dPth / 2.0
    d2S12dfdPth = -d2S11dfdPth / 2.0
    constants[:3, :3, 1, 0] = S12_0
    constants[:3, :3, 2, 0] = dS12df
    constants[:3, :3, 3, 0] = d2S12df2
    constants[:3, :3, 0, 1] = dS12dPth
    constants[:3, :3, 1, 1] = d2S12dfdPth
    for i in range(3):
        constants[i, i, 1, 0] = S11_0
        constants[i, i, 2, 0] = dS11df
        constants[i, i, 3, 0] = d2S11df2
        constants[i, i, 0, 1] = dS11dPth
        constants[i, i, 1, 1] = d2S11dfdPth
    for i in range(3, 6):
        constants[i, i, 1, 0] = S44_0
        constants[i, i, 2, 0] = dS44df
        constants[i, i, 3, 0] = d2S44df2
        constants[i, i, 0, 1] = dS44dPth
        constants[i, i, 1, 1] = d2S44dfdPth
    return AnisotropicMineral(per, cell_parameters, constants)