Need help for Matlab code: Method of Moments: Conducting Plate

Hi, I am trying to write a Matlab code that works for charged conducting Plate. Spesifically, the problem is from R.F. Harringston's "Field Computation by Moment Method" Ch2.2

I have following codes, but that does not work for the cases if I choose number of subsections more than one.

Here are codes:
I think the problem is in the function MatrixL2x2. Because I calculate the value for the case #subsection = 1

Here are the actual results:
N = 1 => 31.5 picofarads/meter
N = 3 => 37.3 picofarads/meter
N = 4 => 38.2 picofarads/meter
N = 6 => 39.2 picofarads/meter

As I said, the code can find N=1 case result.
Please ignore comments.
Thanks.

I solved the problem, the topic can be closed.