Allow passing functions as parameters

Eg. for passing right hand sides to solver code. The alternative is to have the user write code that has direct access to the internal state, which I would rather avoid. @inducer, thoughts on the feasibility of this feature in Fortran?