jlq = 1
if (jnode.eq.jq) jlq = 2
phival(inode,jnode) = Xquad(ilq)*Xquad(jlq)
INTEGER iq, jq, inode, jnode, i, j, bi, bj, ilq, jlq, m, n,Gi,Gj
INTEGER iq, jq, inode, jnode, i, j, bi, bj, ilq, jlq, m, n
uq = Xquad(ilq) * Xquad(jlq)
vq = Xquad(ilq) * Xquad(jlq)