一维线性元有限体积法MATLAB程序 （1D linear finite volume element method）

• We seek the solution between $\Omega =(0,\pi )$;
• Model Problem:
  $$-\frac{d}{dx}(p\frac{du}{dx})+r\frac{du}{dx}+qu=f$$
  MATLAB code:

%一维一次有限体积法解两点边值精度阶测试
%修改版

iter=3;
pde.left=0;%左边界
pde.right=pi;%右边界
pde.Alpha=0;%左边值条件Dirichlet
pde.Beta=-1;%右边值条件Neumann
pde.p=@(x) x;
pde.r=@(x)  x;
pde.q=@(x) x;
pde.f=@(x) 2*x.*sin(x)-cos(x)+x.*cos(x);
pde.realu=@(x) sin(x);
pde.Drealu