function?f?=?Pade（ynx0）
syms?t;

A?=?zeros（nn）;
q?=?zeros（n1）;
p?=?zeros（n+11）;
b?=?zeros（n1）;
yy?=?0;
a（1:2*n）?=?0.0;

for（i=1:2*n）
????yy?=?diff（sym（y）findsym（sym（y））n）;
????a（i）?=?subs（sym（yy）?findsym（sym（yy））?0.0）/factorial（i）;
end;

for（i=1:n）
????for（j=1:n）
????????A（ij）=a（i+j-1）;
????end;
????b（i1）?=?-a（n+i）;
end;

q?=?A\b;
p（1）?=?subs（sym（y）fin