fseipel@hpuxa.ircc.ohio-state.edu (Frank E. Seipel) (03/10/91)
Taylor series at a point. Calling sequence:
4: f(x)
3: x
2: order (integer)
1: point
Note: f(x) may be f(y), f(z), etc, but 3: must contain the name of this
variable.
Example:
4: 'X^3"
3: 'X'
2: 5
1: 3
result:
4:
3:
2:
1: 27+9*(-3+X)^2+(-3+X)^3+27*(-3+X)
STO as TAY or any other descriptive name.
Unrelated question: Is there a partial fractions program which can deal
with quadratic terms which can't be factored (ex: f(x)/(X^2+1) )? Wayne's
routines are great; however, they give the answer in terms of imaginary
numbers which is not desirable for my purposes.
------------------------cut here-----------------------
%%HP: T(3)A(R)F(.);
\<< 0 0 \-> eq v t p x
y
\<< v PURGE 0 t
FOR x eq
IF x
THEN 1 x
FOR y v \.d
NEXT
END p v STO
EVAL v PURGE v p -
x ^ x ! / *
IF x
THEN +
END
NEXT COLCT
\>>
\>>