优化后借助MATLAB用三次样条函数s(x)去模拟汽车门曲线--四川大学--硕士研究生数值分析作业

23.用三次样条函数去模拟汽车车门的曲线，车门曲线的型值点数据如下 (P89页) ：

`%%%
function f=cf(y)
 p=length(y);
 for k=1:p-1;
     f(k)=y(k+1)-y(k);
 end
%%%
 function [lk,mk,ek]=jh(x,h,f)
 n=length(x);
 for k=2:n-1;
     lk(k)=h(k)/(h(k-1)+h(k));
     mk(k)=h(k-1)/(h(k-1)+h(k));
     ek(k)=3*(lk(k)*f(k-1)+mk(k)*f(k));
 end
%%%
 function h=bc(x)
 n=length(x);
 for k=1:n-1;
     h(k)=x(k+1)-x(k);
 end
%%%
x=[0,1,2,3,4,5,6,7,8,9,10];
y=[2.51,3.30,4.04,4.70,5.22,5.54,5.78,5.40,5.57,5.70,