read("simplify.sum");
read("simplify.term");
x:=normal(x):
x1:=numer(x):
x2:=denom(x):
x2:=convert(taylor(x2,s=0,d+1),polynom):
x:=x1/x2;
readlib(coeftayl):
y:=coeftayl(x,s=0,d):
y:=normal(y):
HS:=normal(HS+y):
save HS,"simplify.sum":
quit;
