h_retail: h_retail
In CRC1266-A2/moin: Modelling Interactions

Retail Gravity model

1	h_retail(E, Si = 1, Sj = 1, beta, C, c, gammas, g, delta, x, s)

`E`	a matrix of random values (to be "optimized" by the Hamiltonian)
`Si`	a vectory of node attributes, e.g. carrying capacity or population
`Sj`	a vectory of node attributes, e.g. carrying capacity or population
`beta`	lagrange multiplier
`C`	distance matrix
`c`	distance constraining cost of the system
`gammas`	lagrange multiplier vector
`g`	constraining vector
`delta`	lagrange multiplier
`x`	general node constraining value
`s`	a vector of node characteristics