untitled.m

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

h=1/10;
[t,u] = meshgrid([h:h:(1-h)]);
UTRUE = sin(2*pi*t) .* sin(3*pi*u);
right= (h*h) * 13*pi*pi*UTRUE;   % right of equation 
disp(right);
p=reshape(right,81,1);
disp(p);