Simple FEM analysis, defining boundary conditions

Hello,

I'd like to perform a simple FEM analysis and to estimate deformation of a squared beam. I'm following an example in the Matlab documentation, but I need to define boundary conditions a little bit differently.

Please have a look at the attached image. I want beam to rest on static points (mount points) on lower points V1, V2, V3 and V4. The forces are acting on two points in the middle of a beam, on its top.

I want to estimate a bent deformation.

Problem is how to use functions structuralBC() tp define the static points and structuralBoundaryLoad(). I have seen an example where whole face is defined as fixed, but this is not the case here, since lower face will also bent.

structuralBC(smodel,...
structuralBoundaryLoad(smodel,...