诺伊曼边界条件，帮我改一下程序，我礼拜五有用，用的matlab，

Question

各位网友好，我用了狄利克雷边界条件可以做，但是不知道怎么用诺伊曼边界条件，可以帮我改一下程序吗，我礼拜五有用，
用的matlab，

%使用示例:w=burgers (0,1,0,2,20,40)
function w=burgers(xl,xr,tb,te,M,N)
alf=5;bet=4;D=1.0;
xl=0;
xr=1;
tb=0;
te=2;
M=20;
N=40;
f=@(x)2*D*bet*pi*sin(pi*x)./(alf+bet*cos(pi*x));
l=@(t)0*t;
r=@(t)0*t;
h=(xr-xl)/M; k=(te-tb)/N; m=M+1; n=N;
sigma=D*k/(h*h) ;
w(: ,1)=f(xl+(0:M)*h)';
w1=w;
for j=1:n
  for it=1 :3
%牛顿迭代
    DF1=zeros(m,m);DF2=zeros (m, m) ;
    DF1=diag(1-k+2*sigma*ones(m, 1))+diag(-sigma*ones(m-1,1),1);
DF1=DF1+diag ( -sigma*ones (m-1,1),-1);
DF2=diag (2*k*w1);
DF=DF1+DF2;
F=-w(:,j)+(DF1+DF2/2)*w1 ;
DF(1, :)=[- 3 4 -1 zeros(1,m-3) ] ;F(1)=DF(1,: )*w1;
DF(m, : )= [zeros(1,m-3) -1 4 -3];F(m)=DF (m,:) *w1;
%对于DF的狄利克雷边界条件
F(1)=w1(1)-l(j);F(m)=w1(m)-r(j);%对于F的狄利克雷边界条件
w1=w1-DF\F;
  end  
  w(:,j+1)=w1;
end
x=xl+(0:M)*h; t=tb+ (0:n)*k;
mesh(x,t,w')
%解w的三维图    
end

Answer 1

参考一下呢

利用matlab的PDE工具箱求解Neumann边界的Poisson方程 - liang_l - 博客园 我们求解的问题是上述问题的精确解为 u=sin(3*pi*x+pi/4)*sin(2*pi*y+pi/4);接下来用Matlab工具箱求解首先打开matlab，命令窗口处输入：pdetool将看到如下 https://www.cnblogs.com/kmliang/archive/2012/12/28/2836658.html

边界条件的matlab式子,matlab程序求解含有边界条件的线性方程组_吕鑫明的博客-CSDN博客 编写程序，求解含有边界条件(本质)的线性方程组一、编写程序思路：拟编写一MATLAB 程序函数，来处理边界条件(a x k =), 将原线性方程组转化为形如：????????????????----=????????????????????????????????a A b a A b a A b a a A b x x x x x A A A A A A A A A A A A A A A A ... https://blog.csdn.net/weixin_36001196/article/details/115831897