interface(echo=0);
read fullparfrac;
f:=(2*x^7-7*x^5+26*x^3+8*x)/(x^8-5*x^6+6*x^4+4*x^2-8):
res:=fullparfrac(f,x):
normal(f-eval(subs(Sum=sum,res))):
if "=0 then print('ok') else print(`wrong result`,f,res) fi;
f:=x^3/(x^21+2*x^20+4*x^19+7*x^18+10*x^17+17*x^16+22*x^15+30*x^14+36*x^13+40
*x^12+47*x^11+46*x^10+49*x^9+43*x^8+38*x^7+32*x^6+23*x^5+19*x^4+10*x^
3+7*x^2+2*x+1):
res:=fullparfrac(f,x):
normal(f-eval(subs(Sum=sum,res))):
if "=0 then print('ok') else print(`wrong result`,f,res) fi;
f:=36/(x^5-2*x^4-2*x^3+4*x^2+x-2):
res:=fullparfrac(f,x):
normal(f-eval(subs(Sum=sum,res))):
if "=0 then print('ok') else print(`wrong result`,f,res) fi;
quit