In reply to Daniel Vazquez.

Yes. As I state in the first sentence of the first section, you can use the ODE subroutine to solve systems of differential equations in SAS 9. The ODE subroutine uses the same user-defined function to define a system (such as the SIR model), but has a slightly different interface (API). Sere the SAS 9 documentation for the ODE Call.