function?[S]?=?swanju（?xff1af1b?）
%?swanju??第一種邊界條件的三彎矩方法：
%???此處顯示詳細說明
%?x?為自變量x
%?f?為x對應的函數值f（xi）
%?f1a為f在左端點的一階導數值
%?f1b為f在右端點的一階導數值
syms?X;
n?=?length（x）;
h?=?diff（x）;
ff?=?diff（f）./h;????%ff表示三彎矩方法中的f（x（i），x（i+1））；
mu（11:n-2）?=?h（11:n-2）./（h（11:n-2）+h（12:n-1））;
lambda（12:n-1）?=?1-mu（11:n-2）;
mu（1n-1）?=?1;????????%μ最后一位為1?
lambda（11）?=?1;????%λ第一位為1
d?=?6./（h（11:n-2）+h（12:n-1））.*（ff（12:n-1）-ff（11:n-2））;
d0?=?6/h（1）*（ff（1）-f1a）;??????%計算d0
dn?=?6/h（n-1）*（f1b-ff（n-1））;????%計算dn
d?=?[d0ddn];
r（1）?=?2;
z（1）?=?d（1）;
for?i=2:n?????????????????%追
????l（i）?=?mu（i-1）/r（i-1）;
????r（i）?=?2-l（i）*lambda（i-1）;
????z（i）?=?d（i）-l（i）*z（i-1）;???
end
m（n）=z（n）/r（n）;
for?j=n-1:-1:1??????????%趕
????m（j）?=?（z（j）-lambda（j）*m（j+1））/r（i）;
end
for?j=2:1:n????????%講求得的m帶入書80頁式4.7.2，求得S（x）
????S（j-1）=vpa（（（x（j）-X）.^3）*m（j-1）/6/h（j-1）+（（X-x（j-1））.^3）*m（j）/6/h（j-1）+（f（j-1）-m（j-1）*（h（j-1）.^2）/6）*（x（j）-X）/h（j-1）+（f（j）-m（j）*（h（j-1）.^2）/6）*（X-x（j-1））/h（j-1）2）;
????fprintf（‘s（%d）=%1.2f‘j-1）;
????pretty（collect（S（j-1）））?%輸出方程形式
end
????
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1194??2019-12-11?09:07??三彎矩方法的三次樣條插值\swanju.m

?????文件????????382??2019-12-11?09:43??三彎矩方法的三次樣條插值\useswanju.m

?????目錄??????????0??2019-12-12?21:27??三彎矩方法的三次樣條插值

-----------?---------??----------?-----??----

?????????????????1576????????????????????3