Partial fractions

Note: The EEPRO  8:Gain and Frequency ,Transfer function 
will work in all cases to produce partial fraction expansion
of given f(s)/g(s).

These examples below work fpr simple poles only
Expand() can give partial fractions but fails sometimea 
part(argument,n) changes value of argument with a sto command,

2 Variable Partial Fraction

:pf2()
:Prgm
:ClrIO
:DelVar x, a, b, abc
:(x+3)/(x^2+11*x+30)->y   example
:factor(y)->z
:getDenom(z)->d
:part(d2,2)->p2
:part(d2,1)->p1
:z-(a/p1+b/p2)->xx
:1.1->x    arbitrary
:xx->xx1
:2.3->x   arbitrary
:xx->xx2
:solve(xx1=0 and xx2=0,{a.b})->abc
:part(part(abc,1),2)->a
:part(abc,2)->b
:Delvar x
:a/p1+b/p2->pp
:Disp y
:Disp pp
:EndPrgm
==============
3 Variable Partial Fraction

:pf3()
:Prgm
:ClrIO
:DelVar x, a, b,c, abc
:(x-2)/(x^3+15*x^2+74*x+129)->y   example
:factor(y)->z
:getDenom(z)->d
:part(part(d,1),1)->p2
:part(part(d,1),2)->p1
:part(d,2)->p3
:z-(a/p1+b/p2+c/p3)->xx
:1.1->x    arbitrary
:xx->xx1
:2.3->x   arbitrary
:xx->xx2
:5.31->x  arbitrary
:xx->xx3
:solve(xx1=0 and xx2=0 and xx3=0,{a.b,c})->abc
:part(part(abc,1),2)->a
:part(part(abc,1),2)->b
:part(abc,2)->c
:Delvar x
:ClrIO
:a/p1+b/p2+c/pp3->pp
:Disp y
:Disp pp
:EndPrgm
======
4 Variable Partial Fraction
 
:pf4()
:Prgm
:ClrIO
:DelVar x, a, b, c, dd, abc
:x/(x^4+18*x^3+119*x^2+342*x+360)->y   example

:factor(y)->z
:getDenom(z)->d
:part(part(d,1),1)->k
:part(part(d,1),2)->p1
:part(k,2)->p2
:part(d,2)->p3
:part(k,1)->p4
:z-(a/p1+b/p2+c/p3+dd/p4)->xx
:1.1->x    arbitrary number
:xx->xx1
:2.31->x   arbitrary number
:xx->xx2
:5.31->x  arbitrary number
:xx->xx3
:5.555->x   arbitrary number
:xx->xx4
:solve(xx1=0 and xx2=0 and xx3=0 and xx4=0,{a.b,c,dd})->abc
  value of abc a=-.5 and b=-2.5 and c=1 and dd=2
 
:part(part(abc,1),2)->a    
     value of abc   b=-2.5 and c=1 and dd=2

:part(part(abc,1),2)->b
    value of abc c=1 and dd=2

:part(part(abc,1),2)->c
   value of abc  dd=2

:part(abc,2)->dd
   value of abc true

:Delvar x
:ClrIO
:a/p1+b/p2+c/p3+dd/p4->pp
:Disp y
:Disp pp
:EndPrgm