diff options
| -rw-r--r-- | series.mac | 6 |
1 files changed, 6 insertions, 0 deletions
@@ -167,6 +167,12 @@ mul_series(factors, max_order) := block( γ_expr_val(x, n) := horner_exp_series(x*rhs(main_ansatz[1]), r, n, -n, 2); +Φ_expr_val(x, n) := block( + if x < 0 then ( + exp(x*Φc[0])*horner_exp_series(abs(x)*(Φc[0]-rhs(main_ansatz[5])), r, n, -n, 2)) + else ( + exp(x*Φc[0])*horner_exp_series(x*(rhs(main_ansatz[5])-Φc[0]), r, n, -n, 2))); + matchdeclare(m_any1, lambda ([e], freeof(β, γ, e))); matchdeclare(m_any2, lambda ([e], freeof(β, γ, e))); matchdeclare(m_any3, lambda ([e], freeof(β, γ, e))); |