def METHOD_NAME():
    x = SpatialCoordinate(triangle)
    xx = outer(x, x)
    m = as_matrix(xx[i, j], (i, j))
    s = m[0, 0] * m[1, 0] * m[0, 1] * m[1, 1]
    e = s((5, 7))
    v = 5 * 5 * 5 * 7 * 5 * 7 * 7 * 7
    assert e == v