MATLAB-有限差分格式的验证

Matlab-有限差分格式的验证

提示：离散网格边界的向前向后二阶差分

---

一、代码

clear
clc
x=linspace(-1,1,1001)
y = linspace(-1,1,1001)

h1=0.002;
h2=0.002;
[y,x] = meshgrid(y,x);

syms xx yy z
z = sqrt(10^2-xx^2-yy^2)
fsurf(z)

digits 8

zz = zeros(size(x));
for i=1:101
    for j=