FRODO: a MuPAD program to calculate matrix elements between contracted wavefunctions