[R-sig-dyn-mod] jacobian for lsodes

Hi everybody,

does anybody have a working example for lsodes in combination with
compiled C code *including* jacvec? I have seen examples for lsoda where
the jacobian function returns the jacobian matrix. However, this does
not work for lsodes which requires jacvec instead.

Best regards,
Daniel