def METHOD_NAME(zcashd, deltas, height, flag):
    theoreticalSupply = Network(MAINNET).SupplyAfterHeight(height)
    block = zcashd.getblock(str(height), flag)
    measuredSupply = block['chainSupply']['chainValueZat']
    empiricalMaximum = measuredSupply + deltas.DeviationUpToHeight(height)
    return (theoreticalSupply, empiricalMaximum, block)