This paper presents some FORTRAN subroutines that can perform mathematical operations on polynomial functions including addition, multiplication, differentiation and integration. The subroutines should be very useful to analysts who require high precision solutions to problems involving a large number of mathematical operations of polynomial functions.