untitled.asv

五点差分型多重网格方法：各种插值算子的比较）

nx=3;
ny=3;
N = nx*ny;
[XB,YB] = meshgrid([1:nx]/(nx+1),[1:ny]/(ny+1));
B       = XB+YB;
b       = reshape(B,N,1);
disp(b);
disp(length(b));