## Integrate 1/(1+x^2), which we know is atan(x).
## It is less known that atan(x) = (1/2) * log((x+i)/(x-i)), which is
## the formula that is used in ratint. The exact value is pi/4.
ratint([1], [1, 0, 1], 0, 1) - pi/4

## The integral of a function with multiple poles: 1/((x+1) * (x+2)^3)
## The function convpow is ours (see the M-file)

ratint([1], conv([1,1], convpow([1, 2], 3)), 1, 2)


% ;;; Local Variables: ***
% ;;; mode: octave ***
% ;;; End: ***

