result[i, j] = arr[i:j+1].var() * (j-i)

import numpy

a = numpy.array([1, 2, 3, 4])