[qm,qn]=size(q); Q=spalloc(n,n,2*qm); for ii=1:qm, Q(q(ii,3),q(ii,4))=q(ii,5); Q(q(ii,4),q(ii,3))=q(ii,5); end