样条

d(i)=(y(i+1)-y(i))/h(i+1)-(y(i)-y(i-1))/h(i)

a(i)=(d(i)-c(i)*a(i+1))/b(i)

d(i)=d(i)*6/(h(i)+h(i+1))

k=j

d(0)=(6*(y(1)-y(0))/h(1)-z1)/h(1)d(n)=6*(z2-(y(n)-y(n-1))/h(n))/h(n)

yes

子程序yangtiao

c(i)=h(i+1)/(h(i)+h(i+1))a(i)=1-c(i)

c(0)=1a(n)=1

f=(x(k)-t(i))*(a(k-1)*(x(k)-t(i))*(x(k)-t(i))-h(k)+6*(k-1))f=f+(t(i)-x(k-1))*(a(k)*((t(i)-x(k-1))*(t(i)-x(k-1))-h(k)*h(k)))f=(f+(t(i)-x(k-1))*6*y(k))/(6*h(k))

f=y(j)

a(n)=d(n)/b(n)

a(i)=a(i)/b(i-1)b(i)=b(i)-a(i)*c(i-1)d(i)-a(i)*d(i-1)

return

x(j)＜t(i)

h(i)=x(i)-x(i-1)

x(j)≠t(i)