2546961Plate mechanics lab: Kirchhoff-Love plates solved by numbl in the browserJeremy Magland 1function L = plate_laplacian(wp, dx, dy)
2% PLATE_LAPLACIAN Five-point Laplacian of a padded field.
3%
4% L = plate_laplacian(wp, dx, dy) takes an (ny+2) x (nx+2) array with one
5% ghost layer (see plate_pad) and returns the ny x nx second-order
6% approximation of w_xx + w_yy at every node of the original grid.
7% Rows index y and columns index x.
9L = (wp(3:end, 2:end - 1) - 2 * wp(2:end - 1, 2:end - 1) + wp(1:end - 2, 2:end - 1)) / dy^2 ...
10 + (wp(2:end - 1, 3:end) - 2 * wp(2:end - 1, 2:end - 1) + wp(2:end - 1, 1:end - 2)) / dx^2;
11end