五点差分格式(共9页).doc

微分方程数值解大作业（一） 椭圆型方程编程计算：采用五点差分格式求如下椭圆型方程其中、及边条件为：1 ， 且边条件如下：问题存在精确解为： 2 ，且边条件如下：问题存在精确解为： 3 ，且边条件如下：问题存在精确解为： .代码：主函数1，差分解function g=fivepoints(x1,x2,y1,y2,M,N)%变步长法h=(x2-x1)/M; %横轴步长k=(y2-y21/N; %纵轴步长m=M-1;n=N-1;h1=h2;r=h1/k2; %五点中的上下两个点的系数t=2+2*r; %五点中的中心点的系数x=x1+(x2-x1)*(0:M)/M; %x，y向量表示横纵坐标y=y1+(y2-y1)*(0:N)/