- Fixed "unterminated sed command".

[ginac.git] / check / exam_pseries.cpp

diff --git a/check/exam_pseries.cpp b/check/exam_pseries.cpp
index a5625efa2c011afbc832c476a22e3069e3fd6e86..1857385f912a647d97c0bab918e0f66a6ad22df5 100644 (file)
--- a/check/exam_pseries.cpp
+++ b/check/exam_pseries.cpp
@@ -26,7 +26,7 @@ static symbol x("x");
 
 static unsigned check_series(const ex &e, const ex &point, const ex &d, int order = 8)
 {
-    ex es = e.series(x, point, order);
+    ex es = e.series(x==point, order);
     ex ep = ex_to_pseries(es).convert_to_poly();
     if (!(ep - d).is_zero()) {
         clog << "series expansion of " << e << " at " << point
@@ -107,7 +107,7 @@ static unsigned exam_series2(void)
     unsigned result = 0;
     ex e, d;
     
-    e = pow(sin(x), -1).series(x, 0, 8) + pow(sin(-x), -1).series(x, 0, 12);
+    e = pow(sin(x), -1).series(x==0, 8) + pow(sin(-x), -1).series(x==0, 12);
     d = Order(pow(x, 6));
     result += check_series(e, 0, d);
     
@@ -120,7 +120,7 @@ static unsigned exam_series3(void)
     unsigned result = 0;
     ex e, d;
     
-    e = sin(x).series(x, 0, 8) * pow(sin(x), -1).series(x, 0, 12);
+    e = sin(x).series(x==0, 8) * pow(sin(x), -1).series(x==0, 12);
     d = 1 + Order(pow(x, 7));
     result += check_series(e, 0, d);
     
@@ -153,32 +153,32 @@ static unsigned exam_series5(void)
     unsigned result = 0;
     ex e, d;
     
-    // Gamma(-1):
-    e = Gamma(2*x);
+    // tgamma(-1):
+    e = tgamma(2*x);
     d = pow(x+1,-1)*numeric(1,4) +
         pow(x+1,0)*(numeric(3,4) -
-                    numeric(1,2)*gamma) +
+                    numeric(1,2)*Euler) +
         pow(x+1,1)*(numeric(7,4) -
-                    numeric(3,2)*gamma +
-                    numeric(1,2)*pow(gamma,2) +
+                    numeric(3,2)*Euler +
+                    numeric(1,2)*pow(Euler,2) +
                     numeric(1,12)*pow(Pi,2)) +
         pow(x+1,2)*(numeric(15,4) -
-                    numeric(7,2)*gamma -
-                    numeric(1,3)*pow(gamma,3) +
+                    numeric(7,2)*Euler -
+                    numeric(1,3)*pow(Euler,3) +
                     numeric(1,4)*pow(Pi,2) +
-                    numeric(3,2)*pow(gamma,2) -
-                    numeric(1,6)*pow(Pi,2)*gamma -
+                    numeric(3,2)*pow(Euler,2) -
+                    numeric(1,6)*pow(Pi,2)*Euler -
                     numeric(2,3)*zeta(3)) +
-        pow(x+1,3)*(numeric(31,4) - pow(gamma,3) -
-                    numeric(15,2)*gamma +
-                    numeric(1,6)*pow(gamma,4) +
-                    numeric(7,2)*pow(gamma,2) +
+        pow(x+1,3)*(numeric(31,4) - pow(Euler,3) -
+                    numeric(15,2)*Euler +
+                    numeric(1,6)*pow(Euler,4) +
+                    numeric(7,2)*pow(Euler,2) +
                     numeric(7,12)*pow(Pi,2) -
-                    numeric(1,2)*pow(Pi,2)*gamma -
+                    numeric(1,2)*pow(Pi,2)*Euler -
                     numeric(2)*zeta(3) +
-                    numeric(1,6)*pow(gamma,2)*pow(Pi,2) +
+                    numeric(1,6)*pow(Euler,2)*pow(Pi,2) +
                     numeric(1,40)*pow(Pi,4) +
-                    numeric(4,3)*zeta(3)*gamma) +
+                    numeric(4,3)*zeta(3)*Euler) +
         Order(pow(x+1,4));
     result += check_series(e, -1, d, 4);