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