function?a=ZJZXEC（xym）
if（length（x）?==?length（y））
????n?=?length（x）;??
else
????disp（‘x和y的維數不相等！‘）;
????return;
end??????????????????%維數檢查

syms?v;
d?=?zeros（1m+1）;
q?=?zeros（1m+1）;
alpha?=?zeros（1m+1）;
for?k=0:m
???px（k+1）=power（vk）;
end??????????????????????%x的冪多項式
B2?=?[1];
d（1）?=?n;
for?l=1:n
????q（1）?=?q（1）?+?y（l）;
????alpha（1）?=?alpha（1）?+?x（l）;
end
q（1）?=?q（1）/d（1）;
alpha（1）?=?alpha（1）/d（1）;
a（1）?=?q（1）;
B1?=?[-alpha（1）?1];
for?l=1:n
????d（2）?=?d（2）?+?（x（l）-alpha（1））^2;
????q（2）?=?q（2）?+?y（l）*（x（l）-alpha（1））;
????alpha（2）?=?alpha（2）?+?x（l）*（x（l）-alpha（1））^2;
end
q（2）?=?q（2）/d（2）;
alpha（2）?=?alpha（2）/d（2）;
a（1）?=?a