Numerical solution of urban residential land-use models