MATLAB：龙贝格积分法微积分运算 有限差分方法求边值问题

龙贝格积分法微积分运算 有限差分方法求边值问题
在MATLAB中编程实现的龙贝格积分法的函数为：[I,step]=Roberg(f,a,b,eps)。
功能：龙贝格积分法求函数的数值积分。
调用格式：[I,step]=Roberg(f,a,b,eps)。
其中，I:积分值；
step: 积分划分的子区间次数；
f: 函数名；
a: 积分下限；
b: 积分上限；
eps: 积分精度；
有限差分方法求边值问题
在MATLAB中编程实现的有限差分方法的函数为：
[k,A,B1,X,Y,y,wucha,p]=yxcf(q1,q2,q3,a,b,alpha,beta,h)
功能：有限差分方法求函数的边值问题。
调用格式：[k,A,B1,X,Y,y,wucha,p]=yxcf(q1,q2,q3,a,b,alpha,beta,h)。
其中，q1、q2、q3：方程系数；
a为方程y的初始化值x0；
b为方程y的区间终值xn；
alpha为y(a)对应的值；
beta为y(b)对应的值；
h为求解步长；
k为计算迭代步数；
A为 矩阵系数；
B1为 矩阵系数；
Y为 矩阵向量；
y为求解目标方程的边值；
wucha为求解误差；
p等于[k’,X’,y,wucha’]，为一个矩阵