function?f?=?Atken（xyx0）
%x?y為坐標向量??x0為插值點的x坐標||?f為x0對應的值

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?=?subs（y1（n）‘t‘x0）;????????????%計算插值點的函數值
else
????simplify（y1（n））;???????????????????%化簡
????f?=?collect（y1（n））;????????????????%將插值多項式展開
????f?=?vpa（f6）;??????????????????????%將插值多項式的系數化成6位精度的小數
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-01-22?19:52??插值\
?????文件?????????734??2012-07-17?20:08??插值\Atken.m
?????文件????????1103??2012-07-17?20:30??插值\BSample.m
?????文件?????????449??2012-07-17?16:56??插值\forsee.m
?????文件?????????728??2012-07-17?20:13??插值\Hermite.m
?????文件?????????459??2012-07-17?20:59??插值\insert.m
?????文件?????????858??2012-07-17?20:06??插值\Language.m
?????文件?????????654??2012-07-17?20:09??插值\Newton.m
?????文件?????????722??2012-07-17?21:04??插值\resemble.m
?????文件?????????934??2008-04-09?18:12??插值\SecSample.m
?????文件?????????914??2012-07-17?20:16??插值\SubHermite.m
?????文件???????????0??2015-01-22?19:52??插值\使用說明.txt
?????文件?????????162??2012-07-18?16:40??插值\典型應用.txt