  A Matlab program is presented that computes derivative corrections in the
S-dual invariant formulation for IIB graviton scattering to any order in
perturbation theory. The coefficients of the four-point function are produced,
pertaining to the non-logarithmic terms. The program can be modified to find
coefficients of the higher-point functions. Instantons have not explicitly been
included.
