Method of Solution

The mathematical rules for differentiation of sums, products, quotients, elementary functions in conjunction with the chain rule for compound functions are applied. The function should be expressed as one or more FORTRAN 90 or FORTRAN 77 procedures. A new type of variables is defined and the overloading mechanism of functions and operators provided by the FORTRAN 90 language is extensively used to implement the differentiation rules.