Shear correction factor ( k = 5/6 ) is applied.

For uniform pressure p (N/m²):

% Find center deflection center_x = floor(nx/2)+1; center_y = floor(ny/2)+1; w_center_FEM = W(center_x, center_y);