I am trying to model a packed bed reactor with a system of three or four partial differential equations and Ordinary Differential Equations(momentum, energy and mass continuity equations). Does any know of any subroutines, either in the MKL library or eslewhere free which can do this for me ? Any solutions for C/C++ as well as for Fortran would be welcomed.
Regards