Implementation of the recurrence relations of biorthogonality.General algorithms for computing some generalized Padé-type approximants are given. The quality of the used approximants depends on the choosen functionals for measuring the quality. In the classical example of the series f(b)=log(1+t)/t functionals are found which give better results than the usual ones. (netlib numeralgo na3)