Hello, I am trying to obtain deflection results on a beam made with shell elements

I wrote a MATLAB code. The problem is that kk is a matrix with 3000*3000 full of zeros (it´s a singular matrix), and it is not supposed to beThe code runs well until I call the function feasmbl1.m, which is shown next