I=b*h^3/12;
i=1:12;
j=1:12;
Phi_i=2*pi.*(i-1)./n;
Phi_j=2*pi.*(j-1)./n;
Phi=Phi_j-Phi_i;
QCij=(R^3/(4*pi*E*I)).*(Phi.*sin(Phi)-(0.5*Phi-pi).*Phi.*cos(Phi)+2);
woQ=QCij.*Qi;