function[EF]?=?diophantine（ANd）
%?Computes?polynomials?E（z^-1）?e?F（z^-1）?
%?delta?=?1-z^（-1）
delta?=?[1?-1];
%?A?=?1?+?a1?z^（-1）?+?...?+?a2?z（-na）
AD?=?conv（Adelta）;
%?note?that?nAD?=?n~a?+?1
nAD?=?size（AD）;
nAD?=?nAD（2）;
%?compute?horizons
?N1?=?d?+1;
?N2?=?d?+?N;
%?Compute?F（z^-1）
%?inilialization?vector?f
f（1:）=?[1?zeros（1nAD-2）];
%?i?=?0?...?n?-1
for?j?=?1:?N2;
%?Note?that?for?i?=?1?corresponds?to?f（j0）
for?i?=?1:nAD-2
???f（j+1i）?=?f（ji+1）-f（j1）*AD（i+1）;
end
???f（j+1nAD-1）?=?-f（j1）*AD（nAD）;
end
F?=?f（1+N1:1+N2:）;
%?Computes?E（z^-1）
E?=?zeros（N2）;
e（1）?=?1;???????????????%?for?the?special?case??1/~A
E（11）?=?e（1）;
for?i?=?2:?N2
????e（i）?=?f（i1）;
????E（i1:i）=e;
end
E?=?E（N1:N2:）;

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

?????文件????????740??2009-09-25?11:05??gpc\diophantine.m

?????文件????????450??2009-09-24?09:58??gpc\getG.m

?????文件????????174??2009-09-25?10:47??gpc\getGij.m

?????文件???????1959??2009-09-26?07:59??gpc\gpc.m

?????文件?????????85??2009-09-25?11:10??gpc\rankn.m

?????文件????????836??2009-09-24?09:57??gpc\staten.m

?????文件????????577??2009-09-24?09:57??gpc\tf2AB.m

?????目錄??????????0??2009-09-26?07:58??gpc

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

?????????????????4821????????????????????8