n
rel
ω
This GPU based script draws 10 cross-sections X after every
_{i,n/2}2rel weighted Jacobi iterations.
Fx.
_{i,j} = Fy_{i,j} = sin(πi/n) sin(πj/n) |

where

See also 2D weighted Jacobi relaxation on CPU.

[1] Linear elasticity (Wiki)

Simulations on GPU