This form uses the classical theory of thin rectangular plates in bending as represented in the formulae of the form.
The solution is obtained with a finite difference scheme using a grid 64x64, with subsequent interpolation over the ratio a/b.
For the concentrated load the stresses at center correspond to a load concentrated over an area roughly equal to ab/64²