ABD matrix calculation in Matlab

Hi All

if someone has worked on composites , could you tell me why in this code , in the line 104 , the q33 is doubles and in the line 108 it is devided by 2 again , all the third collumn elements ?

and I think at the end ABD matrixes are not calculated correctly , cause the results overlap the previous line, in the for loop after line 122

[farzad]

sorry here is the link :

http://www.mathworks.it/matlabcentral/fileexchange/25567-abd-matrix/content/ABD_mat_revised