function?f?=?Atken（xyx0）
syms?t;
if（length（x）?==?length（y））
????n?=?length（x）;????
else
????disp（‘x和y的維數不相等！‘）;
????return;
end??????????????????????????????????????%檢錯

y1（1:n）?=?t;?????????????????????????????%符號函數數組要賦初值
for（i=1:n-1）???
????for（j=i+1:n）
????????y1（j）?=?y（j）*（t-x（i））/（x（j）-x（i））+y（i）*（t-x（j））/（x（i）-x（j））;????
????end
????y?=?y1;
????simplify（y1）;
end

if（nargin?==?3）
????f?=vpa（subs（y1（n）‘t‘x0）6）;????????????%計算插值點的函數值
else
????simplify（y1（n））;???????????????????%化簡
????f?=?collect（y1（n））;????????????????%將插值多項式展開
????f?=?vpa（f6）;??????????????????????%將插值多項式的系數化成6位精度的小數
end

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

?????文件????????687??2018-04-08?19:42??插值算法\Atken.m

?????文件???????1297??2018-04-09?21:42??插值算法\Atken1.m

?????文件??????30056??2018-04-08?17:10??插值算法\birdnest-32.tif

?????文件????????527??2018-03-26?20:19??插值算法\followup.m

?????文件???????2219??2018-03-17?20:18??插值算法\Gauss.m

?????文件???????1297??2018-04-09?21:43??插值算法\Gauss1.m

?????文件???????5174??2018-03-22?17:02??插值算法\IM380.bmp

?????文件????????446??2018-04-09?18:09??插值算法\Language.m

?????文件???????1306??2018-04-09?21:43??插值算法\Language1.m

?????文件????????271??2018-03-26?10:04??插值算法\lineint.m

?????文件???????1309??2018-04-09?21:40??插值算法\lineint1.m

?????文件????????430??2018-03-31?09:16??插值算法\m.m

?????文件????????607??2018-03-17?20:12??插值算法\Newton.m

?????文件???????1316??2018-04-10?10:12??插值算法\Newton1.m

?????文件???????1603??2018-03-31?09:59??插值算法\sancichazhi.m

?????文件???????1336??2018-04-09?21:47??插值算法\sancichazhi1.m

?????文件???????1475??2018-04-10?10:25??插值算法\ThrSample1.m

?????文件???????1374??2018-04-10?10:10??插值算法\ThrSample11.m

?????文件???????1536??2018-03-31?15:28??插值算法\ThrSample2.m

?????文件???????1358??2018-04-09?21:46??插值算法\ThrSample21.m

?????文件??????14039??2018-04-09?22:04??插值算法\untitled.fig

?????文件??????10040??2018-04-10?10:53??插值算法\untitled.m

?????目錄??????????0??2018-04-10?22:30??插值算法

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

????????????????79703????????????????????23