def f_ref(x):
    if (x < 2):
        return (np.array([1.0, 2.0]) * x)
    else:
        return (np.array([3.0, 4.0]) * np.sin(x))
def f_ref(x):
    if (x < 2):
        return (np.array([1.0, 2.0]) * x)
    else:
        return (np.array([3.0, 4.0]) * np.sin(x))