function pq = tensprod1(p,q)
pq = p*q';
pq = pq(:);
return