Integration

Again, the integration itself is performed by a single call to Mgfun[int_of_sys] :

> ode:=op(int_of_sys(sys,x=-infinity..infinity,takayama_algo));

However, the justification that the algorithm selected by the option applies to the integral under consideration is rather technical and is beyond this presentation.