untitled.m

来自「五点差分型多重网格方法：各种插值算子的比较）」· M 代码 · 共 7 行

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);