Note: please use Teguia (or Teguia Tabuguia) if you want to call me by my name. My first name is Bertrand.

I have a Ph.D. (Dr. rer. nat.) in Mathematics, in the field of Computer Algebra (2020), a master's degree in Mathematics (2018), and another in Computer Engineering (2016).

I am a Postdoctoral Researcher in the Non-Linear Algebra group of Bernd Sturmfels at the Max Planck Institute For Mathematics in the Sciences, Leipzig, Germany.

Research Interests: my research interests combine Symbolic Computation, Computer Programming, Differential Algebra, Algorithmic Number Theory, Algebraic Geometry, and the applications of all that in other scientific areas.

Check this website menu to learn more about my background.

Address: Inselstra├če 22, 04103 Leipzig, Germany

Office: F3 07


MPI MiS group member webpage:

UniKassel web page:

Tel:  +49 341 9959 768


Download my CV

Bertrand Teguia T.

Bertrand Teguia T.

'Magic' in symbolic computing

I invite you to read an informal presentation on computations based on Mathematics and computer programming in Computer Algebra.

1. Formal power series in Maxima 5.44: PDFHtml

2. Formal power series in Maple 2019: PDF

You can try some computations with my Maple library if you wish to:

  Download my implementations at The new FPS package.


Talk at the Maple conference 2022: Symbolic powers of functions defined by second-order linear ODEs.

Paper with Ait El Manssour, Rida and Sattelberger, Anna-Laura. D-algebraic functions. In preparation. Illustrations, software, and examples.

What is new in FPS (Maple version): HolonomicPDE,(not to be confused with HolonomicDE, which we want to keep for the univariate case only) for computing holonomic partial differential equations from a given D-finite or holonomic expression. There are two main options:

partialwrt: called as FPS:-HolonomicPDE(f, F(x1,...xn), partialwrt=xk);

Computes a holonomic partial differential equation with respect to a variable xk for a given multivariate D-finite (holonomic in a sense) function f in x1, ..., xk, ..., xn. When no option is specified, the output is a list of n partial differential equations that represent a canonical representation of the input expression as a holonomic function.

elimvars: called as FPS:-HolonomicPDE(f, F(x1,...xn), elimvars=S);

Computes a holonomic partial differential equation with the variables in the set S, a subset of {x1,...,xn}, eliminated from the polynomial coefficients. An accompanying paper with the technique used behind is under preparation.

One can use FPS:-HolonomicPDE to compute specific annihilators and perform so-called telescoping in the differential case.