I have implemented a constraint dynamics physics simulation as proposed by Andrew Witkin et al 1990, but I cannot get the initial constraint “snapping” correctly.

I am creating a game where you manage a spaceship, I want a main feature of it to be at least partially accurate gas diffusion that looks something like the diffusion in the game “oxygen not included” except top down.

Good evening to everyone please i need help. I have these functions bellow:

The allocation of 5 plots to 8 individuals is as follows, with the priority and superiority of the plots arranged in this order:
Second plot > Third plot > Fifth plot > Fourth plot > First plot
Also, the first and second priority sections of the second, third, and fourth plots are located in the same position.
The goal is to assign plots to individuals so that each one can obtain their best possible option, with the priority of the plots arranged from top to bottom, and the higher plots having a greater value and price. However, to ensure fairness, the sizes of the plots are varied so that the different values of the plots are equalized.
Each land has 8 plots:
The first and fifth lands have one plot per row,
The other lands have two plots in the first row and one plot in the remaining rows. python code