function?y?=?nbcspline（x0?y0?x）
n?=?length（x0）;???m?=?length（x）;
y?=?zeros（1?m）;
p（n）?=?0;??q（1）?=?0;??d（1）?=?0;??d（n）?=?0;
for?k?=?2:n-1
????h（1）?=?x0（k）?-?x0（k-1）;??h（2）?=?x0（k+1）?-?x0（k）;
????p（k）?=?h（1）?/?（h（1）?+?h（2））;
????q（k）?=?h（2）?/?（h（1）?+?h（2））;
????d（k）?=?6*（（y0（k+1）?-?y0（k））/h（2）?-?（y0（k）?-?y0（k-1））?/?h（1））?/?（h（1）?+?h?（2））;
end
b（1）?=?q（1）?/?2;?
for?k?=?2:n-1
????b（k）?=?q（k）?/?（2?-?p（k）*b（k-1））;
end

temp（1）?=?d（1）?/?2;
for?k?=?2:n
????temp（k）?=?（d（k）?-?p（k）*temp（k-1））?/