Appendix MATLAB codes for the 3D three-point bending beam case