精选优质文档-倾情为你奉上首先说这个程序并不完善，为了实现通用（1，2，n）格式解题，以及为调用追赶法程序，没有针对节点数在三个以下的情况进行分类讨论。希望能有朋友给出更好的方法。首先，通过函数 sanwanj得到方程的系数矩阵，即追赶法方程的四个向量参数，接下来调用追赶法（在intersanwj函数中），得到三次样条分段函数系数因子，然后进行多项式合并得到分段函数的解析式，程序最后部分通过判断输入值的区间自动选择对应的分段函数并计算改点的值。附：追赶法程序 chase%function newv,w,newu,newd=sanwj(x,y,x0,y0,y1a,y1b)% 三弯矩样条插值% 将插值点分两次输入，x0 y0 单独输入% 边值条件a的二阶导数 y1a 和b的二阶导数 y1b，这里建议将y1a和y1b换成y2a和y2b，以便于和三转角代码相区别n=length(x);m=length(y);if m=nerror(x or y 输入有误，再来);endv=ones(n-1,1);u=ones(
