function [bpx,bpy,wfxy] = grule2d (nquadx,nquady)
%usage:  [bpx,bpy,wfxy] = grule2d (nquadx,nquady);
[bpxv,wfxv]=grule(nquadx);
[bpyv,wfyv]=grule(nquady);
[bpx,bpy]=meshgrid(bpxv,bpyv);
[wfx,wfy]=meshgrid(wfxv,wfyv);
wfxy=wfx.*wfy;
