用帕德形式的有理分式逼近已知函数

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