def METHOD_NAME(language):
    def f(a : 'float[:,:,:,:]'):
        return sum(sum(a[i,k,4,2] for i in range(5)) for k in range(5))**2
    x = randint(0, 10, size=(5,5,5,5)).astype(float)
    f_epyc = epyccel(f, language = language)
    assert f(x) == f_epyc(x)