Add Phi_expr_val

1 files changed, 6 insertions, 0 deletions

diff --git a/series.mac b/series.mac
index 7e44e9b..6ee8655 100644
--- a/series.mac
+++ b/series.mac
@@ -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)));