def METHOD_NAME():
    plane = NumberPlane(x_range=[2, 4])
    origin = plane.x_axis.number_to_point(
        plane._origin_shift([plane.x_axis.x_min, plane.x_axis.x_max]),
    )
    def ref_func(*coords):
        result = np.array(origin)
        for axis, number in zip(plane.get_axes(), coords):
            result += axis.number_to_point(number) - origin
        return result
    coords = [[1], [1, 2], [2, 2], [3, 4]]
    print(f"\n\nTesting coords_to_point {coords}")
    expected = np.round([ref_func(*coord) for coord in coords], 4)
    actual1 = np.round([plane.coords_to_point(*coord) for coord in coords], 4)
    coords[0] = [
        1,
        0,
    ]  # Extend the first coord because you can't vectorize items with different dimensions
    actual2 = np.round(
        plane.coords_to_point(coords), 4
    )  # Test [x_0,y_0,z_0], [x_1,y_1,z_1], ...
    actual3 = np.round(
        plane.coords_to_point(*np.array(coords).T), 4
    )  # Test [x_0,x_1,...], [y_0,y_1,...], ...
    print(actual3)
    np.testing.assert_array_equal(expected, actual1)
    np.testing.assert_array_equal(expected, actual2)
    np.testing.assert_array_equal(expected, actual3.T)