def METHOD_NAME(mesh_interval, family, degree):
    V = VectorFunctionSpace(mesh_interval, family, degree, dim=2)
    x = SpatialCoordinate(mesh_interval)
    f = Function(V).interpolate(as_vector(((x[0] - 0.5)*(x[0] - 0.5), x[0]*x[0])))
    assert np.allclose([0.01, 0.36], f([0.6]))
    assert np.allclose([0.25, 1.00], f([1.0]))