Is there any option in Scipy to use "internally generated full Jacobian" with VODE/BDF method for stiff systems ? Any help would greatly be appreciated. -Lex