两点边值问题的差分逼近(共2页).doc

精选优质文档-倾情为你奉上实验名称：两点边值问题的差分逼近 实验目的与要求：掌握构造有限差分法的基本步骤会用有限差分法解两点边值问题实验内容：当a=0,b=1,p=0,q=r=1,时，构造逼近的中心差分格式。实验环境与器材：机房与MATLAB实验过程（步骤）或程序代码：clc;clear;n=input(输入n=);h=1/n;A=zeros(n-1,n-1);f=zeros(n-1,1);u=zeros(n-1,1);U=zeros(n+1,1);syms s;for i=0:n x(i+1)=i*h;endfor j=1:n-2 f(j,1)=(x(j)4-2*(x(j)3+13*(x(j)2-12*x(j)+2; A(j,j)=1-2/(h*h); A(j,j+1)=1/(h*h); A(j+1,j)=1/(h*h);endf(n-1,1)=(x(n-1)4-2*(x(n-1)3+13*(x(n-