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