I am trying to write a finite difference code in Matlab to solve the Kirchhoff's plate equation. So far, I have been able to generate the biharmonic/bilaplacian operator matrix. However, I am struggling with how to impose the boundary condition. I would be glad is someone can help explain how to achieve this.
I have taken the liberty of attaching here an m.file I found on the net prepared by
I. Danaila, P. Joly, S. M. Kaber & M. Postel. The task is to use this to generate the biharmonic operator matrix of the kirchhoff's plate equation (which I also have taken the liberty to attach here as well). Thanks