Added possibility to pass right hand sides containing function symbols.

For now, right hand side expressions are considered constants and do not enter the system matrix.
This, however, is an implementation detail and can be generalized easily later
through symbolic manipulations in the background.